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Abstract 

We give a detailed account of various connections between several classes of objects: Hankel, 
Hurwitz, Toeplitz, Vandermonde and other structured matrices, Stietjes and Jacobi-type continued 
fractions, Cauchy indices, moment problems, total positivity, and root localization of univariate poly- 
nomials. Along with a survey of many classical facts, we provide a number of new results. 
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Introduction 

This survey and research paper offers a gUmpse at several classical topics going back to Descartes, Gauss, 
Stieltjes, Hermite, Hurwitz and Sylvester (see, e.g., [HI [80l [281 [811 [82] ) , sl\ connected by the idea that 
behavior of polynomials can be analyzed via algebraic constructs involving their coefficients. Thus a 
number of inter-related algebraic constructs was built, including Hurwitz, Toeplitz and Hankel matrices, the 
corresponding quadratic forms, and the corresponding continued fractions. The linear-algebraic properties 
of these objects were shown to be intimately related to root localization of polynomials such as stability, 
whereby the zeros of a polynomial avoid a specific half-plane, or hyperbolicity, whereby the zeros lie 
on a specific line. These methods gave rise to several well-known tests of root localization, such as the 
Routh-Hurwitz algorithm or the Lienard-Chipart test. 

In the 20th century, this line of research was developed further by Schur [24] , Polya [78] , Krein [59l [60l 
[34] and others (see, e.g., [43]), leading to important notions of total positivity, Polya frequency sequences, 
stability preservers etc. A very important part of that research effort was devoted to entire functions, 
in particular, the Laguerre-Polya and related classes. However, this area gradually went out of fashion 
and was essentially abandoned around 1970-1980s, with a few exceptions, such as Pade approximation. 
We should note in this connection that the closely related theory of continued fractions, initiated by 
Chebyshev, Stieltjes and Markov and further developed by Akhiezer, Krein and their collaborators in 
connection with problems of mechanics is only now returning to the forefront due to its connections with 
orthogonal polynomials [SH| [BU [SH [ISl [SO] ■ 

On the other hand, various fast algorithms for structured matrices (Cauchy, Vandermonde, Toeplitz, 
Hankel matrices, and more generally quasiseparable matrices and matrices with displacement structure) 
have been developing rapidly in the last few decades due to the efforts of Gohberg, Olshevsky, Eidelman, 
Kailath, Heinig, Sakhnovich, Lerer, Rodman, Fuhrmann and others, along with applications to control 
theory and engineeting: see the collections [TO] [TTJ [72j [53j [TS] [22] and the references therein. Surveying 
these developments is, however, outside the scope of this paper. 

Furthermore, we must stress the fact that the results and algorithms in this paper are not approximate 
but exact (i.e., give exact answers in exact arithmetic) as is the term localization, which is often used 
in the sequel. For instance, we perform root localization with respect to a line or a half plane and are 
able to determine exactly, e.g., how many zeros of a given polynomial lie on a given line. For sure, this 
does not obviate the need for error analysis of the corresponding algorithms that use finite precision or 
floating-point arithmetic. 

We must note that the questions of root localization are returning into the spotlight also due to the 
development of basic algebraic techniques for multivariate polynomials [lU [40] [TH [iTl [161 [83] and entire 
functions [9, 20 , as well as due to newly found applications of results on multivariate stable and hyperbolic 
polynomials to other branches of mathematics [42l [41] I100[ [T6] . Very promising applications of current 
interest include problems of discrete probability, combinatorics, and statistical physics, such as the analysis 
of partition functions arising in classical Ising and Potts models of statistical mechanics j92i i89[ . 

The goal of this paper is to provide a comprehensive and coherent treatment of classical connections 
between three kinds of objects: various structured matrices, representations of rational functions via 
continued fractions, and root localization of univariate polynomials. Our additional goal is to demonstrate 
that, despite the rich history of this subject, even the univariate case is far from being exhausted, and 
that classical algebraic techniques are useful in answering questions about the behavior of polynomials and 
rational functions. 

Our future goals include using this work for generalizations of these classical results to univariate entire 
and meromorphic functions, including questions on Polya frequency sequences, Hurwitz rational functions, 
and entire functions with all real zeros [98, 12 . One particular area of interest is the class of so-called 
generalized Hurwitz polynomials, which is a useful large class containing all Hurwitz stable polynomials. 
This line of our ongoing research is closely related to the theory of indefinite metric spaces, indefinite 
moment problems, and the corresponding eigenvalue problems [23l [76l [771 [47] . 

We now illustrate our main point above by discussing several new results in this paper. 

Our first illustration is provided by the body of work in Section 14.21 These results provide explicit 
criteria for a polynomial to have only real roots in terms of Hankel and Hurwitz determinants made of its 
coefficients. To give the reader an idea of these statements, we quote two sample theorems from Section r4.2l 
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Consider a polynomial 

p{z) — uqz" + aiz"^^ + • • • + a„, ai, . . . , a„ G M, oq > 0. 

and let 

L{z):-- 



p'{z) So Si S2 



p{z) z z^ z^ 
be its logarithmic derivative. Define the associated infinite Hankel matrix 



S:=S(L):=\\s 



oo 

j+j ll ij =0 



and its leading principal minors 



So 
si 



Sl 
S2 



S2 
S3 

Sj + 1 



S2J-2 



J = 1,2,3,. 



Here is a sample result from Scction [4.2l that provides an interesting connection between total positivity 
of a special Hurwitz matrix and polynomials all whose roots are real and negative: 

Theorem (Total positivity criterion for negative zeros). The polynomial 



g(z) = Oo + a\z + • • • + a„z" 
has all negative zeros if and only if the infinite matrix 
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is totally nonnegative. 

This result is based on a new connection between Stieltjes continued fractions representing a rational 
function P = q/p and a special factorization of the infinite Hurwitz matrix Hij), q) associated to the pair 
(p, q) (see Section [3^ . This connection implies a criterion of total nonnegativity of the latter infinite Hur- 
witz matrix H{p, q) fTheorem l3.44p . which in turn implies our criterion for negative zeros of a polynomial 
g via the total nonnegativity of the infinite discriminant matrix I?oo(ff)- One direction of Theorem 13.441 
was essentially developed in the earlier works of Asner and Kempermanpl [57] . while the other direction is 
given here for the first time. This result simultaneously provides a criterion of stability [51 [571 [551 [55] . 

Section [4.21 provides many other alternative criteria for real zeros, including the cases of only positive, 
only negative, or mixed zeros, in terms of Hankel and Hurwitz determinants as well as the coefficients of 
Stieltjes continued fraction expansion of the logarithmic derivative of a given polynomial. The special role 
played by the logarithmic derivative for the analysis of real zeros is clarified earlier in Section 21 together 
with the main ideas behind counting real roots to the left and to the right of the origin. The underlying 
theory of Cauchy indices is presented in Section 13.31 

Another important class of objects behind the results of this paper on polynomial roots is the class 
of so-called i?-functions. These are rational functions mapping the upper-half plane of the complex plane 
either to the lower half-plane or to itself (see Definition [OJ- The theory of i?- functions connects several key 
objects of this work: continued fractions, Cauchy indices, polynomials with real roots, the moment problem, 
and determinantal inequalities. As an example, we quote the following new result from Section [3.31 



Theorem (Generalized Lienard-Chipart criterion). If a real rational function R = q/p, withp and q 
relatively prime, maps the upper half-plane of the complex plane to the lower half-plane, then the number of 
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its positive poles is equal to the number V~(ao, . . . , a„) of strong sign changes in the sequence (oq, . . . , a„) 
of coefficients of its denominator p. In particular, R has only negative poles if and only if aj > for 
j = 1,2, ... ,n, and only positive poles if and only if aj^iOj < for j ~ 1,2, ... ,n. 

Our final illustration is the generalized Orlando's formula from Section 11.21 Orlando's formula per se 
expresses the Hurwitz determinant of order n — 1 associated to a polynomial p of degree n as the product, 
up to a certain normalization, of all possible pairs Zi + Zj of the zeros of p. Orlando's result goes back to 
1911 [75]. The generalized Orlando formula obtained in this paper says the following: 

Theorem (Generalized Orlando formula). Let p be a polynomial of degree n and q a polynomial of 
degree m < n. Then the resultant of these polynomials can be computed as follows: 

R(p,g) = (-l)^a™+" n (^'' + ^^-) 

l<i<j<2n 

where ao is the leading coefficient of p and where zi, i — \, . . . , 2n, are the zeros of the polynomial 

h{z) :—p{z'^) + zq(z^). 

The classical Orlando's formula then follows as a special case by splitting an arbitrary polynomial into its 
even and odd part and applying the generalized Orlando's formula. 

As these examples show, new connections can be found between several classical matrix classes (Hur- 
witz, Hankel, Vandermonde etc.) made of the coefficients of polynomials, different representations of 
rational functions (Stieltjes and Jacobi continued fractions, Laurent series, elementary fractions), and var- 
ious counting notions (sign changes, the number of roots in a specific domain of a complex plain, Cauchy 
indices). Underlying all these is a coherent theory demonstrating how these connections arise. 

The presented theory is completely general in that it does not single out the special case of stable 
polynomials. The importance of stability, historically the first question on root localization pTl[CTll5^[5^ 
157] , is by now quite well understood. Various stability criteria have been studied in great detail, especially 
in the engineering literature. In this survey (Section 14]), we choose instead to illustrate applications of the 
general theory to polynomials with real roots, a less studied but arguably equally important special case. 

To summarize, the point of this work is to provide a uniform, streamlined, algebraic treatment for a 
body of questions centered around the root localization of polynomials. At present, results of this type 
are scattered in the literature, whereas even textbooks and monographs about polynomials sometimes 
fail to provide answers to some basic questions (i.e., tests for real roots using only the coefficients of the 
polynomial). We hope that this work will serve as a useful reference and source of further research for 
mathematicians in various fields interested in polynomials and their zeros. 



1 Complex rational functions and related topics 

Consider a rational function 

p{z) 



z^R{z):^^ (1.1) 



where p and q are polynomials with complex coefficients 

p{z) :— oqz"- + oiz"^^ + • • • + a„, ao, fli, . . . , a„ e C, oq 7^ 0, (1.2) 

g(^) := 6o2" + 6iz"-i + ---+6„, 60, e C, (1.3) 

If the greatest common divisor of p and q has degree I, < I < n, then the rational function R has exactly 
r = n — l poles. Note that zeros or poles of a rational function are always counted with multiplicities unless 
explicitly stipulated otherwise. Thus, in the rest of the paper, the phrase "counting multiplicities" will be 
implicit in every statement about zeros or poles of functions under consideration. 



1.1 Hankel and Hurwitz matrices. Hurwitz formulae 

In this section, we introduce Hankel and Hurwitz matrices associated to a given rational functions and 
discuss the Hurwitz formulae that connect these two classes. 
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Expand the function (jl.ip into its Laurent series at oo: 



Riz) 



Here Sj — for j < n — \ — m and s„_i_ 



S2 



s_i H h ^ + ^ 



(1.4) 



a-o ' 



where m — deg 



The sequence of coefficients of negative powers of z 

So, Sl, S2, ■ • ■ 

defines the infinite Hankel matrix S :—S{R) := ||si+j ||°°^o- This gives the correspondence 

SiR). 

Definition 1.1. For a given infinite sequence {sj)°^Q, consider the determinants 



(1.5) 



DjiS):-- 



so 

Sl 
S.7-1 



51 S2 

52 S3 

Sj Sj+l 



Si 



S2i-2 



J = 1,2,3,. 



(1.6) 



i.e., the leading principal minors of the infinite Hankel matrix S. These determinants are referred to as 
the Hankel minors or Hankel determinants. 

An infinite matrix is said to have finite rank r if all its minors of order greater than r are zero whereas 
there exists at least one nonzero minor of order r. Knonecker |63j proved that, for any infinite Hankel 
matrix, any minor of order r where r is the rank of the matrix, is a multiple of its leading principal minor 
of order r. This implies the following result. 

Theorem 1.2 (Kronecker [63]). An infinite Hankel matrix S — \\si+j\\°°^i has finite rank r if and only if 



DriS) ^ 0, 

Dj{S) = 0, for all j>r. 



(1.7) 
(1.8) 



Let Dj{S) denote the following determinants 



51 S2 

52 S3 



S3 
S4 



Sj + 1 Sj+2 



S2J-1 



j = l,2,3,.. 



(1.9) 



With a slight abuse of notation, we will also write Dj{R) := Dj{S{R)) and Dj{R) Dj{S{R)) if the matrix 
S ~ S{R) is made of the coefficients (|1.4p of the function R. 

The following theorem was also established by Gantmacher in |36) . 



Theorem 1.3. An infinite matrix S ■ 



^j^i has finite rank if and only if the sum of the series 



R{z) = 1" ^ + ^ 



is a rational function of z. In this case the rank of the matrix S is equal to the number of poles of 
the function R. 



Theorems 11.21 and 11.31 have a simple corollary, which will be useful later. 

Corollary 1.4. A rational function R with exactly r poles represented by the series (|1.4|) has a pole at the 
point if and only if 

Dr^i{R)^0 and Dj{R)=0 for j = r, r + 1, . . . . (1.10) 
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Proof. Since the function R represented as a series (11.41) has exactly r poles, the function 

G{z) := zR{z) — zs_i = sq ■ 



z 



£2 

r2 



£3 

r3 



has exactly r — 1 poles if and only if R has a pole at the point 0. If i? does not have a pole at 0, then G 
has r poles. Thus, by Theorems 11.21 and 11.31 the function R has a pole at if and only if 



Dr^i(G) 7^ and Dj{G) ==0, j = r, r + 1, . 



Since 



D,{R) = D,{G), j = l,2, 
the formula ()l.lip yields the assertion of the corollary. 



(1.11) 



□ 



Theorems 11.21 and 11.31 imply the following: if the greatest common divisor of the polynomials p and q 
defined in (|1.2I) - (|1.3I) has degree ^ < / < degq, then the formula (|1.7I) - (|1.8I) hold for r = n — Z for the 
rational function ()1.1|) . where n is the degree of the polynomial \\.2\ . i.e., the denominator of R. 

Denote by W2j{p, q) the following determinants of order 2j: 



ao 


ai 


02 . 


■ aj-i 


aj . 


■ a2j-i 


bo 


bi 


b2 ■ 


■ bj^i 


bj . 


■ b2j-i 





ao 


ai . 


■ 0.0-2 


aj^i . 


■ 0.23-2 





bo 


bi . 


■ bj^2 


bj-i ■ 


■ b2j-2 








. 


ao 


ai 


aj 








. 


■ bo 


bi . 


■ b, 



i = i,2,. 



(1.12) 



constructed using the coefficients of the polynomials (jl.2p - (|1.3p . Here we set ai — for all i > n and 
bi — Q for / > m. The determinants V2j(p, <z) are called determinants of Hurwitz type or just Hurwitz 
determinants or Hurwitz minors. 

In his celebrated work [48] , A. Hurwitz found relationships between the minors Dj {R) defined in ()1.6p 
and the determinants W2j{p,q) defined by (jl.l2p . 



Theorem 1.5 ( [48l [611 [36l E] ) . Let Riz 

Then 



p{z) 

S/2j{p,q)^al'D,{R), J = 1,2,.... 
Proof. From (|1.4p it follows that coefficients Sk,ak,bk satisfy recurrence relations 



, where the polynomials p and q are defined by p.2p - (|1.3p . 

(1.13) 



a,s_i. 



J =0,1, 2,. 



(1.14) 



These recurrence relations imply the formulae (jl.l3p by direct matrix multiplication, once we take into 
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account (fTTi)) (cf. ^) 
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s_i So 
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Sj-2 Sj-1 



Sj-3 Sj_2 


• ■ • S2J-2 



• ■ • S2J-3 






"2j-l 


. . . ao 









s_i So 


■ ■ ■ Sj-i 



















1 






ao 


ai 


a2 ■ 


■ "i-i 


aj . 


a2j 





bo 


bi . 


■ bj-2 


bj-i ■ 


■ b2j-i 





ao 


ai . 


■ aj-2 
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■ a2j-i 








bo ■ 


■ b^-3 


bj-2 ■ 


■ b2j-2 








. 


■ bo 


bi . 


• b, 








. 


ao 


ai 


aj 



aoV2j b, g). 



Remark 1.6. Since the formulae ()1.13p are algebraic identities also valid for polynomials p = X] j=o '"i ' 
ao 7^ 0, and q — X]J=o^J^"'' ^^^y ^^^^ hold when p and q are entire functions and i? is a meromorphic 
function. 



From p.l2p one can see that, whenever 6o 7^ 0, we have 

V2,(p,g) = V2,(-g,p), J = 1,2,.... (1.15) 

This observation, coupled with Theorem II. 5[ yields the following well-known fact (cf. 27^. 

Corollary 1.7. For two infinite sequences S ■.={sj)J!L_i, T ■.—{tj)°^_i with s_i 7^ and i_i 7^ 0, the 
following conditions are equivalent: 

1) s_it_i + l = 0, 

Sit-i + Si-ito -\ 1- so<j-i + s_i<, = 0, i = 0,l,2, 

2) Dk{S) ^ s^_^\Dk{T), k = l,2,..., 

where Dk{S) and Dk(T) are the Hankel minors (11.61) for the sequences S and T , respectively. 
Proof. Let 



R{z) = 44 = s_i + — + — + — 



then 



p{z) 
p{z) 



So , S]^ , S2 
2 ' ^3 



R{z) 



q{z) 



to tl t2 



^It appears that this basic fact was known much earher than the work of Edrei but we do not know of the original reference. 
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The assertion of the corollary follows from Theorem II. 5 [ formula (|1.15p and the identity 



' So Sl S2 

s_i H h — + — 



t2 



t-i + - + ^ + ^ 



1. 



1.2 Resultant formulae and applications: discriminant and Orlando formulae 

This section is devoted to the resultant, an important object that is closely related to the greatest common 
divisor of two polynomials. The discriminant of a single polynomial occurs as a special case of the resultant 
of that polynomial and its derivative. Moreover, the generalized Orlando formula connects the resultant 
of two polynomials, p and q, with the roots of the aggregate polynomial h{z) :=p(z^) + zq{z'^). Although 
the resultant is classically known, this last connection is new. 

Definition 1.8. Let two polynomials p and q be given in (|1.2|) - (jl.3|) with n:=degp and m:=degq. The 
resultant of p and q is the following determinant of order n + m 



R{p,q): 



ao 


ai 


^m— 1 




O-n-l 


an 


^n+m — l 





ao 






O-n-2 




O-n+m-2 








ao 


fli 


^n — rn 


^n — m+1 


an 


n — 711 


bn-m+l 


■ bn^i 


bn 


b2n-rn-l 


b2n-m 


&2n-l 





bn—m 


■ bn-2 


bn-1 ■ 


b2n-7n-2 


b2n-rn-l 


&2n-2 














bn—m 


bn—m-\-l 


bn 



(1.16) 



where we set : = and bi:—0 for all i > n. 
From (|1.16p it can be immediately seen that 

R(p,g) = (-l)"'"R(g,p). 

The resultant of two polynomials is a multi-affine function of the roots of these polynomials, as the 
following formula shows. 

Theorem 1.9. Given polynomials p and q as in (|1.2p - (jl.3p with bo ^ 0, let Xi (i — 1, . . . ,n) denote the 
zeros of p, and let fij (j = 1, . . . , n) denote the zeros of q. Then 

n(n-l) " " " 

R(p,g) = (-l)^^V2„(p,<z) = aZX{q(\^i) = a^obl J] " = {-^TKX{p{n)- (1-17) 

i=l ij = l i=l 

Proof. At first, assume that all roots of the polynomial p are distinct. Then the function R admits the 
representation 

n 

a.; 



R(,z)^P + Y. 



z — Xi 



where 



, j = l,...,n. 



(1.18) 
(1.19) 



Here \j ^ Xi whenever i j, according to the assumption. Comparing the representation (jl.lSp with the 
expansion ()1.4p . we obtain 



Sk 



^a,Af, fc = 0,l,2,. 



(1.20) 
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bo 



and .s_i = — ^ 0. From (|LT^ . (ITI^ and it follows that 

ao 



i{n-l) 



^{n-l) 



R(p,g) = (-l)^^V2„b,g) = (-l)^^af i?„(i?). 



The formulae (IL^ yield 

/ So Sl S2 
Sl S2 S3 

\Sn— 1 Sji-^i 



Sn 
S2n-2/ 



This formula combined with the residue formula (I1.19P implies 
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we have 
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anX, 
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1 A3 
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02^2^^ 


"3^3"^ 


■ ■ anX"^' 
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X2 

A„ . . 


Arv 


This equality implies 
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\ n — 1 
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\n— : 
A„ 



i=l ^' j<z 



p'{z) = ao^W_{z- Afe), 
1=1 fe=i 

n n n 



'i=i fc=i 



j<i i<j 

This product has exactly n{n — 1) factors of the form Ai — Xj, therefore 



]Jp'iX,) = aj}(-l)^^ n(A, - A, 



Now from pT^ - (ITT^ we obtain 



n(n — l) 



(1.21) 



(1.22) 



(1.23) 



(1.24) 



(1.25) 



The formula pTTT)) follows from ([L^ and PT^ . 

If the polynomial p has multiple zeros, we can consider an approximating polynomial p^ with simple 
zeros such that 

Iimpj(z)=p(z) for all z. 

£— >-0 
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Then 



where Re{z) 



V2„(i?e) > ^2n{R), 



. The formula (|1.17l) is vahd for the polynomials q and , so it is also valid at the 
limit, i.e., for the polynomials q and p. □ 

Corollary 1.10. Under the conditions of Theorem ] 1 . 91 let degq = m < n. Then 

R{p,q) = (-l)^^a^-"V2„(p,q) = a",^f[q{X,) = 

n m 7n (1.26) 

i=lj=l j = l 

Proof. If degg = rn{< n), then the first nonzero coeSicient of q is bn-m, according to (|1.3p . Therefore, 



n(n— 1) 
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CLji—fYi—i 
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a2n-l 
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^n — m— 1 
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bn 



Thus, we have 



n(n— 1) 

V2„(p,g) = (-l)^^a^™R(p,g). 



This relation and the formula ()1.13|) . ()1.25|) yield ()1.26p when p has only simple roots. But (|1.26p is also 
valid when p has multiple zeros, which can be proved by an approximation argument as above. □ 



The formulas ()1.17p and (jl.26p now imply the well-known property of the resultant: 
Corollary 1.11. R(p, q) — if and only if the polynomials p and q have common roots. 

Next, we consider a function that allows us to test whether a single polynomial has multiple roots. 
Definition 1.12. Given a polynomial (|1.2p with roots (i — 1,. . . ,ri), its discriminant is defined as 



T>{p):=al-'l[{K~X,f 



(1.27) 



It is clear from (|1.27l) that the discriminant of a polynomial is equal to zero if and only if the polynomial 
has multiple zeros. But multiple zeros of a polynomial are the zeros that it shares with its derivative. The 
following connection between the discriminant of p and the resultant of p and p' should not come as a 
surprise. 



Theorem 1.13. Given a polynomial (|1.2p of degree n, we have 



R(p,p') = (-l)^^aoD(p). 



(1.28) 
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Proof. Indeed, the resultant R(p,p') satisfies (|1.26l) . Together with (|1.24p and (jl.27p . it gives 

" ri(n-l) ri(n-l) 

2—1 j<i 



From (|1.26p we obtain 



n(n— 1) 

V2n(p,p') - (-l)^^anR(p,p') = agD(p). 



Thus, the discriminant of the polynomial p is the following (2n — 1) x (2n — 1) determinant: 



D(p) 



1 

ao 



nao 


{n — l)ai 


{n 


- 2)02 . 


a-n-i 








ao 


ai 




a2 


a-n-i 


a„ 








nao 


{n 




2a„_2 










ao 




ai 


a„_2 















nao 


3a„_3 


2an-2 













ao 


a„_3 


an-2 
















ai 


02 


an 















(ri — l)ai 





(1.29) 



We now give one more formula for the resultant, which is very close to the well-known Orlando for- 
mula [73] (see also [36]). More precisely, the application of this formula to the resultant of the odd and 
the even parts of a polynomial yields exactly the Orlando formula. 

Theorem 1.14. Let the polynomials p and q be given by (jl.2p - (11.31) . with degp — n and deg q = m < n—1. 
Then the resultant of these polynomials can be computed as follows: 



%{n-l) 



R(p,g) = (-l)^^a^'+" n + 



l<i<k<2n 

where Zi (j = 1, . . . , 2n) are the zeros of the following polynomial of degree 2n 

h{z)^p{z')+zq{z'). 



(1.30) 



(1.31) 



Proof. Let Ai (i = 1, . . . , n) be the zeros of the polynomial p, and let fij {j = 1, . . . , m) be the zeros of 
the polynomial q. From ()1.3ip it follows that 



and 



From (I1.32P we obtain 

pi4) 



h{z) + h{-z) 



j = 1, . . . ,m. 



2 _ h{z) - h{-z) 



q{z') 



q{h) - ± 



2z 
fe(±7A") 



I = 1, 



h{-Zk) 



1(4) 



h(-Zk) 
2zk 



fc = 1, . . . , 2n, 



(1.32) 
(1.33) 

(1.34) 



where Zk are the zeros of the polynomial h. Since deg 5 = m by assumption, we conclude that bn-m 7^ 
but bo = ■ ■ ■ = bn-rn-i = 0. Thus, (|1.26l) implies 



R{p,q) = bll_,nYlpiHj) 



(1.35) 
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Using ((Oil), PTMI) and PTM)) . we have 

m m 2n 2n 

(R(p, q))2 = C-™ n va^) - n «o n (a^^- - = n 

2n 2n , / X 2n 2ti 

= ai- n '^(-^) = n ^ = «o" n ^ ik-^ + 

fc=l fc=l A;=l 4=1 



2n 2n 



_ 2m+2n 



a, 







/c=l 1=1 l<i<fe<2n 



Thus, we obtain 

l<i<fc<2n 

To determine the sign, we consider the special case h{z) = (z — 1)^". Then 

, {z-lf^ + iz + lf- {z-lf-_{z + l)^ 

so that degp — n, degg = n — 1. The roots of the odd part q can be determined from the equation 
[z — 1)^" ~ [z + 1)^" except that the zero root should be discarded. This shows that the roots of q are 
{w^ : A; = 1, . . . , TT, — 1} where wj. is defined from the equation 2/{'Wk + 1) = 1 — e"^, A; = 1, . . . , — 1. 
Consequently, 



n— 1 n—1 / ^ 

n^(-^) = n 

k=l k=l \ 1 — e n 



2n 



n — 1 " ^ ri-1 




, , sin"'( — ) sm^" (^) ' 

fe=l V n ; fc=l \ n J 

Thus, according to the last formula in p.26p . the sign of the resultant R{p, q) in our special case is 

n(n— 1) n(n+l) n(n— 1) 

(_l)«("-i)sign6'/(-l)^— = (-1)^— = (-1)"+^—. 
As we already established, the resultant of p and q is a, polynomial in the coefficients of p and q, hence in 

the coefficients of h. The expression (—1) 2 Oq" ni<fc(-2i + ^fe) is a symmetric function of the roots of 
h multiplied by its leading coefficient to the power 2n — 1, and hence also a polynomial in the coefficients 

n(n — 1) 

of h. Since the two polynomials must be identically equal, we conclude that the sign (—1) 2 occurs 
at all times whenever degp = n, deg q = n — 1. 

We now show how to produce a formula for the case m < n — 1 from the already established formula 
for m = n — 1. Given a polynomial q of degree n—1, set bi through bn-m-i to zero, thus obtaining a 
polynomial of degree m. Observe what happens to the determinantal expression (jl.l6p . Since the lower 
left block of size n x (n — m — 1) is now filled with zeros, the upper-triangular submatrix above produces 
the factorization 

R(p,g) = a^'"-lR(p,g), 

where R(p, 5) denotes the "old" resultant of p and q constructed as if degg were equal to {n — 1). Thus 
the "new" resultant R(p, q) satisfies the equation 



l<i<k<2n 

hence R{p,q) = (-l)'^ao+" ni<.<fe<2„(^^ + ^j)- 
Our next statement can be proved analogously. 
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Theorem 1.15. Let the polynomials p and q be given by (jl.2p - (jl.3p . and let degq = m < n = degp. 
Then the resultant of these polynomials can be computed by the formula 



■t(n-l) 



(1.36) 



l<i<fc<2n+l 

where Zi (i = 1, . . . , 2n + 1) are the zeros of the polynomial 

g{z)^qiz^)+zp{z^). 



The famous Orlando formula is a simple consequence of Theorems ll.l4f[LT5l Before proving the 
Orlando formula, we introduce the following determinants for the polynomial (|1.2p 



(1.37) 





oi as 


as 


aj . . 


• a2j-i 








ao 02 


a4 


a& ■ ■ 


• a2j_2 






A,(p) = 


ai 


as 


as . . 


• a2j_3 






ao 


a2 


04 


a2j-4 


, J = 


1, . . . ,n, 










.. 


Oj 






where we set : = for i > n. 














Definition 1.16. The determinants A,(p) {j 


= 1,- 


■ ■ ,n) are called 


the Hurwitz 


Hurwitz minors of the polynomial 


P- 












It is easy to see that the polyn 


omial p 


can be always represented as 


follows 




p{z) = 


Po{z') 


+ zpi(z^) 






where 
















Po{u) = 




+ a^u'^ 


-1 + ...^ 








Pi{u) = 




+ a2u' 


-1 + ...^ 


- a,i-i, 





(1.38) 



(1.39) 

+ a„_i, 

if the degree n of the polynomial p{z) is odd: n — 21 + 1, and 

Poiu) — flQu' + a2U^^^ + • ■ ■ + a„, 
Pi(u) = aiu'""^ + a^u^^^ + • • • + a„_i, 

if n = 2L 

Theorem 1.17 (Orlando, |73[ 136] ). Let the polynomial p of degree n be given by (jl.2p . Then the deter- 
minant A„„i(p) defined by (|1.37p can be computed from the formula 



(1.40) 



n(n— 1) 

A„_i(p) = (-l)^^a^l n (^»+^.^-), 

1<2<J <n 



(1.41) 



where Zi, i = 1, . . . ,n, are the zeros of the polynomial p. 
This equality is known as the Orlando formula. 



Proof. At first, let the degree n of p be odd, n = 2/ + 1. Then ([L38l) - (|09l) show that degpo < I, 
degpi — I, and the leading coefficient of pi is ao. Thus, p.l2p . ()1.26p and (|1.36p imply 



1(1+1) , , „ 
A„_i(p) = (-l)'V2,(pi,Po) = (-l)^a[,-'^°s^"R(pi,po) = (-l)'ag' H + ^^■)' 

l<i<j<2l+l 



which coincides with (|1.4ip since n = 2Z + 1. 
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If n = 21, then (jl.38p and (|1.40p imply degpi < I — 1, degpo = and the leading coefficient of 
the polynomial po is oq. As above, the formulas (I1.12[) . (|1.26l) and (ll.30|) can now be combined to obtain 

i<i<j<2; 

which is exactly the formula (|1.4ip . □ 



1.3 Euclidean algorithm, the greatest common divisor, and continued frac- 
tions: general case 

In the previous subsection, we considered the resultant of two polynomials and observed that it is equal 
to zero if and only if these polynomials have a nontrivial common divisor. The standard way to find their 
greatest common divisor is via the Euclidean algorithm. 

Let us consider polynomials p and q given by (|1.2p - ()1.3p and let us denot^ 

foiz):=p{z), fi{z):=q{z) - —p{z). 



Construct a sequence of polynomials /o, /i, ■ • ■ , //c (k < n) by the following formula 

/j-i(^) = + i = {fk+i{z) = o), 



(1.42) 



where qj is the quotient and fj+i is the remainder from the division of fj~i by /_,-. The last polynomial 
in this sequence, fk, is the greatest common divisor of the polynomials /o and /i (and also all other 
polynomials fj in the sequence). In other words. 



where hk{z) = 1. Now, denote 
Rewriting (|1.42l) . we obtain 



fj{z) ^ hj{z)fkiz), j = 0, 1,. 

_ hj{z) 



R,iz): 



Rj{z) 



1 



(1.43) 
(1.44) 



j = l,...,fc. 



gj(z)+i?j+i(z) 

where Rk+i{z) = 0. Using this equality, we can represent the function Ri{z) as a continued fractioi|f|: 

h{z) hAz) 1 



(1.45) 



Ri{z) 



h{z) ho{z) 



qiiz) + 



1 



(1.46) 



q2{z) 



1 



qsiz) 



qk{z) 



It is easy to see that, for each j = 0, . . . , fc— 1, the polynomial hj is the leading principal minor of order 
k~j of the following k x k tridiagonal matrix 



(qk{z) -1 
1 qk-l{z) -1 
1 qk-2{z) 





V 



\ 






q2{z) -I 
1 qi{^)J 



(1.47) 



^Thus, deg /i < dog/o. Obviously, if degg < degp, that is, if fep = Oi then f\ = q. 
■^The functions R = q/p and _Ri are related via the formula R{z) = R\{z) + bo/ao- 
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In particular, hoiz) = detj'{z). Thus, the determinant of the matrix J^{z) is the denominator of R{z). 
This is a very useful observation, as certain properties of the function R turn out to be connected to the 
location of the eigenvalues of the generalized (matrix polynomial) eigenvalue problem 

J{z)u^O. (1.48) 

Conversely, the behavior of the eigenvalues of the problem (|1.48p can provide information about certain 
properties of the function R. Later in this paper we will give examples of such interrelation. 
More generally, in addition to the fraction (11.46^ we may consider the fractions 

hj-i{z) 1 

+ 1 

qj+iiz) + : 



<lk{z) 



The continued fraction expansion (jl.46p can be found efficiently in terms of the corresponding Hankel 
minors, as the next theorem shows. 

Theorem 1.18. Two rational functions R and G vanishing at oo are related via the identity 

g{z) + G{zy 

where g is a polynomial of degree m(> 1) with leading coefficient a 

g{z) ^ az"" + ■ ■ ■ (1.50) 

if and only if 

Di{R) ^ D2iR) = ■ ■ ■ = Dm-i{R) ^ 0, {when m > 1) (1.51) 

and 

m(m-i) (—^^3^) (n) 
D„,+,(i?) = (-l)^^.L_L_|^, J- = 0,1,2,..., (1.52) 



where Dq(G) : = 1, and the determinants Dj{R) and Dj{G) are defined by (jl.6[) . 

Proof. Since the functions R and G vanish at cx), they can be expanded into Laurent series 

R{z)^— + — + ■■■ + — ^ + + , G(z) = — + — + — + ••• . (1.53) 
Moreover, the conditions (|1.49p ~ (ll.50p hold if and only if 

So = si = • • • = s„i_2 = and Sm-i 0. (1-54) 
In fact, if R satisfies prig)) - (IT3g|) . then 

5, = lim z'+^R(z) = lim ^ ^ = 0, i = 0, 1, 2, . . . , m - 2. (1.55) 

z-s-oo 2->oo g[z) G(Z) 

^ + 

and 

s„,_i = lim z^Riz) = lim ^ ^ ttt^ = - ^ 0. (1.56) 

z^oo z-i-oo 7lZ +---+7m G[z) a 
a H 1 

Now assume that the condition (11.541) holds and that R(z) — ——- -— for some polynomial /. If 

f{z) + G{z) 

deg / = j for some 1 < j < m, that ifQ, if f{z) = jz^ + • • • , where 7 7^ 0, then Sj-i — lim z^ R{z) = — 7^ 0, 

*The degree j of / cannot be 0, since the limit lim R{z) = — is nonzero for a constant nonzero function /(z) = 7. 
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contrary to (|1.54p . On the other hand, if the degree of / is greater than to, then s„i_i = hm z"^R(z) = 0. 

2^00 

Thus, / must be of exact degree to since otherwise a contradiction arises with one of the conditions (|1.54p . 
Also note that the equahties (|1.54p are equivalent to (|1.5ip . Moreover, from (|1.54p we have 

Sm-l 








S2m-4 S2m~3 
S2m-3 S2m-2 



m(m — 1) 



-^m— 1 

For convenience, denote the coefficients of the polynomial g as follows 

g{z) =: az"" + + t_,„+i2'"-2 + . . . + t^2Z + t-i. 

If the functions R and G satisfy (|1.49p - (|1.50p . then, according to (ll.53p - (|1.54p . we have 



Sm + 1 



This identity implies the following relations: 

1 

— 1 1 

a 



if 



J = 0,1, 2,. 



= 1. 



(1.57) 



(1.58) 



Now, the equality (fT32l) for L»™(i?) follows from (fL57l) - (fL58l) . 
For a fixed number j > 1, consider the determinant Dm+j{R)'- 



Dm+j{R) — 



s„_i 



-^m — 1 -^m 

-^m— 1 -^m -^m+l 

S2m-4 S2m-3 S2m-2 

S2m-3 S2m-2 S2m-1 

S2m-2 S2m-1 S2m 

•S2m+j-3 S2m+j-2 S2m+j-l 



-^m+j — 1 

S2m+J-3 
S2m+j-2 
S2)n+j-l 

S2m+2j-2 
t 



Add to each column [i = m-\-j, m+j—1, ...,3,2) columns i— 1, 2, . . . , 1 multiplied by 

tj-3 ii-2 



(1.59) 



, respectively. This eliminates the entries in the upper right corner of the determinant (jl.59p 

using (jl.58p . So, we can rewrite the original determinant as a product of the following two determinants 
of order to and j, respectively: 



where 
















^771— 1 




dii 


di2 ■ 


■ dij-i 


du] 



















d2i 


d22 ■ 


■ d2j-i 


d2.j 




_,(i?) = 





Sm-l 










d3i 


d32 ■ 


■ d3j^i 


d3.,3 






^m — 1 















dj2 ■ 




d]] 




■ — ~Sm- 




^12-1 


-^m+zi— 3 ' 

a 






tii+i2-3 


t 




-11-2 • 


a 




Sm 


a 


^7n — l 


a 



(1.60) 



From ([T38)) and (fOOl we obtain 



An+,(i?) = (-1)- 



Sm— 1 
Sm 





^7n— 1 



Sm.-\-j — 2 ■^m+J — 3 



5m-l 



to 

h 



h 

t2 



tj-2 t]~l 
tj-i tj 



tj-2 
tj-l 



t2j-i t2j-3 
t2j-3 t2j~2 
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Conversely, if the equalities (|1.51l) - (|1.52p hold, then we obtain sg = si = • • • = s„i_2 = from (|1.51l) by 
induction, which was already proved to be equivalent to the fact that R satisfies (jl.49p , with a polynomial 

m(m— 1) 

g of degree at least m and Dm{R) — (— 1) 2 Sm_i. If degg > m, then 

Sm_i = lim z'^Riz) = lim ^ x = 0, 

z-s-oo z^oo g[z) Cj[Zj 

~,ni 

m(m— 1) m(m— 1) ]^ ]^ 

therefore, Dm(R) = (-1) 2 s™_, = (-1) 2 = and - = 0, according to Ol. Thus, 

a™ a 

from (|1.5ip - (|1.52p we get Di{R) — for all i G N, which is impossible, since i? is a rational function, hence 

at least one of the minors Di{R) must be nonzero due to Theorem 11.21 Consequently, the polynomial g 

satisfies ([T30)) with a^O. D 

If a rational function R with exactly r poles is expanded into a continued fraction (|1.46p . then the 
rational functions Rj defined in (|1.44p satisfy the relations (|1.45p . where 

= +••• , aj^O, j = l,2,...,fc. (1.61) 

Here ni + n2 + ■ ■ ■ + rik = r and rii > 1, since deg fi < deg /o, as is remarked above. The other degrees Ui 
are greater or equal to 1 due to the structure of the Euclidean algorithm p.42p . 

From (jl.5ip - (|1.52p we obtain, for a fixed integer j (j = 1, 2, . . . , fc), the following product formulas: 

»i(tii-1) Y1 "i 1 

ni+2 iz 

j j 

ni(ni-l) Ui 1 n2(n2-l) 

= (-l)^^(-l)^- — •(-l)^^(-l)-- X 

ni+2 Tii 

■ -On3+n4 + ---+nj (^s) = • • • 

"2+2 23 

"2 



This chain of equalities results in the formula 



D„,+„,+...+„^.(i?) = n(-l)^^ -(-1)*^" -IT 1 ' J-l,2,...,fc. (1.62) 



i=l i=l "i+2 Yl 

p=i+l 

Remark 1.19. Our discussion above can be summarized as follows: Suppose that a rational function 
R with r poles has a continued fraction expansion ()1.46p with polynomials qj satisfying ()1.6ip . Then 
Theorems II. 2[ [TT51 and [T.18I and formula; (11.621) imply that all Hankel minors Dj{R) are equal to zero, 

except for the minortQ D„j (i?), _D„j+„2 (i?), . . . , Dn-i+n2-\ i-nfc(^)j which are not zero and which can be 

calculated from the formulas (ll.62p . 

Using Theorem 11.181 and the formula (|1.62p , we now describe equivalence classes of rational functions 
whose sequences of Hankel minors are the same. 

Theorem 1.20. Two rational functions 

R[z) = h ^ + ^ H and G(z) 



^Recall that ni + n2 + ■ ■ ■ + rii^ = r. 
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both vanishing at infinity, have equal Hankel minors 

D,{R) = D,{G), J = 1,2,... , (1.63) 
if and only if their continued fraction expansions 

R{z) = and G{z) = ^ 



+ ^ qi{z) + ^ 

q2{z) + ^ q2{z) + Y 

qsiz) H — qsiz) H 



fti(^) qk2{z) 



satisfy 



ki^k2=:k, (1.64) 
and the polynomials qj and qj, for each j (j = 1, 2, . . . , fc), have equal degrees and equal leading coefficients: 

qj{z) = q;,-z"^' H , 

J „ j = l,2,...,fc. (1.65) 

Proof. If the continued fraction expansions of the functions R and G satisfy the conditions (|1.64p - (|1.65p . 
then the equaUties (|1.63l) fohow from Remark 1 1 . 1 9 1 and the formulas (|1.62l) . 

Conversely, let the Hankel minors associated with the functions R and G satisfy (|1.63l) . Therefore, by 
Theorems 11.21 and 11.31 the functions R and G have equal number of poles. Moreover, they can be then 
represented as follows: 

R{z) = Ri{z) = and G{z) ^ Gi{z) = ^ 



qi{z) + R2{z) ' ' ' qi{z) + G2{z)' 

It follows from Theorem 11.181 that the degrees and the leading coefficients of the polynomials qi and qi 
coincide. According to (|1.51l) . the functions i?2 and G2 have equal Hankel minors. So we can apply the 
same argument to them. Thus, Theorem 11.181 shows that the degrees and the leading coefficients of each 
pair of the polynomials qj, qj are equal to each other and, consequently, the number of those polynomials 
must be the same as well. Since the functions R and G have equal number of poles as was proved above, 
the equality (|1.64|) thus follows. □ 

Remark 1.21. Note that the equalities Dj(R) — Dj{G) do not imply the equality of the functions R = G. 
Counterexamples are quite easy to construct. For instar 

^, . 1 1 
R{z) = 



z 

Giz) = - 



z 

1 1 



;, we 


can 


take 


1 


2 


5 








1 


2 


6 


— ^ 


l3 





1 z 
z-1 

z — 3 

Then i? ^ G but Di{R) = D2{R) = Di{G) = D2{G) = 1. 

Finally, let the function R have a continued fraction expansion (|1.46p and also a Laurent series expansion 

= f!ii^ + ^ + f!ii±i + . . . , (1.66) 

where ni — deggi > 1. 
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Consider the functions 



(1.67) 



Qjiz) 



constructed using the polynomials (|1.6ip by the Euclidean algorithm (|1.42p . 

Definition 1.22. The polynomials Qj are called partial numerators, the polynomials Pj partial denomi- 
nators, and the fractions Fj partial quotients of R. 

The denominator Pj{z) of the fraction Fj{z) is the jth leading principal minor of the matrix 





fqiiz) -1 

1 q2{z) -1 
1 q^iz) 





V 



\ 






qk-i{z) -1 
1 qk{z)J 



and Pk{z) = ho{z) (see (|1.46p V Let nij denote the sum of the degrees ni through nj (see (|1.6ip ): 

mj — ni + n2 + ■ ■ ■ + nj, j = l,2,...,k. (1.68) 

Then deg Pj = mj . 

Notice that, for a fixed number j (1 < j < fc), the initial terms of the Laurent series of the function Fj 
coincide with those in the Laurent series ()1.66|) of R up to, and including, the term ^'"^ ^ : 



Fiiz) 



Qjjz) _ Sm~i 



„0) 



v2m 



(1.69) 



In fact, each coefficient Si of the series (|1.66|) can be found from the recurrence relation 

s,^ lim fz*+ii?(0) - s„i_iz*""i+i - s„,z'+"i s,_izl , i = ni - 1, ni, ni + 1, . . . (1.70) 

Using the expansions (11.461) and (|1.67l) of the functions R and Fj, respectively, together with the for- 
mula prrni) . we obtain pTMl) . 

To find an explicit formula for the polynomials Pj that depends only on the coefficients Si, we introduce 
the following notation: 



Pj{z) Poj^™' + Pl^jZ-^'-^ + • • • + Pr,^,-l,JZ + Pra^^j, 

For a fixed number j between 1 and k, the formula (jl.69l) implies 



j = l,2,...,k. 



Qj{z) = P,{z) 



We must require that the coefficients of z^^ , z^^, . . . , be zero, which leads to the system 



(1.71) 



SlPrrijJ + S2Pmj-l,j 



Snij-lPlJ + SrnjPo,j 
SrrijPlJ + Smj+lPo,j 



(1.72) 



^2m4 — 
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where we set so := si := • • • := s„j_2 0. By Cramer's rule, the solution to the system (ll.72p satisfies 

Si 



p. 



50 Si 

51 Si 



Si-1 Si+1 
Si Si+2 



Smj 
Smj+l 



Smj— 1 -^rrij • • • ^mj-^-i—l '^mj+i • • • ^2mj~l 

This formula implies the following representation for the polynomials Pj : 



i = 0, 1, 2, . . . , rrij — 1, 
j = l,2,...,fc. 



50 Si S2 

51 S2 S3 

Smj— 1 SyUj SyUj+I 
1 Z z2 



SrUj 
Smj +1 

S2mj-1 



Here, according to the notation (|1.6ip . 



^'o,j=]Jai, j = l,2,...,fc. 



j = l,2,...,k. 



(1.73) 



This formula, combined with ()1.62p and the notation ()1.68p . yields 

1 



Po. 



2j 



0,2j + l 



"2i-l("2i-l-l) 

(-1) 2 Drn,^_,{R) 



Dm2i-2 {R) 



for J = 1,2,..., 



i=l 



"2i-l("2i-l-l) \ n2i-l 

2 Dm^^-AR) 



»2t("2i-l) 

( -1) 2 



1 

"2i("2i-l) \ n2i 

(-1) 2 



1 

"21 



Dm2^-l{R) 



Dmii-AR) 



"2j + l("2j + l-l) \ n2i+l 

( -1) 2 An2,(i? ) 

-Dm2j + i(^) 



for j = 0,l,2. 



In case the Euclidean algorithm is regular, the notion that will be introduced in the sequel, these formulae 
take a simpler form. 

The coefficients of Qj can be obtained from the coefficients of Pj via the formula (ll.7ip : denoting 



we find the coefficients of Qj by the following 



Qi.j — 



D,n, [R) 





50 Si 

51 S2 



So Si 

^rnj —ni — i+l ^rrij — ni —2+2 
^rrij —ni —2+2 ^rrij —ni —i-\-3 



S2 



m j — n 1 — z 



S2mj —n-i —2+1 



Sni+i — 1 
Sjnj 
Snij+l 

S2m,-1 



for z = 0, 1, . . . , rrij — ni. 



Note that the rational functions Fj{z) = ^'''^^^ are exactly the diagonal elements of the Pade table of 



the function R, see pUl 
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1.4 Euclidean algorithm: regular case. Finite continued fractions of Jacobi 
type 

In this section we discuss the best known case of the Euclidean algorithm, which leads to orthogonal 
polynomials, 3-term recurrence relations and other related phenomena. We give conditions for regularity 
and discuss the form of partial quotients and the generalized eigenvalue problem corresponding to the 
regular case. 

Suppose that all the polynomials qj resulting from an application the Euclidean algorithm (|1.42|) are 
linear: 

^(2) = "j^ + "i) /3j e C, aj y^O, j = 1, . . . , r. (1.74) 
We call this situation the regular case of the Euclidean algorithm. In the regular case, 

degfj{z) = n-j, j = 0,...,r, (1.75) 

and r = n — I, where I is the degree of the greatest common divisor fr of the polynomials /o and fi (and 
of all other polynomials fj in the sequence). 

Thus, the polynomials fj satisfy the following three-terms recurrence relation: 

f,-^{z)^(a,z + P,)f,{z) + f,+,{z), j = l,...,r. (1.76) 

Consequently, the function R expands into the continued fraction 

„/ s f}{z) h-\(z) 1 

R{z) = > 1-77 

/o(z) hfi{z) 1 

aiz + Pi^ 

012Z + /32 H ; 



a^z + /33 + 



UrZ 



where the polynomials /lo and hi are defined by (jl.43p . 

Definition 1.23. Continued fractions of the type (|1.77p are usually called J-fractions or continued frac- 
tions of Jacobi type. 

Remark 1.24. If a rational function G satisfies the condition lim G{z) ~ c, < \c\ < 00, then we will 
say that G{z) has a J-fraction expansion if the function R{z) = G{z) — lim G{z) can be represented as 
in (fTTTl) . 

Thcorcm ll. 181 with some simple modifications implies the following corollary, which will be useful later. 
Corollary 1.25. Two rational functions R and G, both vanishing at 00, satisfy the following condition 

az + p + G(z) 

if and only if 

Dj{R)^^^]^D,.i{G), J = 1,2,..., (1.78) 

where Dq{G) -.— 1. 

Using Corollarv ll.251 we can now prove a criterion when a rational function expands into a J-fraction. 
Theorem 1.26 ([74]). A rational function 

i?(z) = s_i-H^ + ^ + ^-K--- (1.79) 

z z'' z-^ 

with exactly r poles has a J-fraction expansion if and only if 

Dj{R)^0, j = l,...,r, (1.80) 
where Dj{R) are the Hankel determinants defined in (jl.6p . 
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Proof. In view of Remark 11.241 it is sufBcient to consider the case s_i = 0. Suppose that the func- 
tion R has a J- fraction expansion (jl.77p . i.e., a continued fraction expansion (I1.46P with polynomials qj 
satisfying (|1.74p with k — r. Then the formula (|1.62p yield 

D,{R) = {-1)^\{ ^0, j = l,2,...,r. (1.81) 

^=l 

Now suppose that the inequalities (|1.80l) hold. Then Di{Ri) = s'^q'^ ^ 0, where we denot^ 

Ji) Ji) Ji) 
i?i(z) -.^ R{z) = ft^ + fi^ + f2^ + • • • 

Therefore, R\{z) can be represented as follows 

aiZ + /3l + i?2(2) s^q' 
(-l)J-l 

Now, Corollary 11.251 implies that Dj(R\) ~ — 2 -\ ^j-i(^2)i j — ■■ Therefore, the function R2 is 

of the same type as Ri, that is, Ri satisfies the conditiont|3 

DjiR2)^0, j = l,2,...,r-l, 
which are analogous to (|1.80p . In particular, Sg 7^ 0, where 



i?2(^) 



J2) J2) (2) 



If we continue this process, then for each function 



„(i) „(i) Jj) 

we obtain that Sq'' ^ 0, j — 1,2, . . . ,r. Consequently, 

Ri(z) = -— , 7 = l,2,...,r, 

where Rr+i = 0. This means that the function R can be represented as a J-fraction (see (ll.45p - (|1.46p ). □ 

This theorem and Theorem 11.51 imply the following statement, which was proved in |101) (Theorem 
41.1) with inessential differences. 

Corollary 1.27. Let R{z) — ^\ \ , where p and q are defined in (jl.2l) - (|1.3p . and let the function R have 
exactly r poles. The function R can be expanded into a J-fraction (|1.77p if and only if 

V2,(p,g)^0, J = l,...,r, (1.82) 

where \J2j{p,q) are defined in (|1.12p . 

f (^) 

Corollary 1.28. Given a rational function R{z) ~ , where deg /i < deg/o, the Euclidean algo- 

rithm (|1.42p applied to the polynomials fo and fi is regular if and only if the inequalities (|1.80p hold. 



^Recall that we assumed s_i = 0. 

^For j > r, we have = Dj{R2) = 0, and therefore R2 has exactly r — 1 poles. 
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If we make an equivalence transformation of the fraction (ll.77p . as in [1011 p. 166], and set 

do:=— , d,:^ , j = l,2,...,r-l. (1.83) 

■■=- — . j = l,2,...,r, (1.84) 
then the J-fraction (|1.77p takes the form 

R{z) = '^-^ . (1.85) 

di 

z — ei 



2—62 



dr-1 



Moreover, from ([LM|) and we have (see (TUTJ p. 167]) 

i 

D,{R)=Dj_^{R)-\[d,.x. (1. 

In fact, 

jQ-i) ^ 1 j(j-i) 1 i=3 1 1 

1=1 1=1 » i=l 

= (-i)^(-i) -2 i-iy-'- ■ n ^.-i(^) = ^.-i(^) • Il'^-i- 



l—l l—l 



The expression (jl.86D implies an interesting formula: 

Dm 



- '-'':i7^r ' ^ ^■-0,l,...,r-l, (1.87) 



where we set D^i{R) -.^ Do{R) -.^l. 

Remark 1.29. From the formulae (II. 83^ . (I1.84p and (I1.87P it follows that both forms of J-fraction expan- 
sions (|1.77p and (|1.85l) of the function R{z) are unique, that is, all coefficients in (|1.77p and (|1.85l) are 
defined uniquely for the function R. 

Remark 1.30. (cf. [lOlj p. 170]) Suppose that the function R has a J-fraction expansion (|1.77p . Then 

R{-z) = i ^ . (1.88) 

aiz - /3i -I 

a2Z - /32 H 



1 



UrZ — Pr 



According to Remark ll.29[ (a properly normalized) J-fraction expansion is unique. Therefore, if the 
function R is odd, i.e., R{z) = —R{—z), then aU /3j in (|1.88p must be equal to zero. Conversely, if (5j are 
all equal to zero, then R{z) is obviously an odd function of z. 



Thus, the formula p. 861) implies 



i=Q 

In other words, if the function R has a J-fraction expansion (|1.85p . the minors D.j{R) do not depend on the 
coefficients e^, j = 1, . . . , r, whereas the minors Dj{R) obviously do. In order to establish the dependence 
of these minors on the coefficients of the fraction (|1.85p . we first prove the following simple fact. 
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Lemma 1.31. Let the complex rational function 

^^^^ z z'- z-^ 
with exactly r poles, has a J -fraction expansion (|1.85l) . Then 

Dr{R) 



Dr{R) 



n 

fe=l 



(1.90) 



(1.91) 



where iik are the poles of the function R, and the minors Dj{R) and Dj{R), j — 1,2,..., are defined 
by (mH) and (IL^ . 

Proof. If R has a pole at zero, then the lemma holds true, since m this case Dr{R) 7^ by Theorem II. 26[ 
while Dr{R) = by Corollarv ll.4l Now assume that R has no a pole at zero. 

At first, suppose that the function R has only simple poles /ifc 7^ /ij whenever k ^ j, so it can be 
represented as a sum of partial fractions 



k=l 



z- fJ-k 



This formula together with (|1.90[) gives the following well-known formulae 



Si = 



E 

k=l 



i^kfj,i, 



J = 0,1,2,. 



(1.92) 



On the other hand, (|1.22l) implies 



DriG) ^Yl'^k 
fe=i 

Further, recall that Dr{R) = Dr{^), where 



1 1 

Ml 



(1.93) 



$(z) = zi?(z) = E-. + E^^ = 



k=l 



k=l 



z- f^k 



3Q 



r — 1 
Mr 



Si S2 . S3 



So, analogously to ()1.93p . we obtain 

r 

i5,.(i?) = A.(*) = n^^fe^'^) 



fe=l 



1 


1 


1 


Ml 


/i2 


Mr 






r—1 
• Mr 



(1.94) 



Now (fTMl) follows from (fTMl and (fL94l) . 

Let R — — and suppose that the function i? has multiple poles, or, equivalently, that the polynomial p 
P 

has multiple zeros. In this case, we can consider an approximating polynomial p^ with simple zeros such 
that 



Then 



where Re{z) 



lim ^^(z) = p(z) for all z. 
e— >o 



Dr{Re) > Dr{R), Dr{Re) > D,. {R) , 



g(g) 
Peiz)' 



The formula (|1.9ip is valid for the polynomials q and Pe, so it is also valid at the 



limit, i.e., for the polynomials q and p, since the product of all zeros of the polynomial Pe{z) tends to the 
product of all zeros of the polynomial p whenever e — > 0. □ 
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Let us consider the following tridiagonal matrix 









V 



62 








63 






6r-l 



\ 






■\/ 



(1.95) 



constructed using the coefBcients of the J-fraction (|1.85p . 

Suppose that the coefficient rfo in (|l-85l) equals 1 and consider the partial quotients of the J- 
fraction 



z — e\ 



di 



j = l,...,r. 



(1.96) 



Z - 62 



Then — R and the polynomial Pj is the characteristic polynomial of the leading principal submatrix of 
Jr of order j, for any j = 1, . . . , r. 

Using the formula (|1.9ip . it is now easy to prove the following theorem. 

Theorem 1.32. Let the matrix J„ be defined by (jl.95p and let the function R defined by (|1.90p have a 
J-fraction expansion ()1.85|) (with do — 1). Then the leading principal minors \Jr\\ , m — 1, . . . ,r, of the 
matrix Jr can be found by the following formults: 



\Jr\T 



Dm{R) 

D„,{Ry 



ni = 1 , . . 



(1.97) 



Proof. At first, we note that the formula (|1.97p holds for det(Jr) — \ Jr\i- Indeed, on the one hand, the 
determinant of the matrix Jr is the product of its eigenvalues. On the other hand, the eigenvalues of the 

Dr{R) 



matrix are the poles of the function R, so by (|1.9ip . their product equals 



DriR)' 



Let us now turn to the functions Fj introduced in (jl.96p . According to (|1.69p . we have 



Pv^{z) 



So 

z 



£i 

^2 



S2m-1 



Z 



2m 



(m) 
2^2m+l 



m = 1, 



(1.98) 



where the coefficients Si, i — 0,1, ... ,2m — 1 coincide with those of the function R defined in (|1.90p . 
Consequently, 

Dj{Frn) = Dj{R), Dj{Frn) = Dj{R), J = 1, . . . , m. 

Thus, applying Lemma ll.31l to the function Fm and to the leading principal submatrix of J7r of order m, 
we obtain 

Dj{Fm) _ D,n{R) 



\Jr\T 



DjiF^) D„,{R)' 



m 



1, 



as required. 



□ 



Now suppose that the function R has a J-fraction expansion (jl.85p . where do may differ from 1. Then 



we can consider the function G{z) 



R{z) 



. It is clear that 



An(-R) = dS'A„(G), D,{R) = d^W.iG) 
Consequently, according to (|1.97p . we get 

Dm{G) _ DrniR) 



m ~ 1, 



\Jr\T 



D^{G) D^iR)- 



m = 1, . . . ,r. 



(1.99) 
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Thus, if the function R has a J- fraction expansion (|1.85p . then it foUows from (I1.99P and (jl.97p that 
the minors Dj{R) can be found as foUows 



Dm{R) = \Jr\T ■ n < 



1=0 



where |i7r |'" is the leading principal minor of order m of the matrix J7r, m = 1, . . . , r. 

Finally, note that the matrix J^j can be replaced throughtout the discussion above by the matrix 



/ei 


1 


. 


0\ 


di 


62 


1 . 








d2 


63 • 











. 


• 6^-1 1 


Vo 





. 





which is diagonally similar to the matrix via the transformation 



diag(l, \fd[, \Jdid2, 
which also preserves all principal minors. 



Let us now return to the equation (jl.48p . If all the polynomials qi in the fraction (I1.46P are linear, as 
in (|1.74p . then the equation (|1.48p becomes a generalized eigenvalue problem 

(Az + B)m = 0, (1.100) 

for the matrix pair (A, B) of order r where A is diagonal and B is tridiagonal (for a similar setup, see [38]): 



( ar 





VO 





Ctr-l 











. 






(fir 


-1 





. 


0\ 





. 







1 




-1 . 


. 





ar-2 ■ 


. 










1 




. 











, B = 













■ 61^2 
















. P2 


-1 





. 


ai) 




^0 








1 


/3i/ 



(1.101) 



and the function R expands into a J-fraction p.77p . 

The polynomials hj (j = 0, . . . , r— 1) in (I1.77P are known to be the leading principal minors of order r—j 
of the matrix J{z) — zA + B as we mentioned above (see also [11] [34]). In particular, hi^{z) = det(zA + B). 

As we already mentioned in Section 11.31 we can localize the eigenvalues of the problem (jl.lOOp using 



properties of the function R. For example, if aj > 0,/3j > for all j = 1, 



in (|1.74p (or, equivalently. 



in ()1.77p ). then every polynomial qj is a function mapping the closed right half-plane into the open right 
half-plane. In fact, Regj(z) — aj Rez-|-/3j > 0, whenever Rez > 0. Now note: if functions Fi and F2 map 

the closed right-half plane to the open right-half plane, then so do the functions Fi ~\- F2 and — (j 1, 2). 



Consequently, in this case the function R represented by (jl.77p also maps the closed right half-plane to 
the open right half-plane, being a composition of such maps. But since the function R has a positive real 
part in the closed right half-plane and is finite there, then all its zeros and poles and, subsequently, all 
zeros of the polynomials /lo and hi lie in the open left half-plane. In summary, if aj, l3j > 0, j — 1,2 . . . ,r, 
then the eigenvalue problem (|1.100p is stable, that is, all its eigenvalues lie in the open left half-plane. This 
example constitutes the subject of Problem 7.1 in 11 . A similar result was obtained in [37\ . 

Finally, at the end of this subsection, let us discuss the form taken by the partial quotients Fj defined 
by (|1.67l) in the regular case of the Euclidean algorithm. Since all polynomials qj are of degree one 
(see p.74p '). we have k ~ r, the number of poles of R. Moreover, for a fixed integer j {I < j < r) we have 
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deg Pj = rrij = j, where Pj is denominator of the fraction Fj of the form (cf. (I1.73p ) 



P( 



50 Si S2 

51 S2 S3 



Sj + 1 



S2j~l 



1,2,, 



The leading coefficients Pqj can be determined by the formulae (see, e.g., [67] ) 



Pr 



2j + l 



i-iy 



i-iy 



Dl{R)Dl{R)---Dl_,{R) 
Dl{R)Dl{R)---Dl_^{R)D2j{Ry 

Dl{R)Dl{R)---D%{R) 
Dl{R)DUR) ■ ■ ■ Dl_,{R)D2j+i{R) 



j = 0,l,2,.. 



- 1. 



1.5 Euclidean algorithm: doubly regular case. Finite continued fractions of 
Stieltjes type 

Assume, as above, that the rational function R has a series expansion (|1.79p . where s_i — 0, and consider 
the function 

Fiz) zR{z') = f£ + fi + + . . . (1.102) 

z z-^ Z'^ 

The function F can be also represented as a series 



Z Z^ Z'^ 



where 



*2j 

*2i+l 



0, 



i = o,i,. 



(1.103) 
(1.104) 



Remark 1.33. Let r denote the number of poles of the function R (counting multiplicities). Note that the 
function F has 2r poles if and only if R has no pole at zero. Otherwise, F has only 2r— 1 poles. 

Lemma 1.34. The following relations hold between the minors Dj{R), Dj(R) and Dj{F) defined in (|1.6p 
and (fTOl) .- 

D2,{F) = D,{R)-D,{Ky 

J = 1,2,..., (1.105) 



D 



2j-l 



D,{R) ■ D,_i{R), 



where we set Do{R) 

Proof. First interchange the rows and columns of the determinant (F) 



D2j{F) 



to 


h . 


• *2j-2 


^2j-l 




So 





Sl . . 







h 


t2 ■ 


• *2j-l 


t2, 









.. 







t2 


h ■ 


t2j 


hj+i 




Sl 





S2 . . 


S-j 





2j-3 


t2]-2 ■ 




tij-z 







Sj_l 


.. 





S2j- 


2j~2 


t2]-l ■ 


■ t^j^s 


tij-2 









Sj .. 


• S2J-2 





2j-l 


t2, . 


■ t4j-2 


tij-2 









.. 





S2j- 
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so that (2i--l)st row and (2i— l)st column move to the zth position, for each z = 2, 3, . . . , j; this produces 



D2AF) = 



So 


Sl . 













Sl 


S2 ■ 













S]-i 


Sj . 


■ S2J-2 














. 





Sl 


S2 


Sj 





. 





S2 


S3 


Sj+i 





. 







Sj + l . 





= D,{R)-D,{R). 



The formula for the determinants D2j-i{F) can be proved in the same way. 



□ 



This lemma, Theorem 11.261 and Corollary 11.41 can now be combined to derive the following corollary, 
which will be important later. 



Corollary 1.35. The function F defined in (|1.102|) has a J-fraction expansion if and only if 

D,{R)^Q, j = l,2,...,r, (1.106) 

5j(i?)^0, j = l,2,...,r-l, (1.107) 

Dj{R) ^ Dj{R) ^0, j = r+l,r+2,... (1.108) 

where r is the number of the poles of the function R, which generates F . In addition, Dr{R) = if and 
only if the function R (hence also the function F) has a pole at 0. 

Since F is evidently an odd function. Remark 11.301 shows that its J-fraction expansion (if any) has to 
be of the following form: 

F{z)^ , c, ^0, J = 1,2,..., ft. (1.109) 



1 



CiZ 



C2Z + 



1 



C32 



1 

CfcZ 



Here fc = 2r if the function R has no pole at and fc = 2r — 1 if it does, as follows from Remark 11.331 and 
Theorem 11.261 Making the equivalence transformation (|1.83p - (|1.84p in (|1.109p with Ci replacing and 
with Ci — f3i — 0, we obtain by induction 



C2j 



Y[d2^ 

i=0 



1,2, 



(1.110) 



i=0 



Cl 



i-1 

i=0 

II d2^ 

i=0 



J =2, 3, 



(1.111) 



since Cj+i 



Cjdj 



from (frSS]) . The formula; (fOT]) and (|1.105p imply 

D2^-i{F) ■ D2,+i{F) D,+i{R) ■ D,^i{R) 



DliF) 



D,{R)-D,{R) 



J'2i-l 



D2r-2{F) ■ D2r{F) _ A-l(i?) ' A(i?) 



D,{R)-D,_^{R)^ 



i = l,2,...,r. 



i = l,2,...,r. 



(1.112) 
(1.113) 
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Combining (|1.110p - (jl.llll) and (I1.112p ~ (|1.113p . we get very interesting and ultimately very useful relations 
between the coefficients Cj and the minors Dj{R) and Dj{R) (see [SI]): 

C2j^-^ , j = l,2,...,r, (1.114) 

D,MR)-D,iRy 

Remark 1.36. According to CoroUarv 11.41 the function R has a pole at zero if and only if Dr-i{R) ^ 
and Dr{R) = 0. From (I1.114p and (jl.l07p we conclude that C2r = oo in this case. 

Upon making an equivalence transformation in (|1.109p . replacing by z, and removing the factor z 
(as in |1Q11 p. 170]), we obtain 

R{z) = i , c,^0, where T^h \m\ < oo, ^^^^^^ 



CiZ - 



C2r-lZ if R{0) — OO. 



C2 



C3Z + 



1 

T 



Definition 1.37. Continued fractions of type (|1.116p are called continued fractions of Stieltjes type or 
Stieltjes continued fraction. Accordingly, if (|1.116p holds for a function R with r poles, then we say that 
R has a Stieltjes continued fraction expansion. 

Remark 1.38. If R{oo) — cq, where < |co| < 00, then we say that R has a Stieltjes continued fraction 
expansion whenever the function G{z) :— R{z) — cq has one. 

Summarizing all the previous results, we obtain the following criterion for a rational function to have 
a Stieltjes continued fraction expansion. 

Theorem 1.39 ([93l |94l ESI IMl HOH [62]). Suppose that a rational function R is finite at infinity, has 
exactly r poles, and can be represented as a series (|1.79p . The function R has a Stieltjes continued fraction 
expansion (|1.116p if and only if it satisfies the conditions (jl.l06p - (|1.108p . In that case, the coefficients of 
the Stieltjes continued fraction can be found from the formults (|1.114p - (|1.115p . 



Assume that 



F{z) = zR{z') = ^ (1.117) 

and F{oo) = 0. The function F has a J-fraction expansion (I1.109P if and only if the Euclidean algorithm 
applied to the polynomials go and gi is regular. Since F is an odd function, (jl.lOQp implies 

9j-i{z) ^ Cjzgj{z) + gj+i{z), Cj^O, j = l,2,...,fc, 

where k is equal to 2r— 1 or 2r depending on whether or not R has a pole at zero. Performing the 
transformation 

g2i-iiz) 



92i{z) ^ g2t{z), g2i-i{z)^ , j = l,2. 



we obtain 
where 



9j-i{z) = qj{z)g;j{z) + gj+i{z), qj{z)^0, j = l,2,...,fc, 

(1.118) 



q2iiz) = C2j, 1 = 1,2,... 

q2i-i{z) = C2j-iz^, i = l,2,...,r. 

If R and, subsequently, F have a pole at zero, then q2r{z) does not exist, according to Remark 11.361 
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F(z) 

„„ „„„ ^„„^ . ,, ,. _ is even. Therefore, both 

9o[z) 9a{z)_ z 

polynomials 50 and 171 are even, hence so are all subsequent polynomials gj{z), j — 2, 



Since F{z) 



^^^^^ is an odd function (see (|1.117p ). the function 



, k. Equivalently, 



9jiz) = fji^^)- 



From (jl.llSp we see that the polynomials qj{z) are even, so are functions of z^. Denoting qjiz'^) — qj(z), 
we obtain 



f,^,iz^)^q,iz')f,iz^) + f,+,iz^), J 



1,2,. 



Since 



fijz') Uz) Fjz) 

/o(2^) go{z) z 



= R{z^), replacing by z, we get 



R{z)^ 



h{z) 
Mz) 



and 
and 



fj-iiz) = qj{z)fj{z) + /j+i(z), qj{z) ^0, j = 1,2, . 



q2i{z) = C2i, j = l,2,... 

q2i-i{z) = C2i-iz, i = l,2,...,r. 
Recalling that we consider the case i?(oo) = 0, we see that the polynomials fj have fixed degrees 



(1.119) 
(1.120) 



h^{z) = h2^z''-' 

where n = deg fo>r and /ij 7^ for j = 0, 1, . 



1 = 0,1,2,..., 
z = l,2,...,r. 



(1.121) 



So, the function R{z) = 



h{z) 
Mz) 



has a Stieltjes continued fraction expansion (|1.116p if and only if the 



application of the Euclidean algorithm to the polynomials /o and /i has the form (|1.119l) - (|1.120p produces 
their the greatest common divisor fk{z), where k = 2r, if |-R(0)| < (X), and A: = 2r — 1 otherwise. We 
already know from the condition (|1.106p that the Euclidean algorithm (|1.119p - (|1.120p must be regular for 
the function R to have the expansion (jl.llGp . But in this case we obtain one more set of inequalities, 
namely, (jl.l07p . This justified calling such an instance of the algorithm doubly regular. 
Now consider again the rational function 



R[z) = — - = s-i H h ^ + ^ 

p(z) Z Z'^ Z"^ 



(1.122) 



where the polynomials -p and q are defined in (|1.2I) - (|1.3I) . 

We now introduce another (infinite) matrix associated with the function (|1.122l) . This object differs 
significantly from the Hankel matrix constructed in (jl.Sp from the coefficients Sj. In particular, this new 
matrix is made of the coefficients of the polynomials p and q. 

Definition 1.40. Given polynomials p and q from (|1.2p - (|1.3p . define the infinite matrix i/(p, q) as follows: 
if degq < degp, that is, if &o = 0, therjl 



iJ(p,g) 





ai 


0^2 


0-3 


04 


05 


...\ 





&i 


b2 


b3 


bi 









ao 


ai 


0.2 


as 












bi 


b2 


63 


bi 
















■7 



(1.123) 



° Generally speaking, h\ may be allowed to be zero. However, in this section we consider functions expanding into Stieltjes 
continued fractions, and Di(iJ) = so 7^ is one of the necessary conditions for such an expansion to exist by Corollarv ll.351 
At the same time, so 7^ implies b\ = aoso 7^ 0. 
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if degq = degp, that is, 60 7^ 0, then 





bi 


62 


63 


bi 


65 


...\ 





ao 


ai 


0-2 


03 


04 







bo 


bi 




63 


64 










ao 


ai 


02 


03 




V 












■7 



(1.124) 



The matrix (j) is called an infinite matrix of Hurwitz type. 

Remark 1.41. The matrix H{p,q) is of infinite rank since its submatrix obtained by deleting the even or 
odd rows of the original matrix is a triangular infinite matrix with oq ^ on the main diagonal. 

Together with the infinite matrix H(p,q), we consider its specific finite submatrices: 

Definition 1.42. Let the polynomials p and q be given by p.2p - (|1.3p . If deg q < degp = n, let H2n{p, q) 
denote the following 2n x 2n-matrix: 





/61 


62 


63 .. 


bn 





.. 









oo 


Ol 


02 . . 


o„_i 


o„ 


.. 













61 


62 .. 


• bn-1 


bn 


.. 








'H2n{p,q) = 





Oo 


fli . . 


0„_2 


On-l 


o„ . . 
















.. 


ao 


Ol 


02 . . 


a„ 













.. 





61 


62 .. 


bn 












.. 





Oo 


fli . . 


o„_i 


an/ 


= n, let 'H2n+i 


(p, q) denote the following 


(2n+i; 


X (2nH 


-l)-matrix 






/ Oo 


Ol 


02 . 


o„_i 


o„ 


. 










bo 


bi 


62 . 


■ • 6n-l 


6n 


. 













Oo 


ai . 


a„_2 


Ori-1 


o„ . 








'H2n+i{p,q) = 





bo 


bi . 


■ • bn-2 


6„-l 


6n • 
















. 


ao 


Ol 


02 . 


an 













. 


■ ■ bo 


61 


62 . 


bn 







\o 





. 





Oo 


ai . 


o„_i 


On 



(1.125) 



(1.126) 



Both kinds of matrices 'H2n{p,q) and 'H2n.+i{p,q) are called finite matrices oj Hurwitz type. The leading 
principal minors of these matrices will be denoted bjlf] Aj [p, q) . 

The infinite Hurwitz matrix H{p,g) has an interesting factorization property: 
Theorem 1.43. If g{z) = goz'- + giz'-^^ + ■ ■ ■ + 9i, then 

H{p-g,q-g)^H[p,q)T{g), (1.127) 
where T{g) is the infinite upper triangular Toeplitz matrix made of the coefficients of the polynomial g: 

\ 



Tig) = 





51 


92 


53 


54 





50 


51 


52 


53 








50 


51 


52 











50 


51 














50 


\: 











(1.128) 



Here we set gi := for all i > I. 

®That is, Aj(p, g) is the leading principal minor of the matrix 'H2n{p,g) of order j if degg < degp. Otherwise (when 
degg = degp), Aj(p,g) denotes the leading principal minor of the matrix 'H2n+i (p, 9) oforder jr. 
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Proof. Straightforward multiplication of matrices H{p,q) and T{g). 

Denote by rij{p, q) the leading principal minor of the matrix H{p, q) of order j (j = 1,2,.. .). We now 
derive some key connections between these minors and the minors Dj, Dj and V2j we encountered before. 

Lemma 1.44. Let the polynomials p and q be defined by (|1.2p - (|1.3p and let 

„/ ^ q(z) Sn Si 

R{z) = ^ = s_i + ^ + 4 + . . . 

p[z) Z Z'^ 

The following relations hold between the determinants 7]j{p, q) and Dj{R), Dj{R), V2j(p, q) defined by (jl.6l) . 
()1.9|) and ()1.12|) . respectively. 



If deg q < deg i/ien 



J = l,2,.. 



772,+ib,g) = aoV2,(zg,p) = (-lpa^^+'5,(i?), j = 0, 1,2, 
If deg q = degp, then 



V2j+i{p,q) = ba^2j{p,q) = baal^Dj{R) 



i = 0,1,2,. 



(1.129) 
(1.130) 

(1.131) 
(1.132) 



V2j{p,q) = aQbQW2j-2{h,p) ^ {-ly ^boal^ ^Dj^i{R), j^l,2,..., 

where h[z) := zq{z) -zp{z) and Dq[R) :— Dq{R) := 1. 

ao 

Proof. First, we prove the more complicated equalities (|1.13ip - (|1.132p . The formula (|1.13ip follows 
from (|1.124p . pTT^ and ([TTTO)) . To prove (|1.132p . we consider the function 



G(z) :— zR{z) — s_iz = 



Hz) 

p{z) 



Si 



S2 



= So H H ^ 



Here we used the fact that s_i — —. From (|1.6p and (|1.9p it follows that 

ao 



D,{G)=D,{R), 
Next, for a fixed index j — 1,2,..., we have 



J = 1,2,. 



(1.133) 



(1.134) 



V2jip,q) 



bo 


bi 


62 . 


• fe.-l 


b, . 


■ b2j-i 





ao 


ai . 




aj-i ■ 


■ a2j-2 





bo 


bi . 




bj-i ■ 


■ b2j-2 








ao ■ 


• 


O'j-2 ■ 


■ a2j-3 








. 


ao 


ai 


aj 








. 


. bo 


bi . 


■ b, 








. 





ao 





For each i = 1, ... ,j — 1, we now subtract the (2z)th row multiplied by — from the (2i+l)st row to obtain 

ao 



(1.135) 



V2j{p,q) = aobo\72j-2{h,p) = ao^o (-!)■' ^^2j-2{p,h) . 

This proves the first part of (I1.132p . Now from (fTTSl) . (|1.133l) and (|1.134p we have 

^2AP,h)^cil'D,{G)=al'b,(R). 

Combined with (|1.135p . this implies (jl.l32p . The formula (|1.129p - (|1.130p can be proved analogously 
to (|1.13ip - (|1.132p . applying Theorem ll.Sl to the functions R{z) and zR{z) and using Definition 11.401 □ 



We next turn to finite matrices of Hurwitz type and consider some of their properties. 
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Theorem 1.45. Let the polynomials p and q be defined as in (|1.2l) - (jl.3p and let 



Ifdegq < degp, then 



R{z) = —— = s_i H h ^ 

p{z) z 



A2,_i(p,g) = al'-^Dj{R), j = l,2,...,n, (1.136) 

A2,b,g) = {-iyal'D,{R), j = l,2,...,n. (1.137) 



Ifdegq — degp, then 



Ayip,q) - al'D.iR), j = l,2,...,n, (1.138) 

A2,+i(p,g) = (-l)-'a^-'+'5,(i?), J =0,1, 2,..., n, (1.139) 

where Do{R):—l and the determinants Dj{R),Dj{R) are defined by (jl.6p . (jl.9p . respectively. 

Proof. Apply Theorem II. 51 to the functions R{z) and zR{z) if degq < degp or to the functions R{z) and 

zR{z) ~ z— if degg = degp, as in the proof of Lemma [1.441 □ 
ao 

We now summarize ah previous resuhs and add one more fact: 
Theorem 1.46. Given polynomials 

p{z) = aoz" + fliz"^-^ + • • • + a„, ai,...,a„GC, oq 7^ 0, 
q{z) = 6o^" + fei^""^ + ■■• + &«, 6o,...,&„eC, 



let the function 



T?! ^ 9(2) So si 

R[z) = = s-i H h ^ 

p(z) z 



have exactly r poles {r < n), counting multiplicities. The following conditions are equivalent: 

1) the Hankel determinants Dj{R) and Dj{R) defined in p. 61) - (|1.9p satisfy 

D,{R)^Q, j = l,2,...,r, 

D,{R)^Q, j = l,2,...,r-l, 

D,{R)^Dj{R)^Q, j=r + l,r + 2,...; 

moreover, Dr{R) — if and only if R has a pole at 0; 

2) the function R has a Stieltjes continued fraction expansion (jl.ll6p whose coefficients Cj can be found 
by the formula; p.ll4p - (|1.115p ; 

3) the Euclidean algorithm applied to the polynomials goiz) :=p(z^) anff^ 9iiz) '■= zq{z'^) -zp{z^) is 

ao 

regular and deg q > n — 1; 

4) the infinite matrix H{p, q) factors as follows: 

if deg q < deg p, then 

Hip, q) = J(ci) J(c2) • • • J(cfc)i/(0, l)r(ff), (1.140) 

if deg q = degp, then 

H{p, q) = J(co) J(ci) • • • J{ck)HiO, l)T{g), (1.141) 



-'If degq < degp, then bo = 0, according to our convention. 
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where k = 2r ~ I, if R{0) = oo, and k = 2r, otherwise. The polynomial g is the general common 
divisor of p and q, the matrix T{g) is defined in (jl.l28p and 



J(c):= 



/c 


1 















(I 














A 








1 



































c 


1 















1 






















1 




, H{0,1):= 































c 


















1 




V 










■■) 




V 










■■/ 



(1.142) 



The coefficients Cj are defined by the formulas (jl.ll4p - (jl.llSp or, equivalently, as follows: 
if deg q < deg p, then 

vlj-iiP, q) 



C2j-1 

if deg q — degp, then 

C2J-1 = 

C2j -- 



V2j~iip,q) ■ V2jip, q)' 



vljip^ q) 
mj-i{p, q) ■ V2j+i{p,q)' 
Vy+iip,q) 



J = 1,2,. 
1,2,. 



J 



1,2, 



i = 0,1,2,, 



mjip^q) ■ V2j+2{p,q)' 

Here we set ri(){p, q) :— 1 and k :— 2r — 1 if R{0) — oo, whereas k := 2r if |i?(0)| < oo. 
Proof. Theorem 11.391 establishes the equivalence of conditions 1) and 2) 



(1.143) 
(1.144) 

(1.145) 
(1.146) 



By Theorem ll.26[ condition 3) is equivalent to the fact that the function z !-> F{z) — 



go{z) 



has a 



J- fraction expansion. At the same time, this is equivalent to condition 1), according to Theorem II. 261 and 
Lemma [1.341 Thus, we have proved the equivalence of conditions 1) and 3). 

Now we prove that condition 4) follows from condition 3). As was shown above, whenever the Euclidean 
algorithm applied to the polynomials go and gi is regular, the polynomials fa — p and fi~q satisfy ()1.119p - 
(jl.l20l) with the coefficients Cj determined by the formulae (I1.114p - (ll.ll5p . From Lemma [1.441 we see that 
the coefficients Cj can be found by the formula (|1.143p - (|1.144p or (jl.l45l) - (|1.146p . 

At first, let degg < degp and let g — gcd{p,q). The algorithm (jl.ll9|) - (|1.120|) produces a sequence 

p{z) 



of polynomials /o, /i, . . ., /fc, where k — 2r if p(0) 7^ or fc = 2r— 1 otherwise. Note that fo{z) 



and fi{z) 



q(z) 



, that is, gcd(/o,/i) — 1. Consider any four consecutive polynomials /j-i, /j, /j+i, 



g{z) 

(j = 1, . . . , 2r — 3) in this sequence such that deg/j_i > deg fj, i.e., j is odd. Then the matrix H{fj-i,fj) 
(see Definition 1 1.40p satisfies 



Hifj-iJj) = J(c,)i?(/,+i,/,) = J(c,)J(c,+i)i?(/,+i,/,+2). 



(1.147) 



This formula can be easily obtained by straightforward calculation using (|1.119|) ~ (|1.121l) . For the polyno- 
mials /2r -2 and /2r-i, the formula (|1.147p has a different form: 



if 2r - 1, then 

if(/2r-2,/2r-l) = J{c2r-l)H {f2rj2r-l) = J {c2r-l)H [0 ,1) , 

since /2r-i = fk = gcd(/o, /i) = 1, and therefore /2r(z) = fk+i{z) = 0, 



(1.148) 
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if = 2r, then 



H{f2r-2, f2r~l) — J {c2r-l)H {f2r , f2r-l) = J {C2r-l) J {c2r)H {f2r+l , f2r) 
= J{c2r^l)J{c2r)H{0,l), 



(1.149) 



since f2r = fk = gcd(/o,/i) = 1, and /2r+i(z) = fk+i^z) = 0. 
Thus, from the formulas (|1.147p - (|1.149p we obtain 

Hifo, h) = J(ci) J(C2) • • • J(cfe_i) J(cfe)i7(0, 1). (1.150) 
At the same time, Theorem II .431 impHes 

H{p,q) = HifoJi)ng). (1.151) 

The formula? (pl5ni) - (frT5T|l yield (fTTIOll . 

As for the case degp = degg, we denote fo{z) := ^7^, fi{z) := ^^^^^ — '^oP(^) ^ .^^j^gj-g g^j^^j gj^^^ 

5(2;) g[z) ao 

by straightforward calculation that 

H{p, q) = J(co)i/(/o, /i)r(g). (1.152) 

Since deg/o > deg/i, the matrix i/(/o, /i) satisfies (|1.150p . Thus, from (|1.152p and (|1.150p we obtain the 
factorization ()1.14ip . 

Conversely, if condition 4) holds, then we can reconstruct the algorithm (I1.119p - ()1.120p using the 
factorizations (|1.140p or (|1.14ip as follows: the coefficients of the polynomials and fj (if deg/^^i > 
deg/j) are the entries in the first and the second rows, respectively, of the matrix 

H{fj-i,fj) = J(c,) J(c,+i) • • • J[ck-i)J[ck)H{Q, 1), J = 1, . . . , k. (1.153) 

Here we have /o(2;) = and Ji{z) — ^^^^ — cop(^) ^ ^]^g]-jTT| _ Note that 50 (-z) — /o(^^) and 

9{z) g[z) oo 

gi{z) — z/i(z^) in our notation. As was shown above, the algorithm (I1.119p - (|1.120p for the polynomials 
/o and /i is equivalent to a regular Euclidean algorithm for the polynomials go{z) and gi{z). Therefore, 
condition 4) implies condition 3). □ 

Remark 1.47. Regarding equivalence classes of rational functions, we should note the following: Suppose 
that two rational functions R and G with exactly r poles each satisfy the inequalities (|1.106p - (|1.107p . 
Then R{z) = G{z) if and only if 



D,{R)^D,{G), 
D,{R)=D,{G), 



1,2 



; ^1 ■ 



since these equalities guarantee that the corresponding Stieltjes coefficients of R and G coincide by (|1.114p - 
(II. lisp . We would like to remind the reader that, according to Theorem 11.201 the equality of the minors 
Dj{R) and Dj{G), for each j, per se does not guarantee that the functions R and G are equal. 



2 Real rational functions and related topics 

In this section we develop connections among several notions: the Euclidean algorithm and its variant, the 
Sturm algorithm (Section 12. II) . Cauchy indices fSection l2.2p . various representations of rational functions, 
and their associated Hankel minors. Those diverse topics turn out to be connected to the same basic 
question of counting roots or poles of rational functions. Throughout this section, we assume that all our 
rational functions are real. 



^^As was mentioned above, feo = if degp > degg. 
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Thus, consider a real rational function 

i?(z) = 44 = «-i + - + ^ + ^ + --- ' z = -l,0,l,2,.... (2.1) 

p{z) Z Z'^ z'' 

where p and q are real polynomials 

:— uqz^^ + fliz"^^ + • • ■ + a„, ao, ai, . . . , a„ € M, ag 7^ 0, (2.2) 
q{z) := 6o^" + biz^'-^ + ■■■ + 6„, 60, • • • , &n e (2.3) 

In what follows we assume that the function R has exactly r = n — I poles (counting multiplicities) where 
I is the degree of the greatest common divisor of the polynomials p and q {0 < I < n). 



2.1 Sturm algorithm and Probenius rule of signs for Hankel minors 

For real polynomials /o :=p and /i := q, it is convenient to use a modification of the Euclidean algorithm, 
namely, its variant known as the Sturm algorithm. Suppose that we run the Euclidean algorithm (|1.42p 
starting with the polynomials /o and /i. If we denote 

f,{z) :=(-l)^/,(z), j=0,l,...,fc, 
Qiiz) ■■= {-lyqjiz), j = l,...,fc, 

then the polynomials fj and qk satisfy the following relations: 

f,^i{z) ^ q,{z)f,{z)~f,+i{z), j=0,l,...,fc, (2.4) 

where fk+i{z) = 0. 

Definition 2.1. The relations (|2.4p represent the so-called Sturm algorithm. If all the polynomials qj are 
linear, the algorithm is called regular. 

The polynomial fk is the greatest common divisor of the polynomials /o and /i. Thus, the Sturm 
algorithm also produces the greatest common divisor of two initial polynomials, but the Sturm form turns 
out to have advantages over the Euclidean form in the real case, as we will clarify later. Very roughly, in 
the real case signs of some quantities are more easily traced using the Sturm algorithm than the Euclidean 
algorithm. In the complex case, the issue of signs has no comparable significance. 

In connection with the Sturm method, we mention briefiy the so-called Sturm sequences, which, how- 
ever, will not be used much in the sequel: 

Definition 2.2. A sequence of polynomials go, gi, ■ ■ gn is called a Sturm sequence on the interval (a, b) 
if 

1) 5o(^*) — for some G (a, 6) gi{z*) ^ 0; 

2) gj{z^,) = for some € (a, b) (7j_i(z*)(7j+i(z*) < for j = 1, 2 . . . , n — 1; 

3) gn{z)^0 Vze{a,b). 

The sequence /o, /i, . . ., fr from the Sturm algorithm (|2.4p is easily seen to be a Sturm sequence on 
any interval where friz) does not vanish. Moreover, if we denote 

h^iz):=i^ j- = 0,l,...,r, 

fr{z) 

then hr{z) = 1 and the sequence hg, hi, . . hr is a Sturm sequence on the real axis. 

Next, we must introduce several notions of sign changes for sequences of real numbers. 
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Definition 2.3. Given a sequence t :={to, ti, . . . , tn) without zeros, we say that its number of sign changes 
is the number of indices j between 1 and n satisfying tj^itj < 0. The number of sign retentions is the 
number of indices j between 1 and n satisfying tj_itj > 0. 

For a sequence with zeros, the maximum number of sign changes obtainable by an appropriate choice of 
signs of any zero entry is called the number of weak sign changes and is denoted by V^(t) = Y~^{to, . . . , t„). 
The minimum number so obtainable is called the number of strong sign changes and is denoted by (t) = 
Y~ {to, . . . , tn)- The number of weak P~'^(t) and strong P (t) sign retentions can be defined correspondingly. 

Note that the number of weak sign changes does not increase and the number of strong sign changes 
does not decrease under small perturbations of the elements of a given sequence. Also note that the number 
of strong sign changes can be determined simply by discarding all zero elements and counting the number 
of ordinary sign changes in the obtained sequence. Finally, note that the number of sign changes and the 
number of sign retentions (be it ordinary, weak, or strong, respectively) always sum up to n if n + 1 is the 
length of the sequence: 

Y^{to,ti, . . . ,tn) +P^{to,ti, . . . ,tn) = n. (2.5) 

In the sequel, we will need only the notion of strong sign changes, so our discussion of weak sign 
changes above is included for completeness only. We will also need another important method of counting 
sign changes (and sign retentions) specifically introduced by Frobenius for sequences of Hankel minora^n 
{Do{R), Di{R),D2{R), . . . , Dr{R)) and {D^{R), Di{R),D2{R), . . . , Dr{R)) for a real rational function R. 
As usual, we set Do{R) ■.— Dq{R) 1. 

Rule 2.4 (Frobenius [23 BS])- If, for some integers i and j (0 < i < j), 

D,iR)^0, A+i(i?) = A+2(i?.) = --- = A+j(i?)=0, A+j+i(i?)^0, (2.6) 

then the number V'' {Do{R), Di{R), D2iR), DriR)) of Frobenius sign changes should be calculated by 
assigning signs as follows: 

signA+.(i?) - (-l)^^^signA(i?), i^ = l,2,...,j. (2.7) 

The number of Frobenius sign retentions {Dq{R) , Di{R), D2{R), • . . , Dr{R)) is defined accordingly. 

This assignment has an interesting property, which will be useful later. 

Corollary 2.5. // the condition (j2.6|) holds and if integers fi and v {i < fJ. < v < j) have equal parities, 
then 

Y''{D^{R),D^+,{R), D,{R)) = {D^{R), D^+,{R), . . . , D,{R)). (2.8) 

The Frobenius method of counting sign changes is so pervasive in the rest of this paper that we adopt 
the notational convention 

V(t):=V^(t), P(t):=P^(t). 



2.2 Cauchy indices and their properties 

In this section, we introduce a special counter known as the Cauchy index. Consider a real rational function 
F, which may in principle have a pole at oo, unlike the function R. 



Definition 2.6. The quantity 



Ind.(F):= +^ F(.-0)<0<F(. + 0), 

^ ' [-1 if F(cj-O) > > F(u; + 0), ^ ' 



is called the index of the function F at its real pole uj of odd order. 



Since our function R is real, all its minors Dj{R) and Dj{R) and therefore the sequences 
{Dq{R), D-i_{R), D2{R), ■ ■ ■ , Dr{R)) and {E)o{R),Di{R),D2{R),...,Dr{R)) are real. 
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We also set 

Ind„(F):=0 (2.10) 

if w is a real pole of the function F of even order. 

Suppose that the function F has m real poles in total, viz., uji < uj2 < ■ ■ ■ < 

Definition 2.7. The quantity 

Ind::(F):= Ind^.(J^)- (2-11) 

i : a<uJi<b 

is called the Cauchy index of the function F on the interval (a, b). 

We are primarily interested in the quantity lndt^{F) , the Cauchy index of F on the real line. However, 
since the function F may have a pole at the point oo, it is convenient for us to consider this pole as real. 
From this point of view, let us introduce the index at oo following, e.g., To do so, we consider the 

function as a map on the projective line PR^ :=R^ U {oo} into itself. So, if the function F has a pole at 
oo, then we let 

{+1 if F(+oo) < < F(-oo), 
-1 if F(+oo) > > F(-oo), (2.12) 
if sign F(+oo) = sign F(—oo). 

Thus, the generalized Cauchy index of the function F on the projective real line is 

IndpK(F) :=Indt^(F) +Indoo(F). (2.13) 

Remark 2.8. Obviously, if the function F has no pole at oo, then the generalized Cauchy index IndpR(-F') 
coincides with the usual Cauchy index Indl^(F). 

Following [11] , we list a few properties of generalized Cauchy indices, which will be of use later. 
First, note that a polynomial q{z) — cz'^ + ■ ■ ■ {v — degq) can be viewed as a rational function with a 
single pole at oo, hence 

T , / ^ |— signc if V is odd, , .x 

Indoo(g) = <^ / . (2.14) 

I U it IS even. 

The following theorem collects all properties of Cauchy indices that we will need. 

Theorem 2.9 (see HT). Let F he a real rational function. 

1) If d is a real constant, then IndpR((i + _F) = IndpR(i^). 

2) If q is a real polynomial and |i^(oo)| < oo, then IndpR((7 + F) = Indl|^(i^) + Indoo('z)- 

3) If G and F are real rational functions that have no real poles in common, then 

Indt^(F + G) = Indt^(F) + Indt-(G). (2.15) 

4) IndpR (^^^^ ^ IndpR(F). 

Proof. Properties 1), 2) and 3) follow immediately from the definition of Cauchy indices (j2.11|) - p.l2p and 
from (|2.14p . The proof of Property 4) is reproduced from [IT]. The projective line PM^ is divided by two 
points, and oo, into its positive (0,oo) and negative (— oo,0) rays. Suppose that a variable z traverses 
PM^ and returns to its starting point. Clearly, the number of crossings of F{z) from (— oo, 0) to (0, oo) 
must equal the number of reverse crossings. The crossings through oo occur at the poles of the function 
F; they are accounted for in the sum ()2.13p with the appropriate sign. The crossings through occur at 

the zeros of F, i.e., at the poles of the function z i— !■ ^ , and they are accounted for in the analogous 
formula for IndpR f ^ j • As a result, IndpR f ^ j + IndpR(F) =0. □ 
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Let us apply the Sturm algorithm (|2.4p to the numerator and denominator of the fraction R (see (|2.ip '). 
As a result, we obtain another kind of continued fraction, which slightly differs from (|1.46p : 

i?(z) = .s_i + (2.16) 

liiz) ; 



q2{z) - 



93(2) - 



Here the polynomials qi have the form 

gj(z) =a,z"^ +••• , aj^Q, j = l,2,...,fc, (2.17) 

wher^H ni + n2 + ■ ■ ■ + Uk = r and > 1, j = 1, . . . , /c. From Theorem II .181 it is easj0 to specialize the 
formula pTeS]) to the function ([2TT6)) : 

^IL ni(ni-l) _L. 1 

i?„,+„,+...+„^.(i?)-[](-l)^^ •[]—-— , j = l,2,...,fc. (2.18) 

Applying property 1) and then inductively properties 4), 2) and 3) of Theorem 12.91 to the function (I2.16p . 
we get the following result: 

Theorem 2.10 ([H]). If a rational function R is represented by a continued fraction (|2.16p . then 

k 

IndpR(i?) = -^Ind^(9j). (2.19) 
j=i 

This theorem implies an important fact recorded in Theorem 12 . 1 1 1 below . That fact is closely connected 
to the theory of quadratic forms. In fact, it was initially proved for rational functions with simple poles 
by Hermite using quadratic forms [IS] (see also [HI pp. 397-414]) and for arbitrary rational functions by 
Hurwitz "351 (see also [BT, HI])- Our proof differs from those proofs as well as from the proofs of 
Gantmachcr in 36 and Barkovsky in in that it does not use the theory of quadratic forms, only 
Frobenius Rule 12. 4[ some properties of continued fractions, and Theorem 12. 101 about Cauchy indices. 

Theorem 2.11 ([45 1 144 j ). // a rational function R with exactly r poles is represented by a series (|2.ip . 
then 

Indt^iR) = r -2V{Do{R), Di{R), D2{R), . . . , Dr{R)). (2.20) 
where the determinants Dj{R) are defined in (jl.6p and where Dq{R) 



Proof. According to (j2.5p . the assertion of the theorem is equivalent to the following formula 

Indt^(i?) = V{Do{R),D^{R), D-2{R), . . . , DAR)) - YiDoiR), D,{R), D-^iR), . . . , DriR)). (2.21) 

Moreover, according to Remark we have IndpR(i?) = Ind^^(i?). 

First, let us assume that the function R has a J-fraction expansion obtained via a regular Sturm 
algorithm (see Definition 12. ip . that is, where all the polynomials (j2.17p are linear and k = r. Thus, Ui = 1 
(« = 1, 2, . . . , r), and the formula (I2.18P takes a simpler form (cf. (|1.8ip ) 

Dm^il 2,-2.+! ' J = l,2,...,r. (2.22) 

z=l 



^^Recall that r is the number of poles, counted with multiplicities, of the function R. 

'^^One should replace G by —G in the formula 1 11.491 1 and apply the formula 111.521 1 inductively to the function 1 12.161 1. taking 
into account that (-l)JDj(G) = Dj{-G). 
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This implies the foUowing formula: 

i-i 

a 



In our (regular) case, from (|2.14p and (|2.19p we obtain 



a,^ i, .1. n— ' J = l,2,...,r. (2.23) 



Indi^(i?) = ^signa,. (2.24) 

i=i 



But ((2?23| implies 



i?,_i(i?) 



signa, = sign-^-^ = P(2?,_i (i?), I?, (i?)) - Y{Dj^i{R), D,{R)), j - 1,2, . . (2.25) 

Combining ^^14^ and dSSS]), we obtain ([23T|) . 

Now let all polynomials in (|2.16p be of odd degrees, that is, let all the numbers nj be odd. Then p.l4p 
and (|2.19p imply 

k 

Indt^(i?) = ^signa,. (2.26) 



Using notation (|1.68p from (|2.18p . we see that 



«j =(-1) 2 ■ ^ -11^^ j = l,2,...,fc. (2.27) 

Take into account that the numbers m.j-i and rrij — 1 have equal parities for every j, 1 < j < fc, since 
their difference nj — 1 is even. Hence, according to Corollarv l2.5l fsee (|2.8p '). we have 

P(An,_i(i?), i?m,_i+i(i?), . . . , D,n,-i{R)) - Y {D„,,_,{R) , D„,^_,+i{R) , . . . , Z?„^_i(i?)) = 0. (2.28) 
The Frobenius Rule (I2.7P gives 

K-l)("j-2) 

signAn,-i(-R) = signi:)„^_i+„^_i(i?) = (-1) 2 signZJ^^. (i?), j = l,2,...,fc. 
And from ()2.27p we obtain 

D^._,{R) 

signckj = signa/ = (-1) 2 sign— — — — = 

n^in^ K-l)K-2) . An,-l(i?) . An,-l(fl) (2.29) 

= ( — 1) 2 •( — 1) 2 . sign — ; — — sign — - — — — 

^ ' ^ ' ^ D„,,{R) ^ D,n,{R) 

= P(i?„^._i(i?),An,(i?))-V(An,-l(i?),i?™,(i?)), J = l,2,...,fc. 

Here we used the condition that all (j = 1, 2, . . . , k) are odd. From (|2.28p and (j2.29p we obtain 
signaj = P(i?™^._,(i?),i^™^_,+i(i?),...,An,-i(i?),An,(i?)) - 

(2.30) 

- Y{Dra,_, (R), A„,_i+l(i?), . . . , -Dm,-l(i?), D„,^ (R)) , j = 1, 2, . . . , /c. 

Substituting this formula into (|Z!26)) yields ((Z!2T|) . 

Let us now consider the general case and let 1 < ji < j2 < • • • < jri < k he the indices of those 
polynomials q^^, qi^, . . ., qi^ in (|2.16l) that have odd degrees. Then (|2.14p . (I2.17P and (|2.19p imply 

Indt^(i?)=^signa,,. (2.31) 
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Let i and j [1 < i < j < k) be integers such that and Uj are odd and ni_|_i, 71^+2, . . . , are all even. 
Then the numbers m^, to^+i, . . . , rnj^i, rrij — 1 (see (|1.68p ) have equal parities and Corollary 12.51 yields 

P(An.(i?), . . .,Dm,-iiR)) - V(A„,(i?), . . . = 0. (2.32) 

On the other hand, the formula (12.29^ is yalid in the present case, since rij is odd by assumption. Therefore, 
signaj = P(i?„.(i?),An.+i(i?),...,i?m,-i(i?),An,(i?)) - 

(2.33) 

- V(i?„.(i?), i^™,+i(i?), . . . , D^^.i{R),Dm^ (i?)). 
Applied to the indices ii, 12, . . ., ir] , the formula (I2.33P yields 

sign a,, - P(A„,,_,(i?),i?m,,_,+i(i?),...,i?m,,(i?)) 

(2.34) 

- YiD^,^_^ (R), D^^^_^+iiR), Dra,^ [R)), ^ = 1, 2, . . . , 77. 

If in < fc, then the indices n^-^+i, rij^4_2, . ■ ■ ,nk are all even, so all indices mj^ , mj^+i, mj^+2, . . . , irik haye 
equal parities, hence Corollary 12.51 implies 

P{D^^^ (i?), An,„+i(i?), . . . , An. (R)) - V(An,„ (i?), (R), . . . , An. (i?)) = 0. (2.35) 

The formulae ([234)) - (l235l) with (|23T|) yield ([2^, as desired. □ 

Now we show how Theorem 12 . 1 II can help in calculating the Cauchy index of a rational function if that 
function has a J-fraction or Stieltjes continued fraction expansion. 



Theorem 2.12. // a real rational function ()2.ip with exactly r poles has a J-fraction expansion 

1 



R{z) = s_ 



1 

aiz + /?i 



1 

a2Z + ^2 Y 

asz + (^3 



1 



arZ + j3r 

and m (0 < m < r) is the number of negative coefficients otj (j = 1, 2, . . . , r), then 

Indl:;^(i?) = r-2TO. (2.36) 

Proof. From (|2.25p we get that is negative if and only if there is a sign change in the se- 
quence {Dj_i{R),D^{R)}. Consequently, m = V(A(^), A(i?), A(i?), • ■ • , A(i?)), and ([^35)1 follows 
from (11201) • Also (^3^)) follows from ^^TI^ . □ 



Theorem 2.13. If a real rational function (|2.ip with exactly r poles has a Stieltjes continued fraction 
expansion (|1.116p and if ni (0 < m < r) is a number of negative coefficients C2j-i (j ~ 1, 2, . . . , r), then 
the Cauchy index Indlj^(_R) can be found by the formula (|2.36l) . 

Proof. The formula (|236l) follows from (|1.115p and (IQ}]) . □ 

Using Theorem 12. IH we can also produce formula; for calculating the Cauchy index of a rational 
function on the intervals (— oo,0) and (0, A-oo). 

Theorem 2.14. Let a real rational function R with exactly r poles have a series expansion (j2.1l) . Then the 
indices of the function R on the intervals (0, +00) and (—00, 0) can be found from the following formulae: 
• if |i?(0)| < 00, then 

Ind+°°(i?) = r- [v(l,A(i?),A(i?),...,A(i?))+V(l,A(i?),A(i?),...,A(i?))] , (2.37) 
Ind*L^(i?) = V(l,A(i?),A(i?),...,A(i?))-V(l,A(i?),A(i?),...,A(i?)), (2.38) 
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if z = is a pole of R of order v, then 

l+(Tl 



lndP_^{R) = i— ^ + V(l, D2{R), . . . , Dr-i{R)) - Di{R), D^iR), Dr{R)), 



V(l, Di{R), D2{R), Dr{R)) + V(l, Di{R),D2{R), . . . , Dr-i{R))\ , (2.39) 

(2.40) 



where ai = sign I lim R{z) ) , and 



ax if Indo R{z) 7^ 0, 
— CTi if Indo Ri/^) = 0. 

Proof. At first, let the function R have no pole at 0, that is, suppose that |i?(0)| < 00. Represent R{z) 
in the form 

R(z) = i?(i)(z) + i?(2)(z) + Giz), (2.41) 
where the function G has no real poles of odd order, that is, where Indl|^(G) = 0, and 



A 



A 



U) 
U-l 



A'i 



U) 



■^\{z+ U/jYi [z + CJj)'j 



U-l 



Z + UJj 



i?(^)(z) = 



A 



U) 



(2.42) 



A 



if) 



where all Ij are odd and Wj > 0. Here m_ ( m+) is the numbers of negative (positive) poles of odd order 
of the function R. It follows from Definitions 12.61 and 12.71 and from the formula; (|2.4ip - (|2.42p that 



Ind°_^(i?) ^ Ind°_^(i?«) = EsignKf), 

m_+m+ , , 

IndJ°°(i?) = Ind^°°(i?(2)) = ^ sign(^|^M. 
Now consider the function z H' F{z) — zR{z). From (12. ip and (I2.4ip we have 



(2.43) 



F{z) = zR^^\z) + zR^'^\z) + zG(z) = s^^iz 



So 



(2.44) 



where 



(i) 



+ 



^,A'(' - A 

Z + LOj 



i=i 



0) 



4(3) 



0) _ 4 0') \ 

m_ 



(2.45) 



w., Ar' - A 



From Definitions 12.61 and 12.71 and from the formulae (I2.43p - (|2.45p we obtain 



0) 



Ind"_^(i^) = Ind'i^(zi?(i)) = -^sign(a.,4f) = -^sign(A|f) = -Ind°_(i?), (2.46) 



Ind+-(F) ^ Ind+-(zi?(2)) ^ ^ sign(^,4f) = ^ sign(A«) = Indo+°°(i?), (2.47) 
since all ujj are positive. 



j—ni- +1 
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Theorems and and the formulaae (PTiSt - d^TTl) yield 

Indt^iR) = Ind"_^{R)+lnd+°^{R) = r-2Y{l,DiiR),D2iR),...,DriR)), 
Indl ^ (F) = lnd°_^{F) + Ind+- (F) = - Ind^ ^ (i?) + Indo+- (R) - 

= r - 2 V(l, i?2(i^), . . . , = r - 2 V(l, 52(i?), . . . , 

These formulae imply (P371) - (I05|) . 

If the function i? has a pole of order v (> 1) at 0, that is, if -R(O) = oo, then instead of (|2.4ip and (|2.44p 
we obtain 

R{z) = i?(i)(z)+i?(2)(z) + ^ + G(z), 

F{z) = zi?(i)(z) + zi?(2)(^) + _£_ + ;2G(z), 

where C = lim z^Riz), the functions R^^^ and i?*^) are the same as in (Em, and Indl;^(G') = 0. The 

formulae p.46p - p.47p remain the same. 

If is odd, then Indo(i?) = ai ^ but Indo(F) = 0, Dr{F) = Dr{R) = (see Corollary [Ol) , and 

Indt^(i?) ^ Ind"^(i?)+Ind+°°(i?)+(Ti =r-2V(l,i?i(i?),I?2(i?),...,i?r(i?)), 

Indt-(^^) - Ind" ^(F) + Ind+-(F) - - Ind"_^(i?) + Ind+-(i?) = 

= r-l-2V{l,D^{F),D2{F),...,Dr-i{F)) - 

= r-l - 2 V(l, Di{R),D2{R), . . . , Dr-i{R)). 



(2.48) 



(2.49) 



If is even, then Indo(i?) = but Indo(F) cri 7^ 0, Dr{F) = Dr{R) 0, and 

Indt^(i?) = lnd'L^{R)+lnd+°°{R) = r -2Y{1, Di{R),D2{R),...,Dr{R)), 

Indt^(F) = Ind'i^(F) + Ind+-(F) + a, = - Ind" ^(i?) + Ind+°^(i?) = 

- r-l-2V(l,i^i(F),i^2(F),...,i?r-i(F)) = 

= r-l - 2 V(l, 5i(i?), 52(i?), . . . , Dr-i{R)). 

The formulae (gigi-dSlHl) give ^X^-^L^. □ 

If a rational function R has a Stieltjes fraction expansion, then using Theorems 12.141 and for- 
mulae (|1.114p . we can establish relations between Ind" ^^(i?), Ind(J'°°(_R) and the signs of the coefficients in 
the Stieltjes fraction of R. 

Theorem 2.15. Suppose that a real rational function R with exactly r poles has a Stieltjes continued 
fraction expansion (|1.116p . Then the indices of the function R on the intervals (0, -|-c») and (— oo,0) can 
be found by the following formulas: 

if \R{0)\ < 00, then 

IndJ°°(i?) = n^-Uo, (2.50) 
Ind°^(i?) = r-K + 71,], (2.51) 

where Uq is the number of negative coefficients C2i-i, z = 1, 2, . . . , r, and Ug is the number of negative 
coefficients C2i, i — 1, 2, . . . , r. 

if z = is a pole of R of order v, then 



Ind+°°(i?) = i-^+n,-n„, (2.52) 
Ind"_^(i?) = r-i±^-K + n,], (2.53) 
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where is the number of negative coefficients C2i_i, i = l,2,...,r, is the number of negative 

' 6i if lndo(i?)^0, 



coefficients C2i, i — 1,2, . . . ,r ^ 1, 5i = sign ( lim z'^R{z) ) , and (5; 



2 — 



if Indo(i?) = 0. 



Proof. From the formula (|1.114p - (|1.115p we obtain 

no = Y{l,Di{R),D2iR),...,DriR)), 

n,^Y{l,Di{R),D2{R),...,Dk{R)), 

where fc = r if |-R(0)| < oo, and fc = r — 1 if R{0) — oo. The assertion of the theorem follows from these 
formula and from Theorem 12 .141 □ 

3 Rational functions mapping the upper half-plane to the lower 
half-plane 

3.1 General theory 

Now we specialize general properties of complex and real rational functions stated in the previous sections 
to the following very important class of functions: 

Definition 3.1. A rational function R is called an R-function of negative type (respectively, positive type) 
if it maps the upper half-plane of the complex plane to the lower half-plane (respectively, to itself) : 

Imz > =4> Imi?(z) < — negative type; 

Imz > Imi?(z) > — positive type. 

The name R-function appears first in the works of M.G. Krein and his progeny in connection with the 
theory of Stieltjes string and Stieltjes continued fractions (see, for example, [HH [52]). Below we discuss 
several well-known and some new relationships between i?-functions and continued fractions of Stieltjes 
type (see Definition ll.37p . By now, these functions, as well as their meromorphic analogues, have been 
considered by many authors and have acquired various names. For instance, these functions are called 
strongly real functions in the monograph |90j due to their property to take real values only for real values of 
the argument (more general and detailed consideration can be found in [19], see also Theorem 13.41 below) . 

Remark 3.2. In the sequel, we will deal with i?-functions of negative type only. But if an i?-function F 
is of negative type, then the function ~F is evidently an _/?-function of positive type. Hence all results 
obtained for i?-functions of negative type can be easy reformulated for i?-functions of positive type. 

At first, let us prove one necessary condition for a rational function to be an i?-function. 
Lemma 3.3. // a real rational function 

p{z) 

is an R-function of negative type, where p and q are real polynomials, then 

|degp-deg(j| < 1, (3.1) 

and 

R{z) = -az + (i + R{z), a>0, /3 € R, R{oo) = Q. (3.2) 
If a = 0, then R is an R-function of negative type. 
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Proof. Indeed, 

R{z) = h{z) + R{z), R{oo)=0, 

where h{z) = cz^ + ■ ■ ■ is a real polynomial of some degree j. Obviously, for 2 G C such that Imz > and 
sufficiently large, we have 

sign (Im i?(2:)) — sign (Im /i(z)) = sign{c sm{j ip)) , 

where ip = argz. If j > 2, then Im_R(z) takes both positive and negative values for some z from the upper 

complex half-plane, say, for argz = — and argz — — . Therefore, the degree of h is necessarily at most 

1: h{z) — ~az + /3, where a,/3 G M. But Imh(z) = — almz, thus, a must be nonnegative. Moreover, if 
a — 0, then lmR{z) = Imi?(z) < for z such that Imz > , and R{z) is an i?-function of negative type. 

Thus, if deg q > degp and the function i? = — is an i?-function of negative type, then 

P 

deg(7<degp+l. (3.3) 

q p 
Let degq < degp and let the function i? = - be an i?-function of negative type. Then the function — is 

p q 

an i?-function of negative type too and, therefore, 

degq < degp + 1. (3.4) 
From (ISa-IHa we obtain (O). □ 



The following theorem lists the most important properties of i?-functions. Parts of this theorem can 
be found in [75l |6ll [IS ESI [HI El El [90] . 

Theorem 3.4. Let p and q be real and coprim^^ polynomials satisfying p.ip . For the real rational 
function 

p{z) 

with exactly n — degp poles, the following conditions are equivalent: 

1) R is an R-function of negative type: 

Imz > ^ Imi?(z) < 0; (3.5) 

2) The function R can he represented in the form 



i?(z) = -az + /5 + ^ — a > 0, /? e R, e M, j = l,2, (3.6) 
j=i J 

where 



7. = #T>0, J = (3.7) 

3) The index of the function R is maximal: 

IndpR(i?) = max (degp, deg q) ; (3.8) 

4) The function R has a J -fraction expansion: 

R{z)^-az + /3 + , (3.9) 

aiz + I3i 

Q!2Z + (32 



asz + 

where a, > 0, g M, a > and (3 g M; 



1 



^^This condition is introduced for simplicity and just means that the number of poles of the function R equals the number 
of zeros of the polynomial p. 
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5) The polynomials p and q have only real roots and satisfy the inequality 

p{uj)q' {uj) - p' {uj)q{uj) < Q for all a; G R; (3.10) 

6) The roots of the polynomials p and q are real, simple and interlacing, that is, between any two 
consecutive roots of one of the polynomials there is exactly one root of the other polynomial, and 

3 w e M : p{uj)q{uj) - p{uj)q{uj) < 0; (3.11) 

7) The polynomial 

g{z) ^ \p{z)+ ^iq{z), (3.12) 
has only real zeros for any real A and /i, + /i^ 7^ 0, and the condition p.lip is satisfied; 

8) The function R{z) has real values only for real z: 

R(z) e M =^ z e M, (3.13) 

and 

3we]R: -00 < i?'(a;) < 0; (3.14) 

9) Let the function R be represented by the series 

R{z) = -az + l3+^ + ^ + ^ + --- (3.15) 

with a > and /? G M. The following inequalities hold 

D,{R)>0, i = l,2,...,n, (3.16) 
where the determinants Dj{R) are defined in (jl.6p . 

The equivalence of conditions 1) and 2) is usually called the Chebotarev theorem [97l[T9]. The equiva- 
lence of conditions 1) and 9) in case of meromorphic functions is the famous Grommer theorem [5^ HI [T5]. 
The equivalence between 2) and 4) was also proved by Grommer [3S| (see also [1011 Finally, the 
equivalence of 1) and 6) is a modification of the famous Hermite-Biehler theorem (for example, [5T1 

Proof. The scheme of our proof is as follows: 

1) =^ 2) =^ 3) ^ 4) ^ 9) 

8) ^ 7) ^ 6) ^ 5) 

1) => 2) Let the function R satisfy (|3.5p . Lemma 13.31 guarantees that R can be represented in the 
form p.2p . Thus, we have to prove that the function R in the representation (|3.2p has the following form 



wherj^^l q{z) = q{z) — {az + /3)p{z) and 



i?(z) = ^=^^^, u;,€R, j = l,2,...,n, (3.17) 



q(uji) q(uj-i) 

= ^TT^ ^ ^TT^ > 0' ^■ = 1.---"- 3-18 

p'{uj) p'iujj) 

At first, assume that the function R and, therefore, the function R have a nonreal pole A of some 
multiplicity j(> 1). Then for z = A + e, |e| — )■ 0, we have 

sign (Im i?(z)) — sign ^Im R{z)j = sign ^Ini — sign (sin(arg c — j ■ arge)) . 



^^Here a and are the numbers from the representation 1 13. 2| I. 
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From these equalities one can see that we can obviously choose such complex numbers ei and Si that 
zi = \ + £i and Z2 — \ + 62 are from the upper half plane, but sign (Imi?(zi)) = — sign (Im i?(2;2)). It 
contradicts with (13.51) . Thus, the functions R and R have only real poles. 

Now let us assume that R and R have a real pole ^ of some multiplicity j > 2. And let z — fi + e. 
Since /i is a real number, we have Imz = Ime. If we take z from the upper half-plane of complex plane 
and sufficiently close to /z, that is, < arge < n and |e| — ?> 0, then 

sign (Im R{z)) = sign ^Im R{z)^ = sign ^Im — = — sign {A sin(j • arg e)) . 

These equalities show that if we choose Si and £2 such that argei = — and arg £2 = — , then zi = /i + ei 

and Z2 = fJ, + £2 both belong to the upper half-plane, if j > 0, but sign (Imi?(zi)) = — sign (Imi?(z2)). 
Thus, the function R (and the function R) has only simple real poles, and R satisfies p.l7p . where the 
residues jj are determined by the formula 

7,= hm {R{z){z-a))^ lim '\ . ' ^ ^ J = l,...,n. 

For z sufficiently close to a pole LUj {j — 1,2, ... ^n), we have 

sign (Im R{z)) = sign (im R{z) \ ~ sign ( Im — — — j = — sign (7^ Im z) . 



These equalities combined with p.Sp yield the positivity of all jj. Consequently, the function R satis- 
fies ([3T7l) - (|3?T8ll . but R has the representation ((3?6l) - ([3T71) . 



2) => 3) Let the function R satisfy p.6p . If degp > degg, then a = and max(degp, degg) — degp = n. 
In this case, R has no pole at 00, so from (|2.9p we obtain 



IndpR(i?) = Indl!^(i?) = n = max(degp, degq). 
If degq — degp + 1, that is, a > 0, then R has a pole at 00, so according to (|2.9|) . and ()2.12|) - (j2.14p . 
IndpR(i?) — sign a + Indt|^(i?) = 1 + n — dcg q = max(degp, deg q). 

3) 4) Since the function R satisfies p.ip . it can be represented as follows: 

R{z)^-az + l3 + R{z), (3.19) 
where a, ^ G M and -R(oo) 0. Then from (|2?T2l) - (|2?T5l) we have 

IndpR(i?) = signa + Indt:^(i?) =signa + Indt^(^), (3.20) 
therefore, according to (|2.9p and (|2.12p . 
if deg q — degp + 1, then 

-n- 1 < IndpR(i?) < 71 + 1, (3.21) 

if degg < degp, then 

-n< IndpR(i?) < n, (3.22) 

The condition p.Sp is equivalent to 

T , /m jn + 1 if degq = degp +1; 

IndpK(i?) = <^ (3.23) 
I 71 II deg q < degp. 

From (IS2ni)-(E2Sl) we obtain that a > in ([XTO]) and Indt^(i?) = Ind+;^(^) = n. 
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Let us expand the function R into a continued fraction (I2.16p . Then Theorem 12.101 yields 

fe 

IndpR(i?) = -y^Indoo(gj), 

where k < n. Since IndpR(i?) must be equal to n, we see that k — n and that all polynomials qj 
(j = 1, 2, . . . , n) are linear with positive leading coefficients: 

qj{z)^ajz + (3j, > 0, /3j G R, j = l,2,...,n. (3.24) 

fc 

This follows from the fact that degq^ — n and from (I2.14p . 

Thus, if the function R satisfies (|3.8I) . then we obtain from (13.19^ (where a must be nonnegative) 
and p. 241) that R has a J- fraction expansion p.9p . 

4) => 5) If the function R{z) has a J- fraction expansion p.9|) . then 

i?(z) = 44 = + ^44: a>0> (3.25) 
H^) fo{z) 

fi 

where the function — has a J- fraction expansion (j3.9p and vanishes at cxi. Therefore, starting from the 
J" 

polynomials /q, /i, we can construct a sequence /o, /i, ■ • ■ , /n by the Sturm algorithm: 

/,_i(2) = (a,2 + /3,)/,(z)-/,+i(z), a, >0, /3j € R, i = 0,l,...,n, (3.26) 

where /„+i(z) = 0, and gcd(/o, /i) = fn{z) = 1, since the polynomials p and g (and, therefore, /o and /i) 
are coprime by assumption. From (|3.26p we have that, for all w S R, 

foHfiiu:) - f{,{io)h{oj) - -aif^ioj) + /i(c.)/^(c.) - /{(^)/2(c^) = 

(3.27) 

= -ai/2(c.) - a2.fi{u) + /2(c.)/^(c.) - f^{io)f3{oj) = = - EU «./'(^) < 

since all aj > and all /J (w) > and /^(w) > as well. 

Now from (I3.25|) we obtain that q{z) — {~az + l3)fo{z) + fi{z) and p = /q. Thus, p.27p implies 

p{uj)q'{u)-p'{uj)qiuj) = -a/2(a;) +/o(c.)/{(a;) - /^(c.)/i(a;) < 0, 

for any a; G R, as required. 

Let z be a complex number such that Imz 7^ 0. It is easy to see that 

sign ( Im , , ) = Q;sign(Imz) - sign(Im/(z)) , a > 0, /3 € R, 

where / is any function of a complex variable. Thus, if the function R has a J- fraction expansion (|3.9[) . 
then we obtain 



sign (Im R{z)) = sign ( Im J^J~^ ] ^ ~ sign 



Im 2(0; + ai 



^0, 



whenever Im z 7^ 0. Thus, R necessarily has only real zeros and poleS 

5) 6) If the polynomials p and q are real-rooted and satisfy the inequality p.lOp . then the condi- 
tion p. lip holds for them. From p.lOp it also immediately follows that all zeros of p and q are simple. 
Otherwise, there must be a real number oj such that p{uj)q'{uj) — p'{uj)q{uj) = 0, which contradicts p.lOp . 



''At the same time, this proves the imphcation 4) =^ 1). 
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Let real numbers Wi,aj2 (i^i < 1^2) be two consecutive (simple) zeros of p. By Rolle's theorem, we have 

p'{ui)p'{uj2) < 0. (3.28) 

Then from ((XTU)) it follows that p'{uJi)q{uji) < and ^'(^2)^(^2) < 0. Together with ([X^ . these 
inequalities imply q{iOi)q{uj2) < 0. Thus, in the interval (cji,aj2), the polynomial q has an odd number of 
(simple) zeros. In the same way, one can prove that between any consecutive zeros of q, there is an odd 
number of zeros of p. Therefore, all zeros of p and q are simple, and between any two consecutive zeros of 
one of the polynomials there is only one zero of the other polynomial, as required. 

6) =4> 7) Now let the polynomials p and q have simple real interlacing zeros. If numbers wi and uj2 
(wi < u)2) are some two consecutive zeros of the polynomial p, then by interlacing we have q(uJi)q{uj2) < 
and the polynomial g defined p.l2p satisfies the inequality g{uJi)g{uj2) < 0. Therefore, all zeros of the 
polynomial g are real and simple and interlace the zeros of p since | degp — degg\ < 1. 

7) => 8) Let the polynomial g defined by (|3.12[) have only real zeros and suppose, without loss of 
generality, that /i 7^ 0. Then the equation 

p{z) (a^^ + a) = (3.29) 

has only real solutions for any real A and ^ ^ 0. Therefore, the function R = — cannot take real values 

P 

for nonreal z. Otherwise, the equation p.29p would have nonreal solutions for some real A and fi. Thus, 

the function R = — satisfies (I3.13|) . 
P 

The condition p.lip implies (|3.14l) . Indeed, if a real u in p.lip is not a zero of p, then 

p// N P(^)g'(^) -y(t^)g(^) ^ 
U (cj) = J— < U. 

P (w) 

If ui from (I3.1ip is a zero of the polynomial p, then w is a pole of the function R'. However, (|3.1ip shows 
that, for sufficiently small e > 0, 

p{lo + e)q'{Lu + e) - p'{uj + e)q{uj + e) < 0, 

since the function pg' — p'q is continuous and preserves its sign in some vicinity of the point uj. Therefore, 
R'{uj + e) < 0, as required. 

8) 1) If the function R satisfies (|3.13l) . then it has no complex zeros, since is a real value. Thus, if 
Imz > 0, then Imi?(z) =/= 0, i.e., R is an i?-function of positive or negative type. Suppose that Imi?(z) > 

whenever Imz > 0, that is, R is an i?-function of positive type. Then the function F = —R = — is an 

P 

i?-function of negative type, i.e., it satisfies (|3.5p . But we have already proved that 1) 2) 3) 
4) => 5), therefore, we have 

p{uj)q'{uj) -p'{uj)q{uj) < for all G M. 

Consequently, 

for any real uj such that F'{uj) exists. This means that R'{uj) ~ — F'(cj) > for any real lo such that i?'(w) 
exists. This contradicts p.l4p . Therefore, R is an i?-function of negative type. 

4) -4=4' 9) Theorem ll .261 guarantees that the function R has a J-fraction expansion p. 91) with all dj nonzero 
and real, j = 1, 2, . . . , n, if and only if the Hankel minors Dj(R), j = 1, 2, . . . , n, are nonzero. But for 
J-fraction (|3.9p . the formula (|2.22p holds (see also (|2.23p ). This formula implies that the inequalities p.l6p 
hold if and only if all aj > for j = 1, 2, . . . , n. □ 
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Remark 3.5. Comparing representations p. 61) and (I3.15p . one can obtain the following formula: 

n 

= ^7ja;;., 0,1,2,... (3.30) 



Remark 3.6. From the proof of Theorem 13.41 one can see that a sum of two i?-functions of negative type 
is an i?-functions of negative type. Also if R{z) is an i?-fimction of negative type, then the functions 

z I— > — z) and z — are also i?-functions of negative type. 

R{z) 

Using the equivalence of conditions 1) and 7) of Theorem 13.41 one can obtain the following simple fact. 

Corollary 3.7. Let p and q be real coprime polynomials satisfying (|3.1I) . degp > 2. // the function 
R = q/p is an R-function of negative type, then the functions Rj = q^^^ jp'^^^ ,7 = 1,..., degp — 1, are also 
R-functions of negative type. 

Proof. It suffices to prove that the function R\ is an iZ-function. First, consider the case degq — degp— 1: 

n n 

p{z) = J2 = ^i^""^"^ ^here ao > 0, bi ^ 0, n > 2. 

3=0 j=l 

Then R{z) — >■ as z — > +oo, and, for sufficiently large positive z, signi?(z) = signg(z). Suppose that 
R is an i?-function of negative type. By Theorem 13.41 the polynomial g defined by p.l2p has only real 
zeros. Therefore, the polynomial g' also has only real zeros, so the function i?i = q' /p' is an iZ-function of 
positive or negative type according to Theorem 13.41 and Remark 13.21 

Since R is an iZ-function of negative type, R is decreasing between its poles (this follows, for example, 
from (|3.6|) '). Therefore, if a real number ^ is the largest zero of the polynomial p, then R must be positive 
in the interval (^, +c»), so R{z) — > +0 as z — >■ +oo. Consequently, the polynomial q{z) is positive for 
sufficiently large positive z, i.e., 6i > 0. Now it is easy to see that 

p' {z)q" {z) ~ p" {z)q' {z) n(n - l)aofeiz2""'* < as z ^ ±oo 

Hence there exists a real w such that p'{uj)q"{ijj) — p" {uj)q' [lo) < 0, so Ri is an i?-function of negative type 
by Theorem! 



Let now degq = degp. Then 

^W-44 = ^ + ^' h{z) = q(z) - Pp{z), deg/i<degp, (3.31) 
p(z) p(z) 

where /? e R and the function h/p is an iZ-function of negative type by Thcorcm l3.4l (see (|3.6p ). Therefore, 
the function h' /p' is also an iZ-function of negative type. But from p.3ip it follows that 

p'{z) p'{z) 

Thus, the function Ri is an i?-function of negative type by Theorem 13.41 

Finally, if degg — degp + 1, then we may apply the previous result to the functions F — ~p/q and 
Fi = -p'/q'. □ 

This theorem and Theorem 13.41 immediately imply the following result due to V.A. Markov (see [l9l 
Theorem 9, Chapter 1]). 

Theorem 3.8 (V.A. Markov). // the zeros of two real polynomials p and q are simple, real and interlacing, 
then the zeros of their derivatives p' and q' also are real, simple and interlacing. 

If a rational function R is an i?-function, then, using Theorem l2.141 one can easily find the numbers of 
negative and positive poles of this function. 
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Theorem 3.9. Let a rational function R with exactly r poles be an R-function of negative type and let R 
have a series expansion p.l5|) . Then the number r_ of negative poles of R eoualrn 

r_ = V(l,5i(i?),52(i?),...,5fe(i?)), (3.32) 

where k = r — I, if R(0) — oo, and k = r, if |i?(0)| < oo. The determinants Dj{R) are defined in (fTT 



Proof. In fact, Theorem [33] states that Y{1, Di{R), D2iR), . . . , DriR)) = if (and only if) R is an R- 
function (see p.lSp ). Moreover, all poles of R are real and simple, and all residues at those poles are 
positive (see ^^), therefore we get r_ = Ind°_^(i?). Thus, the formula (^^^ follows from (l05)) - (p:ig)) . 
where (72 = fi = 1- D 

It is convenient for us to consider separately the extreme cases of Theorem 13.91 

Corollary 3.10. Let a rational function R with exactly r poles be an R-function of negative type. All 
poles of R are negative if and only if 

5,_i(i?)5,(i?) <0, j = l,2,...,r, (3.33) 

where Dq{R) :=1, and the determinants Dj{R) are defined by (|1.9p . 

Proof. Indeed, Theorem 13.91 implies 

V(l, Di{R),D2{R), . . . , DriR)) = r. (3.34) 

According to the Frobenius rule (12.71) . if Dj{R) ~ for some j < r— 1 but Dj^i{R) ^ 0, then sign Z?j(i?) = 
signl)j_i(i?), that is, Y{Dj^i{R), Dj{R)) = 0. Therefore, if there are zero determinants Dj{R) in the 
sequence (1, Di{R), D2{R), ■ ■ ■ , Dr{R)), then the equality p. 341) cannot hold. Consequently, all minors 
Dj{R) for j = 1, 2, . . . , r are not equal to zero, and from (|3.34p we obtain p.33p . □ 

Remark 3.11. The inequalities p.33|) are equivalent to the following inequalities 

{-l)W,{R)>0, j = l,2,...,r, (3.35) 

Corollary 3.12. Let a rational function R with exactly r poles (counting multiplicities) be an R-function 
of negative type. All poles of R are positive if and only if 

D,{R)>0, j = l,2,...,r, (3.36) 
where the determinants Dj{R) are defined in (jl.9p . 

Proof. If the function R is an i?-function of negative type with positive poles and has a series expan- 
sion p.lSp . then the function 

F[z) = -zR[-z) = -az - pH o + ^ 1 + , 

z z z'^ z^ 

where a > and /3 € M, is an i?-function of negative type with negative poles. In fact, since the function 
R has the form (|3.6p with all positive ujj and 7j, then F can be represented as follows 



R{z)^-az-(3 + J2^^, «>0, Wj,7j>0, i = l,2,. 



'^^Recall that the number V(l, bi(R), D2{R), ■ ■ ■ , -Dfc(-R)) of Frobenius sign changes must be calculated according to Frobe- 
nius Rule [Ml 
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Therefore, according to Theoreni l3.4[ the function F is an i?-function of negative type with only negative 
poles. On the other hand, for a fixed integer j (j > 1) we have 



-si 

-S3 



S2 
-S3 
S4 



-S3 

Si 
-S5 



51 S2 

52 S3 



S3 
S4 

S3 Si S5 



Sj Sj + i Sj+2 



Sj+1 
Sj+2 

S2j-1 



1 

0-10 

1 





{-iy+'s,+2 









^10 0... 
1 ... 
0-1 ... 

-S2J-1 ... (-l)J 

■ D,iR) ■ i-l)^- 



= i-iyD,{R), 

Consequently, from these formula and the inequalities p.35p we obtain p.36p . 



3.2 Some classes of infinite Hankel matrices and the finite moment problem 
on the real axis 

Here we connect the characterization of i?-functions from Section 13.11 to a particular moment problem, 
viz., a discrete moment problem on the real line with a measure supported at finitely many points. This 
problem is quite well known (see, e.g., pl|91j): 

Problem 3.13 (Finite moment problem on M). Given an infinite sequence of real numbers 
it is required to determine number^^ 

71 > 0, 72 > 0, . . . , 7„ > 0, LOl < UJ2 < ■ . . < UJk <0 < Wfe+l < OJk+2 < ■■■ <UJ„ 

so that the equations (j3.30p hold: 

n 

= ^7jwi, z = 0,l,2,... (IS3D1) 
From Remark lXSl it follows that the equalities p.30p are equivalent to the following series representation 

n 

N 7j So Sl S2 S3 
F{z) := > — = — + — + — + — + ••• (3-37) 

^ — ' Z — CJo Z Z'^ z^ 

j=l J 

In this case, the infinite Hankel matrix S = ||si+j||g° is of finite rank n. Thus, our moment problem has 
a solution if and only if the function F determined by the series (j3.37p is an i?-function of negative type 
with r poles and with exactly k (< n) negative poles. The solution of the moment problem is unique, since 
the positive numbers jj and ujj (j — 1,2, . . . ,n) are uniquely determined from the expansion p.37p . We 
will see that Theorems 13.41 and 13.91 provide a solution to this problem in Theorem 13.201 below, with two 
important special cases provided in Theorem 13.181 and Corollary 13.191 However, before proceeding, we 
must introduce and discuss some relevant matrix notions of positivity/nonnegativity and sign regularity. 

Definition 3.14. An infinite symmetric matrix A of finite rank is called r-positive definite if all its 
principal minors up to order r(> 1) (inclusive) are positive. If the rank of the matrix A equals r and A is 
r-positive definite, then A is called positive definite. 



^^For A: = 0, it is understood that uii > 0. Analogously, if A: = n, we have LUn < 0. 
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Definition 3.15. An infinite matrix A of finite rank is called r-strictly totally positive if all its minors up 
to order r (inclusive) are positive. If the rank of the matrix A equals r and A is r-strictly totally positive, 
then A is called strictly totally positive. 

Definition 3.16. An infinite matrix is called totally nonnegative if all its minors are nonnegative. 

If A is a matrix (finite or infinite), then its minor of order j(> 1) whose rows are indexed by ii, 12, . . . , «j 
and whose columns are indexed by li,l2, ■ ■ ■ ,lj is denoted as 



I. 



Definition 3.17 ( |34j ) . An infinite matrix A of finite rank is called r-sign regular if all its minors up to 
order r (inclusive) satisfy the following inequalities: 



(^-l)k=i k=i A 



> 0. 



(3.38) 



I2 ... Ij^ 

If, in addition, the rank of the matrix A equals r, then A is called sign regular. 

With these definitions in place, we are now ready to state and to prove the following theorem. 
Theorem 3.18. A function 

is an R-function of negative type and has exactly r poles all of which are positive if and only if the matrix 



^j=o strictly totally positive of rank r. 



Proof. This proof follows the presentation from [36, p. 238]. 

At first, let us assume that the function R is an i?-function of negative type with r poles and all poles 
are positive. Then by Theorem 13.41 we have 



7,>0, J = 1,2,.. 



(3.40) 



where 



According to formula p.30p (with n 
represented as follows 



< Wi < W2 < • • ■ < l^r- (3-41) 

r), an arbitrary submatrix of the matrix S of order k (< r) can be 



^ik+jkj 



'l2 t2 



\7i^r 72^2" 



4' 



7r<V \wf 



Therefore, 



«1 «2 •■ 
Jl j2 ■ ■ 



Ik 
jk 



^ik+jk 



E 



7<Tl7<T2 ■■■I'Jk 



l<(Ti <(T2 <...<(Tfc <r 



(3.42) 



(3.43) 
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Both determinants here are generahzed Vandermonde determinants. Their positivity follows from p.4ip . 
as was proved in [33] (see also (SHI P-99 and p. 239] or [THl Part V, Chapter 1, Problem 48]). Consequently, 
any minor of S of order fc(< r) is positive, since all 7^8 are also positive. 

Thus, in our case, the matrix S is strictly totally positive and has rank r (the number of poles of the 
function i?), according to Theorem II .31 

Conversely, if the matrix S is strictly totally positive of rank r, then its leading principal minors Dj{R) 
are positive up to order r. Since the matrix S^^^ ~ ||si+i+i||r7=o is ^ submatrix of the matrix S*, it is also 
strictly totally positive and, in particular, its leading principal minors Dj(R) are positive up to order r. 
From Theorem 13.41 and Corollary 13.121 we obtain that the function R is an i?-function of negative type 
with exactly r poles, which are all positive. □ 

From this theorem it is easy to obtain the following corollary. 
Corollary 3.19. A function 



\ So Si S2 S3 

R{z) = — + — + — + — + ■■• 

z z-^ 

is an R-function of negative type and has exactly r poles, all of which are negative, if and only if the matrix 
S — \\si+j\\°°^Q is sign regular of rank r. 

Proof. The function R is an i?-function with only negative poles if and only if the function 



F{z)^-R{-z) = --- + ---, + --- = - 



tj_ 

r2 



t2 



r4 



(3.44) 



is an i?- function of negative type with only positive poles (with the same number of poles as the function R). 



Let T be the infinite matrix defined by T := \\t. 



^Q. Since tj = {—lysj {j = 1,2,.. .), we have 



T = ESE, 



(3.45) 



where the infinite matrix E has the form 



E:^ 



(I 
0-10 

10 

-1 

V ; : : : 



All minors of this matrix except for its principal minors are zero and 



(3.46) 



E 



H «2 



(-1)- , 



since E{j,j) = (—1)^ ^. From the Binet-Cauchy formula (see, e.g. 
T 



Ji h 



jk 



= E 



E E 

Tl<T2<...<Tm ei<e2<---<£ 

h 12 ■ ■ ■ ik 
h 12 ■ ■ ■ ik 



E 



s 



£1 £2 



s 



«1 «2 
Jl J2 



jk 



E 



]l J2 
jl J2 



, ,E'i + Ej!-2fc 

(-1)' = ! ' = 1 S 



n «2 

jl J2 



jk 
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This formula and Theorem 13.181 imply the assertion of the corollary. D 

Now we are in a position to formulate the solution to Moment Problem 13. 131 posed initially. Combined 
with the results in this section, Theorems 13.41 and 13.91 provide the following solution to this problem: 

Theorem 3.20. The infinite moment problem 

n 

St = ^7j^j, i== 0,1,2,... 

7i > 0, 72 > 0, . . . , 7„ > 0; wi < a;2 < • ■ ■ < i^fc < < t^k+i < < ■ ■ ■ < ujn, 

where Si, i = 0, 1, 2, . . ., are given real numbers and 7^ and ujj , j — 1,2, . . . ,n, are unknown real numbers, 
has a solution if and only if the infinite Hankel matrix S — ||si+j||f^=o ^'^^ rank n, the determinants Dj{S), 
j = 1,2, . . . ,n, defined in (jl.6p . are positive, aw(F"l 

k^V{l,Di{S),D2{S),...,Dr,{S)), 

where the determinants Dj{S) {j = 1,2, . . . , n) are defined by (|1.9|) . In that case, the solution is unique. 

3.3 i2-functions as ratios of polynomials 

In this section, we develop determinantal criteria for i?-functions involving the Hankel and Hurwitz minors 
formed from the coefficients of their numerator and denominator. 

Suppose that a rational function R is an i?-function of negative type, and write 

= + + + + (3.49) 

p[zj Z Z'^ Z'^ 

where p and q are real polynomials 

p{z) = aoz" + aiz"-^^ + ■ ■ ■ + On, ao, ai, . . . , a„ S M, qq > 0, (3.50) 
q{z) = 6oz" + + ■ ■ ■ + 6„, bo,bi,...,bn& M. (3.51) 

Since R is an i?-function, we know that degg > degp — 1, that is, b^ + bly^ 0. 

Remark 3.21. Generally speaking, the polynomials p and q may have a non-constant greatest common divi- 
sor g = gcd(p, q) of degree (n — r) for some natural r(< n). In this case, one should consider the function R 
as a ratio of polynomials p = p/g and q = q/g. Then the number of poles of R equals r. 

At first, we describe i?-functions in terms of the coefficients of the polynomials p and q. More precisely, 
we will use the infinite matrix of Hurwitz type H{p, q) (see Definition ll.40[) and the finite matrices of 
Hurwitz type (Definition II .421) constructed using the coefficients of polynomials p and q. 

Our first two theorems cover the case when degg < degp. 

Theorem 3.22. The function p.49p . where degq < degp, is an R-function of negative type with exactly 
r (< n) poles if and only if 

mj{P,q) > 0, j = l,2,...,r, (3.52) 
7j,{p,q) =0, i>2r + l, (3.53) 

where r)i{p,q) are the leading principal minors of the matrix H{p,q), defined by (jl.l23p . 

Theorem 3.23. The function p.49p . where degq < degp, is an R-function of negative type with exactly 
r {< n) poles if and only if 

A2,_i(p,g) > 0, j = l,2,...,r, (3.54) 
A,{p,q) = 0, i = 2r + l,2r + 2,...,27i, (3.55) 

where Ai{p,q) are the leading principal minors of the Hurwitz matrix 'H2n{p,Q) defined in (I1.125p . 

^"Recall that V is the number of Frobenius sign changes. The determinant D„(S) may be equal to zero, but then 
D„-i{S) ^ in that case by Corollary 111.41 1. 
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The next two results cover the case of equal degrees degq = degp: 

Theorem 3.24. The function p.49p . where degq — degp, is an R-function of negative type with exactly 
r (< n) poles if and only if 

r;2,+i(p,g) >0, j = 0,l,2,...,r, (3.56) 
r/,(p,g)=0, i>2r + 2, (3.57) 

where r]i(p,q) are the leading principal minors of the infinite Hurwitz matrix H{p,q) defined by (jl.l24p . 

Theorem 3.25. The function p.49p . where degq = degp, is an R-function of negative type with exactly 
r {< n) poles if and only if 

A2,(p,g) > 0, j = l,2,...,r, (3.58) 
A,(p,q) = 0, i = 2r + 2,2r + 3,...,2n+l, (3.59) 

where Ai{p,q) are the leading principal minors of the Hurwitz matrix 'H2n+i{p,q) defined in (I1.126p . 

All these results can be easily obtained from Theorenis ll.2l and l3.4l and from the formulae (|1.129p . (|1.131|) . 
(HIISSl), (dSSl). Then Theorem inn and the formula ^JW^), ^A^, ([LTTTI) . ([11^ imply 

Theorem 3.26. If the function p.49p . where degq < degp, is an R-function of negative type with exactly 
r {< n) poles, then the number of its positive poles, r+(< r), is 

r+ = V{f]i{p,q),'q3{p,q),r]5{p,q),.. . ,'i]2k+i{p, q)) = V(l, A2(p, g), A4(p, g), . . . , A2fc(p,q)), 

where k = r, if |-R(0)| < oo, and k = r — 1, if R{0) = oo. 

Theorem 3.27. If the function (j3.49l) . where degq = degp, is an R-function of negative type with exactly 
r {< n) poles, then the number of its positive poles, r+(< r), is 

r+ = Y{'q2{p,q),m{p,q),V6{P,q),-- ■,V2k+2{p,q)) V(Ai(p,(j), A3(p,g), A5(p,g),. . . , A2fc+i(p, g)), 
where k = r, if |-R(0)| < oo, and k — r — I, if R{0) — oo. 

In these theorems, the numbers of sign changes in the sequences of Hurwitz minors must be calculated 
by the Frobenius Rule because of the equalities (|1.129|) - (|1.132|) and (jl.l36l) - (|1.139p between Hurwitz 
and Hankel minors. 

We next address separately two extreme cases of Theorems 13 . 26f[3 .271 when all poles of the function R 
are either negative or positive: 

Corollary 3.28. The function p.49p is an R-function of negative type with exactly r (< n) poles, all of 
which are negative, if and only if 

77,(p,g)>0, j = l,2,...,fc, 
77,(p,g)=0, i>k, 

where k — 2r -\- 1 if degq < degp, and fc = 2r + 2 if degq = degp, but r^i{p,q) are the leading principal 
minors of the infinite Hurivitz matrix H{p,q) defined in (|1.123l) . 

Corollary 3.29. The function (|09)) is an R-function of negative type with exactly r (< n) poles, all of 
which are negative, if and only if 

A,(p,g)>0, i = l,2,...,fc, 
Ai(p,g)=0, i = fc + l,fc + 2,..., 

where k = 2r, if degq < degp, and k — 2r + 1, if degq = degp, but Ai{p,q) are the leading principal 
minors of the Hurwitz matrix 'H2n{p,q) or the matrix 'H2n+i{p,q) defined in (jl.l25p - (|1.126l) . 

Corollary 3.30. The function p.49p is an R-function of negative type with exactly r (< n) poles, all of 
which are positive, if and only if, 
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• for degq < degp, the inequalities p.52p " (|3.53p hold and 

V2]-i{p,q)mj+i{p,q) < o, j = i,2,...,r, 

where T]i{p,q) are the leading principal minors of the matrix H{p,q) defined by p.l23p . 

• for degq = degp, the inequalities p.56p - (|3.57p hold and 

V2j{p, q)V2j+2{p, q) <0, j = 1, 2, . . . , r, 

where r]i{p,q) are the leading principal minors of the matrix H{p,q) defined by (jl.l24l) . 

Corollary 3.31. The function p.49p is an R-function of negative type with exactly r (< n) poles, all of 
which are positive, if and only if, 

• for degq < degp, the inequalities p.54p - (j3.55p hold and 

A2j-2{p,q)^2j{p,q) <0, j = l,2,...,r, (Ao(p,g):=l) 
where Ai{p,q) are the leading principal minors of the matrix 'H2n{p,q) defined in (jl.l25p . 

• for degq — degp, the inequalities p.58p - (|3.59p hold and 

A2j-i{p,q)A2j + iip,q) < 0, i - l,2,...,r, 
where Ai{p,q) are the leading principal minors of the matrix 'H2n+i(p, 9)7 defined in (I1.126P . 

At last, there is one more way to find the number of negative (or positive) poles of an i?-function. 
This method turns into the famous Lienard-Chipart theorem when applied to the theory of Hurwitz stable 
polynomials (see [36l Chapter XV, Section 13] Theorem 11). But first, let us introduce (without proof) the 
remarkable but little-known consequence of the famous Descartes Rule of Signs (see, for example, 79, ) . 

Recall that we denote by V~(ao, ai, . . . , a„) the number of weak sign changes in a real sequence 
(ao, fli, . . . , a„), i.e., the number of sign changes with zero elements of the sequence removed. 

Theorem 3.32 ( [99[ 179)). // a real polynomial aoz" + oiz"^^ + • • ■ + a„ has only real roots, then the 
number of its positive zeros, counting multiplicities, is equal to V^(ao, ai, . . . , a„). 

This fact, together with previous results, implies the following theorem. 

Theorem 3.33. If the real rational function p.49p is an R-function of negative type with exactly n pole^^, 
then the number of its positive poles is equal to SC^ {ao,ai, . . . ,a„). In particular, R has only negative 
poles if and only i/E3 o-j > for j = 1,2, ... ,n, and R has only positive poles if and only if Oj^iOj < 
for j = l,2,...,n. 

Proof. Indeed, by Theorem l3.4l if the function R defined in p.49p is an i?-function, then the polynomial p 
has only real roots which are poles of the function R. The number of positive zeros of p (poles of R) equals 
V^(ao, ai, . . . , a„), according to Theorem 13.321 □ 

For two extreme classes of i?-functions, i.e., those with only positive or only negative poles, we can 
obtain a few more criteria. 

Theorem 3.34. The function p.49p . where degq < degp — n, is an R-function of negative type and has 
exactly n negative poles if and only if one of the following conditions holds 

1) On > 0, a„_i > 0, . . . , ao > 0, Ai{p,q) > 0, A3{p,q) > 0, . . . , A2„-i(p,g) > 0; 

2) a„ > 0, 6„ > 0, bn-i > 0, bi > 0, Ai{p,q) > 0, A3{p,q) > 0, A2n-i{p,q) > 0; 

3) a„ > 0, a„_i > 0, . . . , ao > 0, A2(p, q) > 0, Aiip, q) > 0, A2„(p, q) > 0; 

4) an >0, bn> 0, bn-i > 0, 61 > 0, A2{p,q) > 0, A4(p,g) > 0, A2„(p,g) > 0. 

•^'^The latter means that gcd(p, g) = 1. 

^^In fact, the coefRcients must bo simply of the same sign, but we already assumed that ao > (see | |3.50| |V 
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Proof. Condition 1) contains the inequalities p.54p (with r = n). By Theorem 13. 231 this is equivalent to 
the function R being an i?-function of negative type with exactly n poles. Now Theorem 13 . 331 implies that 
the positivity of the coefficients of p is equivalent to the negativity of poles of R. 

Condition 2) also contains the inequalities (I3.54[) (with r — n), which are equivalent to the function R 
being an i?-function of negative type with exactly n poles. By Theorem 13. 4[ all zeros and poles of R are 
real and simple. 

If the function R has exactly n negative poles, then from Theorem l3.4l it follows that it also has negative 
zeros (since the zeros and the poles of R are interlacing) . Now Theorem 13.321 yields the positivity of all 
coefficients of the polynomials p and q. 

Conversely, if the inequalities a„ > 0, 6„ > 0, 6„_i > 0, . . ., 5i > hold, then R has only negative 
zeros, according to Theorem 13.321 The interlacing of zeros and poles of R implies that R may have at 
most one nonnegative simple pole, that is, the polynomial p may have at most one nonnegative simple 
zero. But aq > by assumption, therefore, p{z) +oo as z — +oo. Since all zeros of p are simple and 
real, we have a„ = p(0) < whenever p has one nonnegative zero. This contradicts the inequality a„ > 0. 
Consequently, p has only negative roots, and R has only negative poles. 

If the function R is an i?-function of negative type with exactly n negative poles, then Corollarv 13.291 
and Theorem 13.321 implv condition 3). 

Now suppose that condition 3) holds. It contains the inequalities A2j{p, q) > 0, j — 1,2, . . . ,n, which 
are equivalent to the inequalities {—l)^Dj{R) > 0, j = 1,2, . . . ,n, according to (|1.137p . But we know that 
{—lyDj(R) — Dj{F), where F{z) = —zR{z), and the positivity of these determinants is equivalent to the 
function F being an i?-function of negative type by Theorem 13.41 From the same theorem it also follows 
that all zeros of the polynomials zq{z) and p(z) are real simple and interlacing. From Theorem 13.321 and 
from the inequalities aj > 0, j = 0, 1, . . . , n, we obtain that p has only negative roots. This fact implies 
that all zeros of q are negative too and that they interlace the zeros of the polynomial p since otherwise 
the zeros of zq{z) would not interlace the zeros of p{z). Since the zeros of p and q interlace, the function 
R is an i?- function of negative type by Theorem 13.41 It has only negative poles (the zeros of p). 

As in the case of condition 3), condition 4) holds if the function R is an i?-function of negative type 
with exactly n negative poles. 

If condition 4) holds, then, as before, the function F{z) = —zR{z) is an iZ-function of negative type 
and, therefore, the roots of the polynomials zq{z) and p{z) are real and simple and interlace each other. 
The inequalities hj > Q, j — 1,2, . . . ,n, imply that the polynomial q has only negative zeros. Consequently, 
p{z) has at most one simple nonnegative zero because of the interlacing with the zeros of zq{z). As above, 
p cannot have a nonnegative zero since otherwise the inequality a„ > cannot hold. Since the polynomial 
zq{z) has nonpositive zeros but p{z) has only negative zeros and since their zeros interlace each other, the 
zeros of q and p interlace too. Now Theorem 13.41 implies that R is an i?-function of negative type with 
exactly n negative poles. □ 

In the case degg = degp we obtain similar criteria. The next theorem addresses the situation when 
the degrees degp and degg are equal and all poles of R are negative. 

Theorem 3.35. The function (13.491) . where degq = degp = n, is an R-function of negative type with 
exactly n negative poles if and only if one of the following conditions holds 

1) a„ > 0, a„_i > 0, ao > 0, A2{p,q) > 0, A4{p,q) > 0, A2„(p,g) > 0; 

2) a„ > 0, h„ > 0, bn-i > 0, 6o > 0, A2{p,q) > 0, A4{p,q) > 0, A2n{p,q) > 0; 

3) a„ > 0, a„_i > 0, . . . , oq > 0, Ai{p,q) > 0, A3{p,q) > 0, A2„+i(p,g) > 0; 

4) a„ > 0, 6„ > 0, bn-i > 0, &o > 0, Ai{p,q) > 0, A3{p,q) > 0, A2„+i(p,(j) > 0. 

Proof. This theorem can be proved in the same way as Theorem 13.341 using Theorem 13.251 instead of 
Theorem 13.231 and the equalities ()1.139|) instead of the equalities ()1.137|) . □ 

To address the situation when all poles of R are positive, it is enough to switch to the function 
z i~> ~R{~z) and apply the same methods as in the proof of Theorems 13.341 and 13.351 This yields the 
following two results, one for the case degg < degp and the other for the case degg — degp. 
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Theorem 3.36. The function p.49p . where degq < degp = n, is an R-function of negative type and has 
exactly n positive poles if and only if one of the following conditions holds 

1) (-l)"a„ > 0, (-l)"-ia„_i > 0, . . . , ao > 0, 
Ai {p, q) > 0, A3 {p,q)>0, A2„_i (p, q) > 0; 

2) (-!)"«„ > 0, > 0, (-l)"~26„_i > 0, . . . , &i > 0, 
Ai (p, q) > 0, A3 (p, <z) > 0, . . . , A2„-i (p, q) > 0; 

3) (-l)"a„ > 0, (-l)"-ia„_i > 0, . . . , ao > 0, 
-A2(p,(z) >0, A4(p,<z) >0, (-l)"A2„(p,g) >0; 

4) (-l)"a„ > 0, > 0, > 0, . . . , &i > 0, 

-A2(p,(z)>0, A4(p,g)>0, (-l)"A2„(p,<7) >0. 

Theorem 3.37. The function p.49p . where degq — degp, is an R-function of negative type with exactly 
n positive poles if and only if one of the following conditions holds 

1) (-l)"a„ > 0, (-l)"-ia„_i > 0, . . . , ao > 0, 
A2ip,q) > 0, Aiip,q) > 0, A2n{p,q) > 0; 

2) (-l)"an > 0, > 0, (-l)"-25„_^ > 0, . . . , 61 > 0, 
A2{p,q) > 0, Ai{p,q) > 0, A2„(p,(7) > 0. 

3) (-l)"a„ > 0, (-l)"-ia„_i > 0, . . . , ao > 0, 
Ai(p,g)>0, -A3(p,g)>0, (-l)"A2„+i(p, g) > 0; 

4) i-lTon > 0, (-l)"-i6„ > 0, > 0, . . . , 61 > 0, 
Ai(p,<7) >0, -A3(p,g) >0, (-l)"A2„+i(p,<z) >0; 

3.4 i?-functions with a Stieltjes continued fraction expansion 

We now provide criteria for R functions based on the coefficients of their Stieltjes continued fractions. 
Assume that our function (I3.49|) with exactly r (< n) poles has a Stieltjes continued fraction expansion 

R{z) = co + , Cj-eR, Cj^O, (3.60) 

ciz H 



C2 



1 

C3Z H 



1 



where 



T.^- if |i?(0)|<oo, ^3^^^ 

I C2r-lZ if -R(O) = 00, 



and Co = s_i, co 7^ if and only if degp = degq(= n). From Theorem 11.461 it follows that the inequali- 
ties p.l06p - (ll.l07p hold for the function R. At the same time the coefficients Ci in (I3.60|) can be found by 
the formulae (I1.114p - (ll.ll5p (see (|1.116p V 

Theorems EH and \2J5\ yield 



60 



Theorem 3.38. Let the function (|3.49p with exactly r {< n = degp) poles have a Stieltjes continued 
fraction expansion (j3.60l) - (|3.6ip . The function R is an R-function of negative type if and only if 

C2j-i>0, i = l,2,...,r. (3.62) 

Then the number of negative poles is equal to the number of positive coefficients C2j , j — 1,2, . . . ,k, where 
k ^ r, if \R{0)\ < oo, and k = r - I, if R{0) = oo. 

Note that i?-functions with only positive or only negative poles have Stieltjes continued fraction expan- 
sions, according to Corollaries l3.10l and l3.12l and Theorem ll.461 Theorem I3.38l has the following corollaries. 

Corollary 3.39. A rational function R with exactly r poles is an R-function of negative type with all 
positive poles if and only if R has a Stieltjes continued fraction expansion (I3.60|) - p.61l) , where inequali- 
ties (I3.62p hold and where 

C2j < 0, j = l,2,...,r. 

Corollary 3.40. A rational function R with exactly r poles is an R-function of negative type with all 
positive poles, except for one at 0, if and only if R has a Stieltjes continued fraction expansion p.60p - 
p.6ip . where inequalities (|3.62l) hold and where 

C2j<0, i = l,2,...,r-l. 

Corollary 3.41 (Markov, Stieltjes [HZl [231 HH [SSI US])- A rational function R with exactly r poles is 
an R-function of negative type with all negative poles if and only if R has a Stieltjes continued fraction 
expansion p.60p - (|3.6ip . where 

Q>0, i = l,2,...,2r. (3.63) 

Corollary 3.42. A rational function R with exactly r poles is an R-function of negative type with all 
negative poles, except for one at 0, if and only if R has a Stieltjes continued fraction expansion p.60p - 
p.6ip . where 

c^>0, i = l,2,...,2r~l. (3.64) 
In the sequel we use the following well-known result. 
Theorem 3.43 (Aisen, Edrei, Schoenberg, Whitney, [U [2j [25l [54]). The polynomial 

g{z) ^ gnz^ + giz^-^ + ■ ■ ■ + gi 

has only nonpositive zeros if and only if its Toeplitz matrix T{g) defined by (|1.128p is totally nonnegative. 

The functions whose Taylor coefHcients generate totally nonnegative Toeplitz matrices of the form T{g) 
was introduced by Schoenberg [85l [86] and also studied by Edrei [27l [26l [5] . 

We next prove a criterion of total nonnegativity of infinite Hurwitz matrices. Previously, only one 
direction, i.e., the fact that the infinite Hurwitz matrix of a quasi-stable polynomial^ is totally nonnegative 
(see [6l [57l [46l [55] ) , was known. The necessary and sufficient condition was known only for finite Hurwitz 
matrix of Hurwitz stable polynomials. 

Theorem 3.44 (Total Nonnegativity of the Hurwitz Matrix). The following are equivalent: 

1) The polynomials p and q defined by p.50p - (|3.5ip have only nonpositive zero^^. and the function 
R = q/p is either an R-function of negative type or identically zero. 

2) The infinite matrix of Hurwitz type H{p, q) defined by (|1.123p - (|1.124p is totally nonnegative. 

^^The quasi-stable polynomials are polynomials with zeros in the closed left half-plane of the complex plane. 
^^Here we include the case when q{z) = 0. 
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Proof. If condition 1) holds and R{z) ^ 0, then, according to Corollaries l3.41| -[ 3.42[ the function R has 
a Stieltjes fraction expansion (|3.60p - (|3.61l) and satisfies the inequalities p.63p or p.64p . where r (< n) 
is the number of poles of the function R — q/p. According to Theorem ll.46[ the matrix H{p,q) has a 
factorization (|1.140p or (|1.14ip . where all matrices J{cj) are totally nonnegative by inspection since all cj 
are positive. However, by assumption, all zeros of p and q are nonpositive, therefore all zeros of g = gcd(p, q) 
are also nonpositive. Now from Theorem l3.43l we obtain that the matrix T{g) is totally nonnegative. Since 
the matrix H{0, 1) is trivially totally nonnegative, the matrix H(j), q) is totally nonnegative as a product 
of totally nonnegative matrices (see [34] ) . 

Let condition 1) hold and R{z) = 0. Then instead of the factorizations (|1.140p - (|1.14ip . we obtain 
H{p, q) = H{0, l)T(p). So in this case H{p, q) is also totally nonnegative. 

Conversely, let the matrix H{p, q) be totally nonnegative and q{z) = 0. In this case, T{p) is also totally 
nonnegative as a submatrix of H(j), q). By Theorem 13. 43( p has only nonpositive zeros, as required. 

Let now the matrix H{p,q) be totally nonnegative and let q(z) ^ 0. All submatrices of H{p,q) are 
also totally nonnegative. Since matrices T(p) and T{q) are submatrices of H{p, q) constructed using all 
its columns and all even or all odd rows, they are totally nonnegative. Thus, according to Theorem 13.431 
the polynomials p and q have only nonpositive zeros. Let g be the greatest common divisor of p and q, 
so that p = pg and q = qg. Then Theorems 11.431 and 13.431 implv that H{p,q) = H{p,q)T{g), where the 
matrix T{g) is totally nonnegative since g has only nonpositive zeros. 

First, assume that degp < degg, that is, i?(oo) = 0, and use notation 



p{z) 



fo{z) 

/i(^) 



af'>z^' 



,(0) 



where a^^ = > by assumption (sec (13.5011 ). From the total nonncgativity of the matrix H(p, q), which 



is the same as H{fQ, /i), we have af'' > 



-,(1) 



> 0, i = 1,2, 



,(1 



Show that deg/i = rt — 1. To this end, let us introduce notation Hq :=iJ(/o,/i) and suppose that 

> for some integer 1 < j < n — 1. In this case, we have 



,(1) 



= am 



J + 2 



(0) 


(0) 







a\ 















(0) 





< 0. 



This inequality contradicts the total nonnegativity of the matrix Hq. Thus, we obtain that a'^'^ > and, 
therefore, the polynomial /i is of exact degree n — I. 

Now we can perform the first step of the algorithm (|1.119p - (ll.l20p to get the next polynomial /2: 



/2(^) = foiz) - ciz/i(z) 



(2) 



(3.65) 



JO) 



where ci = -j-r > 0. As in the proof of Theorem 11.461 we obtain the factorization 
% 

Ho = J{ci)Hi, (3.66) 
where Hi :=i?(/2,/i) and the matrix J(ci) is defined as in (|1.142p . If f2{z) = 0, then /i — gcd(/o,/i) 
and the function R{z) ~ is an i?-function of negative type. So, in this case, the theorem is proved. 

Cl z 

If f2{z) ^ 0, then from (I3.65P and from the total nonnegativity of the matrix H{fo, fi) it follows that 



,(2) _ 



1 



*0 "i+l 



1 



(0) (1) 



Hn 



3 

i + 3 



> 0, i = 0,l,. 



1. 



Suppose that deg/2 < n — 1. 



,(2) 



(2) 



,(2) 



= ij^i = and 



Then there exists an integer j, 1 < j < n 

(2) 

> 0. Then from (I3.66P we obtain 



1, such that the coefficients 



Hn 



4 

J+3 



= ci-Hi 



3+3 



= Cl a, 



,(1) 



^^This is possible since q{z) ^ by assumption. 



,(1) 



(2) 
.(1) 



< 0, 
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which contradicts the total nonnegativity of the matrix Hq. So, the coefficient Oq must be positive and, 
therefore, deg/2 = n — 1. Thus, we can perform the next step of the algorithm (|1.119l) - (|1.120p to obtain 
the next polynomial f^: 



(3) 



-2' 



(3.67) 



where C2 = > 0. The formulae ^M^-^Mi imply 



,(3) 



1 



(2) (2) 
"0 "i+l 



-,(2; 



Hi 



13 4 5 
12 3 i + A 



1 



(0) (1) (2) 

Oq Cq 



■Hn 



12 3 4 
12 3 i + 4 



> 0, 



4^) (a«)' 
i = 0,l,...,n-2. 



ci 



12 3 4 

12 3 i + A 



For the matrix -ff (/o, /i), we have the factorization 

Ho = J{CI)J{C2)H2, 

where i/2 = i?(/2,/i). 

Now let us assume that we have a sequence of polynomials /o, /i, . . . 
from the polynomials /o and /i by the algorithm ()1.119p - ()1.120p . that is, 

hi{z) = f2i-2{z) - C2i^izf2i-l{z), i = l,2,.. 

hi+l{z) = f2i-l{z) - C2if2i{z), i = l,2,.. 



fm (1 < m < 2r) constructed 



(3.68) 



where deg/2i-i = deg/2i = n ~ i. Let us also assume that all the coefficients ap"* of each polynomial / 
are nonnegative. Then the coefficients Cj in (13.681) are positive and determined by the formula 



,0-1) 



a = 



>o, 



j = l,2,...,m- 1. 



(3.69) 



Moreover, if we set 

Tj J H{Jm-i,fm) if m is odd, 

[ H{fm, fm-i) if m is even, 

then the matrix Ho has the following factorization 

Ho = J(ci)J(c2) • • • J{Crn-l)H,n-l, 



(3.70) 



as it follows from the proof of Theorem 11.461 

Let us perform the next step of the algorithm (I1.119l) - (|1.120p and obtain the next polynomial fm+i- 



fm+liz) = 



fm-l{z) - CmZfm{z) if 771 is odd, 
/m-l(z) - Cm.fm{z) if 771 is CVCU, 



(3.71) 



and denote the coefficients of fm+i by a'f^^^K If fm+iiz) = 0, then = 2r — 1 or tjt, = 2r — 2 (r is the 

/i 

number of poles of the function R ~ —) and /,„ = gcd(/o, /i). Now the formula p.69|) and Corollaries l3.41l 

and l3.42l show that R is an i?-function of negative type with nonnegative poles, which completes the proof 
in this case. If fm+i{z) ^ 0, then it follows from (I3.70p - p.7ip and from the total nonnegativity of the 
matrix Ho that 



(m+l) 



1 



a 



"o "o 



(m) 



Ho 



1 2 ... 7n+l 777 + 2 

1 2 ... 777+1 7 + 777 + 2 



> 0, i ==0,1,2,..., n- 



-1. 
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We will show that a^Q^~^^^ > 0. In fact, if we suppose that it is not true, then there exists a number j 

{l<j<n- - -1) such that a[,'"+'^ = a^™+'^ = • • • = a^.™+') = and a^"+'^ > 0. Then yields 
L 2 J 

„ fl 2 ... m + 2 m + 3\ 2 m f {m)Y'+ 

2 ... rn + 2 , + m + 2 ^ ^^^^ ' ' ' ^™ ("o j 



4"+^' 

(m) (m) 



< 0. 



(ni I 1) 

This inequality contradicts the total nonnegativity of the matrix Hq. Consequently, Oq > and we 
can run the next step of the algorithm (|1.119p - (|1.120p to obtain the next polynomial /m+2- 

Thus, step by step we construct a sequence of positive numbers ci, C2, . . . , Cfc, where fc = 2r if |i?(0)| < 
00, and k = 2r—\ otherwise. These numbers are exactly the coefficients of the Stieltjes continued fraction 
expansion of the function R. In this case, R is an i?-function of negative type with nonpositive poles 
according to Corollaries 13.411 and 13.421 as required. 

If degp = degg, then R{oo) — cq — — > 0, since H{p,q) is totally nonnegative. If we denote /_i = p 

and fo = q and run the algorithm (|1.119p - (|1.120l) as before, then we obtain that the function R = p/q is 
an i?-function of negative type with nonpositive poles. □ 

Remark 3.45. Note that the easier direction 1) =^ 2) of Theorem l3.44l was proved in [351 Proposition 3.22] 
as a generalization of results due to Asner, Kemperman and Holtz [51 [571135], but the more complicated 
implication 2) =^ 1) appears to be new. 



For finite matrices of Hurwitz type, there is no criterion analogous to Theorem 13.441 However, the 
following two theorems can be derived as straightforward consequences of Theorem 13.441 

Theorem 3.46. // the polynomials p and q defined in (j3.50p - (|3.51l) have only nonpositive zeros, and the 
function R = q/p is either an R- function of negative type or identically zero, then the finite matrix of 
Hurwitz type Hkip, q) is totally nonnegative. Here k = 2n ifdegq < degp, and fc = 2n+ 1 if degq — degp. 

Theorem 3.47 (Total Nonnegativity of the Finite Hurwitz Matrix). Given the polynomials p and q 
defined in (j3.50p - (|3.5ip . the function R ~ q/p is an R-function with exactly n negative poles if and only 
if the finite matrix of Hurwitz type Hk [p, q) is nonsingular and totally nonnegative. Here k ~ 2n if 
deg q < degp, and k = 2n + I if deg q = degp. 

Proof. Indeed, if the function R is an i?-function with exactly n negative poles poles, then the function 
zR{z) also has exactly n poles, so the corresponding Hankel minor Dn{R) is nonzero according to Theo- 
rem [TT3l By Theorem 11.451 this means that defHfc(p, g) 7^ 0, so the matrix Hk{p,q) is nonsingular. But 
HkiPi q) is totally nonnegative as a submatrix of the totally nonnegative matrix H{p, q) (see Theorem l3.44p . 

Conversely, let the matrix Hkip, l) be nonsingular and totally nonnegative. The nonsingularity of the 
matrix T-LkiPiQ) implies that the function R ~ q/p has exactly n — Aegp poles, that is, the polynomials p 
and q are coprime. Now by the same methods as those used in the proof of Theorem 13. 44| we can show 
that the total nonnegativity of the matrix "Hfc {p, q) implies that the function R has a Stieltjes continued 
fraction expansion p.60p - (l3.61l) with positive coefficients, which is equivalent to the function R = q/p 
being an i?-function with negative poles, according to Corollarv 13.411 □ 

In the particular case when p and q are the even and odd parts of some polynomial. Theorem 13.471 was 
first established by Asner [5]. 

Let the polynomials p and q be defined in (j3.50l) - (|3.5ip and let fc = 2n if deg q < degp, and fc = 2rt + 1 
if deg q = degp. In the same way as in Theorem II. 431 one can show the following: if the polynomials have 
a common divisor g of degree I such that p = pg and q = qg, then 

Hk{p,q)=Hk{p,q)Tk{9), (3.72) 

where the matrix T-Lk{p^q) is the fc x fc principal submatrix of the infinite matrix H{p, q) indexed by rows 
(and columns) 2 through fc + 1 , and the matrix Tk{g) is the fc x fc leading principal submatrix of the matrix 
T{g) defined in (111281) . 

If the matrix Hkip, q) is singular and totally nonnegative, then "Hkip, q) can be represented as in p.72p . 
where the polynomials p and q have only nonpositive zeros and R = q/p is either an _R-function or 
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R{z) = 0, but the polynomial g = gcd(p, q) has no nonpositive zeros or g{z) = const 0. This factorization 
of the totally nonnegative matrix Jikip, q) is possible, for example, if all minors of order < fc of the infinite 
matrix T{g) are nonnegative. 

Remark 3.48. If all minors of order < A: of the infinite matrix T{g) are nonnegative, then the sequence of 
the coefficients of the polynomial g is called k-times positive or k-positive. If T{g) is totally nonnegative, 
then the sequence of the coefficients of the polynomial g is called totally positive [871 HB] ■ The functions 
generating fc-positive (totally positive) sequences are usually denoted by PFk (PFoo)- 

Based on the results above, we make the following conjecture. 

Conjecture 3.49. Given two polynomials p and q defined in p.50p - (|3.51|) . the finite matri^^ Hkip, q) 'is 
totally nonnegative if and only if p and q have only nonpositive zeros and R = q/p is either an R-function 
or R{z) = 0, and the polynomial g — gcd(p, q) has no real zeros and belongs to the class PFk-degg- 

Theorems 13.441 and 13.461 implv the following corollaries. 

Corollary 3.50. The following conditions are equivalent: 

1) The polynomials p and q defined by (j3.50p - (|3.5ip are coprime and have only negative zeros, and the 
function R = q/p is an R-function of negative type. 

2) The infinite matrix of Hurwitz type H{p,q) defined by (|1.123p - (|1.124p is totally nonnegative and 
Vkip,q) > 0, where k ^ 2n ifdegq < degp, and k = 2n + 1 ifdegq = degp. 

Corollary 3.51. The following conditions are equivalent: 

1) The polynomials p and q defined by p.50p - p.51l) are coprime and have only nonpositive roots, and 
the function R — q/p is an R-function of negative type. 

2) The finite matrix of Hurwitz type Hk{p,q) is totally nonnegative of rank k — 1, where k = 2n if 
deg q < degp, and fc = 2n + 1 if deg q — degp. 

Sometimes it is convenient to use the inverse indexing of polynomial coefficients. We now state and 
prove a result analogous to Theorem 13. 44( using this alternative ordering of coefficients. 

Corollary 3.52. For the polynomials 



g{z) = ao + aiz + 022^ + • • • + a„_iz" ^+a„z", a„ > 0, oq 7^ 0, 
h{z) = bi+b2Z + --- + 6„_iz"-2 + 6„z"-i, 6„ > 0, 



(3.73) 
(3.74) 



the following conditions are equivalent: 

1) The polynomials g and h have only negative zeros, and the function R — h/g is an R-function of 
negative type. 

2) The infinite matrix 



Hoo{g,h) := 



is totally nonnegative. 

Proof. Let the polynomials g and h have only negative zeros such that R — h/g is a,n i?- function of 
negative type with exactly m (< n) poles. Therefore, by Theorem 13.41 R can be represented as follows 





ai 


0-2 


0-3 


04 


05 


...\ 





bi 


b2 


b3 


^4 









ao 


ai 


02 


0-3 


04 










bi 


b2 


b3 


b4 




V : 












■■) 



R(z) 



bl + 62Z + ■ ■ ■ + bn-iz'^-^ + b„z"-l 

ao + aiz + a2z'^ + • • • + a„_iz""i + a„z" 



E 



z + Ao 



aj,\j > 0. 



•^^Here fe = 2n if deg q < degp, and fc = 2n + 1 if deg q = degp. 
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All zeros of this function are also negative. Therefore, bi 7^ 0. Consider the function 



R{z) = -R (-) = , 



aiz'' 



a2Z 



n-2 



1/A, 



aj,Aj>0. (3.75) 



Thus, if i? is an i?-function of negative type with negative poles, then i? is also an i?-function of negative 
type with negative poles. It is easy to show that the converse statement is also valid. Now by Theorem l3.44l 
and by p.75p . we obtain the equivalence of the conditions 1) and 2) of the theorem. □ 

In the same way, one can prove the following corollary. 
Corollary 3.53. For the polynomials 



g{z) = ao + aiz + a2Z^ + 
h{z) = bo + biz + b2Z^ + ■ 

the following conditions are equivalent: 



_n— 1 1 n 
-\Z + a„z , 



a„ > 0, 

6„>0, 60^0, 



(3.76) 
(3.77) 



1) The polynomials g and h have only negative zeros, and the Junction R = g/h is an R-function of 
negative type. 



2) The infinite matrix 



Hooig, h) 



/&o bi 62 &3 bi 65 

ao ai 02 03 04 

60 bi 62 &3 bi 

flo fli 02 03 

v; ; ; ; ; ; 



is totally nonnegative. 



Corollaries I3.52H3.5^ imply the following two theorems, which, in fact, define a class of interlacing 
preservers (see |28l and references there) and imply one theorem originally proved by Polya. 



Theorem 3.54. Let the polynomials g and h defined in (j3.73p - (|3.74p have only negative zeros, and let 
the function R = h/g be an R-function of negative type. Given any two positive integers r and I such that 
rl < n < (l + l)r, the polynomials 



hr,l{z) 



flo + OrZ + a2rZ^ + O.'irZ^ + ■ 



+ ariz , 



+ briz'-^ 



br + b2rZ + fesrZ + • • • + Orli 

have only negative zeros, and the function Rr_i = hr^i/gr,i is an R-function of negative type. 



(3.78) 
(3.79) 



Proof. Indeed, if the polynomials g and h have only negative zeros, and if the function R = h/g is an 
_R-function of negative type, then by Corollarv l3.53[ the matrix Hoo{g, h) is totally nonnegative. Then all 
its submatrices are totally nonnegative. In particular, the following infinite submatrix whose columns are 
indexed by 1, r + 1, 2r + 1, 3r + 1, . . . and rows are indexed by 1, 2r + 2, 4r + 3, 6r + 4, . . . 



H 00 {gr, I, hrj 



br 

ao 




^2r ^3r ^br 

b2r b^r &4r &5r 

Or f^2r ^3r 047. 

br b2r &3r ^4r 



is totally nonnegative. Now by CoroUarv 13.531 the polynomials p.78p - p.79p have only negative zeros, 
and Rr^i = hrj/gr,i is an i?-function of negative type, as required. □ 



The following theorem can be proved in the same fashion. 
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Theorem 3.55. Let the polynomials g and h defined in (j3.76l) - p.77p have only negative zeros, and let 
the function R — g/h be an R-function of negative type. Given any two positive integers r and I such that 
rl < n < {I + l)r, the polynomials 

gr,i{z) := ao + arZ + a2rZ^ + azrZ^ + ■ ■ ■ + ariz\ 

hr,l{z) ba + brZ + b2rZ^ +b^rZ^ + ■ ■ ■ +WlZ^ 

have only negative zeros, and the function Rr,i gr.i/hr,i is an R-function of negative type. 

Theorems 13. 54I45.55I imply the foUowing result of Polya [78, p. 319] (also implied by Theorem 13.431) . 
Corollary 3.56. // the polynomial 

g[z) = flQ + aiz + 022:^ + ■ ■ ■ + a„_iz"^"'" + a„z", a„ > 0, (3.80) 

has only negative zeros, then for any positive integers r and I satisfying rl < n < {I + l)r, the polynomial 

9r,i{z) = ao + a,.z + a2rZ^ + ■ ■ ■ + a^z' 

also has only negative zeros. 

An analogous result can be obtained for polynomials with only positive zeros. 

Corollary 3.57. Let the polynomial (|3.80p have only positive zeros. Then for any positive integers r and 
I satisfying rl < n < {I + l)r, the polynomial (|3.56p also has only positive (negative) zeros whenever r is 
odd (even). 

4 The number of distinct real zeros of polynomials. Polynomials 
with all real zeros 

In this section we present a sample application of the theory developed in the previous sections. Using 
those methods, we analyze the distribution of zeros of real polynomials with respect to the real and the 
imaginary axes. 

We should note that polynomials with all real roots have been studied in control theory, where this 
property is referred to as aperiodicity. Among the many relevant papers we note the work of Jury, Meerov, 
Fuller and Datta [331 1311 1311 1301 EEl HIl 113] containing special cases of our results, albeit derived using 
mostly different methods. 



4.1 The number of distinct positive, negative and non-real zeros of polyno- 
mials. Stieltjes continued fractions of the logarithmic derivative 

Consider a real polynomial 

p(z) = floz" + fliz""-^ H |-a„, oi, . . . , a„ e R, ao > 0, n > 1. (4.1) 

We denote its logarithmic derivative by L{z): 

^^^^ d\og{p{z)) _ p'{z) _ napz'^'^ + (n - 1)qiz"~^ H h On-i 

dz p{z) aoz" + aiz"-i + • • • + a„ 

If 

p{z) = ao(z - Ai)"^(z - A2)"^ . . . (z - A„i)"", ni + n2 H h n,„ = n, 

where uj is the multiplicity of the zero Xj (j — 1,2, ... , m), then the logarithmic derivative of the polyno- 
mial p has the following form 

m 

j=i J 
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Moreover, if we expand the function L into its Laurent series at oo 

L{z) = — + — + — + — + ■■■ , (4.3) 

z z-^ Z'^ 

then the coefficients Sj are the Newton sums of the polynomial p (see, for instance, |36)): 

m 

Sk = Y,n,\], fc = 0,1,2,.... (4.4) 

It is easy to see from (|4.2I) that the number of poles of the function L equals the number of distinct zeros of 
the polynomial p. But the Cauchy index Indl^(L) equals the number of distinct real zeros of p, since all 
poles of L are simple and since the residue of L at each pole is positive, [STJ[35]. Also from (|4.2p - (l4.3p and 
from Theorem ll.2l it follows that the rank of the matrix S = ||si+j||[f consisting of the Newton sums (j4.4p 
is finite and is equal to m (< n), the number of distinct zeros of p. 

As before, we denote by Dj{L) {j — 1,2,..., m) the leading principal minors of the matrix S{L) — 
||sj+j||g° (see (mH)). Then from Theorem [2TTT] we obtain 

Indt^(L) = m - 2 V(l, D,{L),D2{L), . . . , D„^{L)). (4.5) 

The above facts and the formula (|4.5p imply the following theorem. 

Theorem 4.1 f |361 161)). The number of distinct pairs of non-real zeros of the polynomial p equals 
Y{l,Di{L),D2{L),...,D^{L)). 

Corollary 4.2. 

V(l,Z?i(i),Z?2(i),...,i?™(i)) < [j\ ■ 

We also consider the determinants Dj(L) {j ~ 1,2, ... ,m) defined by (|1.9p and introduce the following 
notation for counting zeros: 

Definition 4.3. Let r denote the number of distinct real zeros of the polynomial p and let r+ and r^ be 
the numbers of distinct positive and negative zeros of p, respectively. 

Theorem 12.141 implies the following simple fact. 

Theorem 4.4. Let the numbers k and I be defined as follow^j\: 

k =. V(l, Di{L),D2{L), Drn{L)), I = V(l, Di{L),D2{L), D^{L)), 

Then the number of distinct pairs of non-real zeros of the polynomial p equals k and 



(4.6) 

P(0) ^ 0, 
if p(0) = 0. 

Proof. By Theorem 14. 4[ k is the number of distinct pairs of non-real zeros of the polynomial p. Since 
Indl|^(i) = r as mentioned above, we conclude 

r = m-2fc. (4.7) 

Now we observe that 

r+=Ind+°°(i), r- ^lnd'_^{L). (4.8) 
If p(0) ^ 0, then r = r+ + and |L(0)| > 0. So, from (gl]) and ^^-^^ we obtain 

r^ = I — k, r^ = m — k — I. 



r = 


m — 2k; 


r^ = 


l-k; 




m — k — I 


r+ = < 






m — fc — / — 1 



2'' Recall that the numbers V{1, Di{L),D2{L), . . .,Dm{L)) and V(l, Di{L),D2{L), . . .,Dm{L)) must be calculated accord- 
ing to Frobcnius Rule 12.41 
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Now let p{0) = 0. Then 

r = r+ + + 1 = m - 2k, (4.9) 

In this case, the number cti equals 1 in the formulae (I2.39p - (|2.40p since all the residues of the function L 
are positive (see (|4.2p '). whereas the number fT2 equals (Ti since all poles of L are simple. Thus, from (|4.8p 
and (|2.39p - (|2.40l) we obtain = I — k, r"*" = m — A; — ^ — 1, as required. □ 

Corollary 4.5. 

V(l, Di{L),D2{L), . . . , D„,{L)) < V(l, i5i(L), D2{L), . . . , 5„(L)). 

Remark 4.6. Since sq = n, we have Di{L) — sq — n > Q. This fact means that the polynomial p has at 
least one zero. If Dj{L) = for j > 2, then p has exactly one zero of multiplicity n. 

Our next statement is a slight modification (we use another continued fraction) and generalization (we 
cover the case p{0) — 0) of Theorem 3.5 from [65] (see also [84]). However, [6^ uses different methods. 



Theorem 4.7. Let the polynomial p be defined by (|4.ip . Then its logarithmic derivative L has a Stieltjes 
continued fraction expansion 

L{z) = ^ e M, Cj^O, (4.10) 

ciz H J 

C2 + ^ 



C3Z 



1 
T 



T = I ^ °' (4.11) 



^iz if p{0) = 



zf and only if L satisfies the inequalities 



D,{L)^0, j = l,2,...,m, 

D,{L)^0, j = l,2,...,m-l. 

Dj{L) = Dj{L) ^0, j > m. 

where m < degp. 

In that case, p has exactly m distinct zeros. Moreover, if the number of negative coefficients C2j_i 
equals k, and the number of positive coefficients C2j equals I, then the number of distinct pairs of nonreal 
zeros of p and the number of its distinct real, positive and negative zeros are given by the formulas (j4.6p . 

Proof. The theorem follows immediately from Theorems ll.46[ 12.151 and 14.41 □ 

From (|1.114|) - (|1.115|) it follows that the coefficients q in (|4.10|) - (|4.1ip can be found as follows: 

where in = m if p(0) ^ 0, m = m — 1 if p(0) — 0, and Do{L) = 1, Dq{L) = 1. 

Thus, Theorem 14.41 expresses the numbers of positive, negative and non-real zeros in terms of the 
number of sign changes in the sequences of the minors Dj{L) and Dj{L), and Theorem 14.71 does the same 
in terms of the Stieltjes continued fraction of L, provided that L has such a continued fraction expansion. 
Now we will obtain formula for those numbers in terms of the coefficients aj of the given polynomial p. 
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Consider the following 2n x 2n matrix. 





(n — l)ai 


{n 


- 2)a2 . 


Oti-I 








0\ 


ao 


ai 








a„ 











nao 


(n 


- l)ai . 


2a„^2 













ao 




ai 


a„_2 


















nao 


3a„_3 


2an-2 
















ao 


an-3 


an-2 



















nao 


(n — l)ai 


a-n-i 





V 










ao 




a-n-i 


On/ 



(4.14) 



which is a Hurwitz-type matrix constructed with the coefhcients of the polynomials p and with the coeffi- 
cients of its derivative p'. According to Definition II. 421 the matrix 'D2n{p) is 'H2n{p, q) (see (jl.l25|) '). where 
q = p' . Denote the leading principal minors if the matrix 'D2n{p) by Sj{p), j = 1, 2, . . . , 2n. We remind the 
reader that 62n-i{p) = aoD(p), where D(p) is the discriminant of p (see (|1.29p ). 



Theorem 4.8. The polynomial p has exactly m < n distinct zeros if and only if 

^2m-i(p)7^0, and 5j{p)^0 for j > 2m. 

At the same time, if 

k = V{l,Si{p),S3{p),.. . ,(52m_l(p)), / = V{l,d2{p),S4,{p), ■ ■ .,S2m{p)), 

then the number of distinct pairs of non-real roots of the polynomial p equals k and 

r = m — 2k ; 



(4.15) 




1 



P(0)^0, 
if p(0)-0, 



where r is the number of distinct real roots of p, and r and r+ are the numbers of distinct negative anc 
positive roots of p, respectively. 



1,2,. 



,n; 
,n. 



(4.16) 
(4.17) 



Proof. By Definition [L42l we have 

6,{p)^A,{p,p'), J 
Therefore, from Theorem [OS] (see (I1.136|) - (ll.l37p ') we obtain 

S2j-iip)^al'-'D,{L), .7 = 1,2,.. 

S2,{p) = {-iyal'D,{L), j = l,2,.. 

Then, from Theorem 11.21 and ()4.16p it follows that the rank of the matrix S = llsj+fcHo" equals m if 
and only if the condition (I4.15P holds. At the same time, we have D,n-i{L) ^ 0, D,n{L) = if and only 
if p(0) = 0, and i5™(i) 7^ if and only if p(0) ^ 0, and Dj{L) = for j > m. 

It suffices to note that (|4.16|) implies that 

k = V(l, h{p)Mp). ■ • ■ , hm-M) = V(l, D,{L),D2iL), D^(L)) 
since ao > 0. But from ()4.17p we have 

I = V(l, 52{p).5i{p), . . ■,S2rn{p)) = V(l, -5i(L), D2{L), . . . , (-l)™A„(i)) 

= m- V(l, Di{L),D2{L), . . . , An(i)), 



(4.18) 



(4.19) 



which gives Y{l,Di{L),D2iL), 
Y{l,DiiL),D2{L),...,DrniL)) = 
and from Theorem 14.41 

^®In this case, Dm—i{L) 0, according to Corollarv ll.4l 



..,Drn{L)) = m - Mf DmiL) ^ 0, and i@ Dm{L) = 0, then 
m — 1 — L Now the assertion of the theorem follows from (j4.18p . (|4.19p 

□ 
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Remark 4.9. For the polynomial p we have 6i{p) = ao-Di(-L) = a^n > 0. 

Using the formula: (|4.16p - (|4.17l) . we can represent the cocfhcients Ci of the Stieltjes continued fraction 
expansion (|4.10p - (|4.1ip of the function L in terms of the determinants 6i{p): 



Si-2{p) -Siip)' 



i = 1,2, ...,2m, (4.20) 



where = — , So{p) = 1. 

ao 



4.2 Polynomials with real zeros 

We now provide several explicit criteria for a polynomial to have only real zeros. These criteria, just as 
those developed before, use the Hankel and Hurwitz minors made of the coefficients of a given polynomial. 
Let us again consider the polynomial 

p(z) = aoz" + aiz""^ H h a„, ai, . . . , a„ G K, oq > 0. (4-21) 

and let 



= + § + ^ + (4.22) 



p'jz) _ £o _^ £i _^ £2 
p{z) z z^ 

be its logarithmic derivative. Theorem 14.41 directlv implies the following results. 

Theorem 4.10. The polynomial p has exactly m (< n) distinct zeros, all of which are real, if and only if 

Dj{L)>0, j = l,2,...,m, 
D,iL)^0, j>m. 

Theorem 4.11. Let the polynomial p have exactly m (< n) distinct zeros all of which are real. Then the 
number of distinct negative zeros of the polynomial p equals V(l, Di{L), D2{L), . . . , Dm{L)). Moreover, 

D,{L)=0, j>m. 

and 

(Dm{L)^0 if p(0)^0, 

\ D^{L) = and D^_^{L) ^0 if p{0) = 0. 

Proof. The polynomial p takes value at if and only if the function L has a pole at zero. Therefore, 
the assertion of the theorem follows from CoroUarv 11.41 and Theorem 14.41 □ 



Remark 4.12. According to Theorem 13.321 if the polynomial p has only real zeros, then the number of its 
positive zeros equals V~ (oq, ai , 02, . . . , a„), i.e., the number of strong sign changes in the sequence of the co- 
efficients of the polynomial p. Obviously, the number of its negative zeros is equal to P~(ao, 01,02, ... ,a- 
i.e., the number of strong sign retentions in the sequence of its coefficients. 

Also from Theorem 14.71 we obtain the following simple corollary. 



Corollary 4.13. Let the polynomial p have exactly m (< n) distinct zeros, all of which are real, and let 
its logarithmic derivative L have the Stieltjes continued fraction expansion (|4.10p - (|4.1ip . Then 

C2j-i>0, j = l,2,...,m, (4.23) 

and the number of positive coefficients C2j equals the number of negative distinct zeros of the polynomial p. 
The coefficients Ci can be found by the formults (|4.12p - (|4.13p or (j4.20l) . 

Let us consider the Hankel matrix 

SiL) = ||,s,+,||o-, (4.24) 

made of the coefficients of the series (|4.22p . 

From Theorems [XTSl KM lilTl Coroharv lXTOl and Remark BT^ we obtain the following straight- 
forward consequences, the first two addressing the case when all zeros of p are positive, and the next two 
when all zeros are negative. 
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Corollary 4.14. The polynomial p has only positive zeros and exactly m (< n) of them are distinct if and 
only if the matrix S{L) defined by (|4.24p is strictly totally positive of rank m. 

Corollary 4.15. The polynomial p has only positive zeros and exactly m (< n) of them are distinct if and 
only if the matrix S{L) is positive definite of rank m and aj-iaj < (j — 1,2, . . . ,n). 

Corollary 4.16. The polynomial p has only negative zeros and exactly m (< n) of them are distinct if and 
only if the matrix S{L) is strictly sign regular of rank m. 

Corollary 4.17. Let all the coefficients of the polynomial p be of the same sign. Then all zeros of the 
polynomial p are negative and exactly m (< n) of them are distinct if and only if the matrix S{L) is positive 
definite of rank m. 

li m — n in Theorems l4.10l - l4.lT] and Corollaries I4.13H4.171 then we obtain criteria of reality (negativity, 
positivity) and simplicity for all zeros of a given polynomial. 

Now we give criteria of reality for all zeros of a given polynomial in terms of the determinants 6j{p), 
which are the leading principal minors of the matrix 'D2n{p) defined by (|4.14p (see e.g. [5^ 1102) ). 

Theorem 4.18. The polynomial p has exactly m (< n) distinct zeros, all of which are real, if and only if 

6iip)>0, S3ip)>0,..., 52m-i(p)>0. (4.25) 

Sj{p)^0, j>2m, (4.26) 
Theorem 4.19. Let the polynomial p have exactly m(< n) distinct zeros, all of which are real. If 

l^Y{l,S2ip),Siip),...,S2rnip)), 

then the number of distinct positive zeros of the polynomial p equals I (when p(0) ^0) or I — I (when 
p{0) ~ 0). The number of all positive zeros of p, counting multiplicities, is equal to Y^{ao,ai,a2, . . ■ ,a„). 

Proof. The theorem follows from Theorem 14. 8| formulae ()4.16p - (|4.17|) and Remark 14.121 □ 

From Theorems 14.191 and 14.181 and our previous results we obtain the following evident consequences: 

Corollary 4.20. The polynomial p has exactly m (< n) distinct zeros, all of which are positive, if and 
only if (j^25l) - (j4:26l) hold and 

(52^-2 (p)'52j (p) < 0, j = l,2, ...,m, {5q{p):=1). 

Corollary 4.21. The polynomial p has exactly m{< n) distinct zeros, all of which are positive, if and 
only if (|^25)) - (|4:26l) hold and 

a-j^ittj < 0, j = l,2,...,n. 

Corollary 4.22. The polynomial p has exactly m (< n) distinct zeros, all of which are negative, if and 
only if the equalities ()4.26|) hold and 

Slip) > 0, S2{p) > 0, . . . , 52,n-l{p) > 0, S2m{p) > 0. 

Corollary 4.23. Let all the coefficients of the polynomial p be positive. The polynomial p has exactly 
m(< n) distinct zeros, all of which are negative, if and only if (|4.25|) - (|4.26|) hold. 

For m = n, Theorems I4.181fl.l9l and Corollaries I4.201f4.23l become criteria of reality (negativity, posi- 
tivity) and simplicity for all zeros of a given polynomial. 

Corollary 14.221 with m = n is per se the Hurwitz stability criterion for polynomials whose zeros are 
real and simple; it is analogous to the standard Hurwitz criterion of polynomial stability |361 Chapter XV, 
Section 6] . In turn. Corollary 14.231 with m = n is an analogue of the Lienard and Chipart criterion j36l 
Chapter XV, Section 13]. Just as the criterion of Lienard and Chipart has a number of versions [361 
Chapter XV, Section 13], we can also obtain other versions of Corollarv 14.231 e-g-j as follows: 
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Corollary 4.24. Let all the coejficients of the polynomial p be positive. Then all roots of the polynomial 
p are negative and distinct if and only if the following inequalities hold: 

52{p) > 0, J4(P) > 0, . . . , 52n{p) > 0. (4.27) 

Proof. The theorem follows from Theorem 13.341 applied to the pair p,p'. □ 

Note that all criteria of Hurwitz stability require n inequalities on the coefhcients of a polynomial of 
degree n 136) . while the criteria of simplicity and negativity of zeros require 2n inequalities on the 
coefficients of a polynomial of degree n. 

The similarity between Hurwitz stable polynomials and polynomials with simple and negative zeros is 
also displayed by the following theorem, which is an analogue of the famous fact that the Hurwitz matrix 
of a Hurwitz stable polynomial is totally nonnegative [6l \57\ [46l [55] . 

Theorem 4.25. The polynomial p of degree n has only nonpositive zeros if and only if its matrix I?2n(p) 
defined in (|4.14p is totally nonnegative. 

Proof. If n = 0, then the assertion is evident. Let n > 1. From Theorems I4.101fl.lll and 13.41 it follows 
that the polynomial p has nonpositive zeros if and only if its logarithmic derivative is an i?-function of 
negative type with nonpositive poles and negative roots. Since I?2n(p) — ^2n(p,p'), the necessity direction 
of the theorem follows from Theorem 13.461 applied to the pair {p,p'). The sufficiency can be proved as in 
Theorem 13.441 using the factorization (|3.72|) where Tk{g) is totally nonnegative whenever g — gcd{p,p'). □ 

From Theorem 13.441 we derive the following corollary. 

Corollary 4.26. The polynomial p of degree n has only nonpositive zeros if and only if the infinite matrix 
H{p,p') defined in (I1.123P is totally nonnegative. 

Now, the following corollaries about Stieltjes continued fractions of logarithmic derivatives can be 
derived from Theorems 13. 4[ 14.101 and 14.111 and from Corollary 14.131 

Corollary 4.27. The polynomial p has exactly m (< n) distinct zeros, all of which are positive, if and 
only if its logarithmic derivative L has a Stieltjes continued fraction expansion (j4.10p " (|4.11l) where the 
inequalities (|4.23p hold and where 

C2j < 0, j = l,2,...,m. 

Corollary 4.28. The polynomial p has exactly m (< n) distinct zeros, all of which are negative, if and 
only if its logarithmic derivative L has a Stieltjes continued fraction expansion (|4?T0l) - (|4?TT|) where 

Ci > 0, i — 1,2, ... , 2m. 

Note that, in Corollaries I4.27ff4.281 all zeros of the polynomial p are simple if and only if m = n. 
At last, we present the following fact, which is a simple consequence of CoroUarv 13.521 

Theorem 4.29. The polynomial 

g(z) = ao + aiz + • • • + a„z", oq ^ 0, > (4.28) 
has all negative zeros if and only if the infinite matrix 





ai 
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3o3 


4a4 




[■■ 

















(4.29) 



is totally nonnegative. 
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Proof. In fact, consider the logarithmic derivative L of the polynomial g 

g'{z) fli + 2a2Z + • • • + nunz"^^ 



L{z) 



ao + aiz + a2Z^ + • • • + a„0'' 



z 



^2 



£2 

^3 



Suppose that the polynomial g has exactly m{< n) distinct zeros, all of which are negative. Then by 
Theorems [Ol 14.101 and 14.111 we have 



D,{L)>Q, j = l,2,...,m, 

(-iy5,(i) >0, j = l,2,...,m, 
D,{L) = D,{L)=Q, j>m. 



(4.30) 
(4.31) 
(4.32) 



According to Theorems 13.41 and 13.91 from (I4.30p - (|4.32p we obtain that the logarithmic derivative L of 
the polynomial g is an i?-function of negative type with all poles and zeros negative. Moreover, all common 
zeros of the numerator and the denominator of the function L are also negative. The converse statement 
is obviously true. 

If now we take h — g' in Corollary 13.521 then R = L and Hoc{g,g') — T^ocig)- So, the assertion of the 
theorem immediately follows from Corollarv 13.521 □ 

Remark 4.30. The necessary condition in this theorem is essentially known (see [28l Remark 3.23]) but 
the sufficient condition is most likely new. 

For the polynomial (|4.28|) with negative zeros, the total nonnegativity of the matrix (j4.29|) implies two 
properties of independent interest (cf. |28[ p. 66]): 



J'j+i 



a,_i 



>0, j = l,2,...,n-l 



(log-concavity) 



and 



2 3 + 1 1 
^ ] 3 



juj {j + l)aj+i 



> 0, J = 1, 2, . . . , n — 1 (weak Newton's inequalities) 



Remark 4.31. From Theorem 14.291 one can also prove Corollarv 13. 561 



At last, we present a couple of concrete examples. 

Example 4.32. Consider the following real polynomial 

f{z)=z^ + az + b, a, &eM. 

By our methods, we obtain a well-known fact 1641 that the polynomial f has only real zeros if and only if 
the coefficients a and b satisfy the ineaualit^^ 



4a-' 



2W < 0, 



At the same time, we show that the zeros of f cannot all be of the same sign. 
In fact, construct the matrix Ae{f) corresponding to the polynomial f 



and find its leading principal minors 6j{f) {j — 1,2, ... , 6); 



/3 





a 










1 





a 


b 











3 





a 











1 





a 


b 











3 





a 





^0 





1 





a 





<5i(/)=3; 62{f) = 



3 
1 



= 0; S^if)^ 



3 0a 
10a 
3 



= —6a; 



•^®We note that 4a^ + 276'^ is the discriminant of the polynomial /. 



(4.33) 



(4.34) 



74 



-4a2; 



-4a^ - 276^ 





'^6(/) = p6(/)| = 6<^5(/). 

According to Theorem \4-18\ the polynomial f has only real zeros if and only if the following inequalities 
hold: 

Siif)>0, <53(/)>0, (55(/)>0. (4.35) 

Moreover, the equality 63(f) ~ must imply 65(f) — 0. But it is easy to see that the case 63(f) = 65(f) = 
is possible if and only if a — b — 0. 

The first inequality (j4.35p holds automatically (see Remark \4-9\ l. The second inequality (j4.35p gives 
a < 0. From the third inequality it follows that the necessary and sufficient condition for the polynomial 
f to have only real zeros is the inequality (|4.33l) . This inequality also holds for a ^ b ~ and implies the 
inequality a < 0. 

If the polynomial f has only real zeros, then, according to Corollaries \4. 22 \4 . 2S\ all its zeros are negative 
or positive if all its coefficients are nonzero. But the coefficient of f(z) at z vanishes, therefore, f cannot 
have all zeros of the same sign. 

Thus, we proved that f has only real zeros if and only if a and b satisfy the inequality (j4.33p . We also 
proved that f cannot have only positive or only negative zeros for any real a and b. 

Example 4.33. From Theorem \4.18\ it follows that if the polynomial p defined by (j4.2ip is of degree n > 3 
with ai — 02 = and 7^ 0, then p cannot have only real zeros since we have 

S3(p) = 0, 

and 

65(p) ~ — OnaQflg < 0. 

Therefore, 

sign 65(p) = - signao = - sign(5i(p). 
This contradicts the inequalities (|4.25p . 
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